Morphological Characteristic Evaluation of Road Aggregates

Morphological Characteristic Evaluation of Road Aggregates: Comparison

Please note this is a comparison between Version 2 by Jue Li and Version 1 by Jue Li.

The sustainable performance of asphalt pavement depends on the quality and mix design of road aggregates. Identifying aggregate morphology and size is a prerequisite step for material design and numerical modeling of asphalt mixtures.

road aggregate
morphological feature
mechanical property
aggregate
morphology
asphalt mixture
shape
angularity
texture

1. Introduction

Asphalt mixtures are used in over 90% of highway construction in China. Currently, asphalt pavement generally has a short service life, and there is a gap in the design life. Aggregates account for more than 90% of the total mass of asphalt mixtures. Mineral aggregate, a vitally important road material, forms a good load-bearing capacity after gradation combination and sufficient compaction and can be used directly in the paving of the subgrade or subbase of road structures [1]. At the same time, multiphase composites are structurally characterized by a coarse aggregate skeleton and fine aggregate filling through binders such as asphalt and cement [2]. These materials have long dominated the pavement structure. Using aggregate properties to improve their mixes’ mechanical properties and service life has been an important area of research in road engineering.

2. Morphological Characteristic Evaluation

The engineering properties of road materials (e.g., asphalt mixes and cement concrete) are influenced by the morphological characteristics (shape, angle, and texture) of the aggregates [15]. Therefore, suitable evaluation methods and indicators can better quantify the morphological characteristics of aggregates. Specifically, the shape profile and angularity of coarse aggregates belong to the macroscopic category, while the surface texture belongs to the microscopic category [16]. The relationship among shape, angularity and surface texture is shown in Figure 1.

Figure 1. Relationship among shape, angularity and surface texture.

2.1. Shape

The shape of the aggregate can be divided into two dimensions and three dimensions, and the analysis of the aggregate shape can be divided into sphericity and shape factor, roundness, formation index, Fourier series, flat length ratio, aspect ratio, width ratio, symmetry, and three-dimensional shape factor [17]. As for the contour characteristics of particles, the concept of two-dimensional planes is mainly involved. With the development of digital image technology, the quantitative indexes of contour shape have progressed from shallow to deep, from the simple fine length to the development of relatively fine angular roundness. The fine length, also known as the axis ratio, is generally defined as the ratio of the long axis to the short axis of the equivalent ellipse of the planar projection of the aggregate.

The indicators characterizing the shape of the aggregates can be distinguished into 2D and 3D. Among them, most of the 2D metrics use length and area for shape characterization, such as form index (2D), roundness, the ratio of breadth to width, etc. It is easy to lose a lot of shape detail information. The aggregates’ shape, angularity and texture are effectively separated by the Fourier series method. Three-dimensional indexes include shape factor (SF), flat and elongated ratio (FER), and sphericity (ψ). These indexes can characterize the aggregate morphology more accurately. The shape index of aggregate with calculation equations and main characteristics is shown in Table 1.

Table 1. Shape index of aggregate with calculation equations and main characteristics.

Shape Indexes	Equations *	Characteristics	Literatures
Form index (2D)	$\sum_{θ = 0}^{360 - Δ θ} \frac{\|R_{θ + Δ θ} - R_{θ}\|}{R_{θ}}$	Quantifying the 2D image of fine aggregates of relative form and using incremental changes in the particle radius.	[18]
Aspect ratio (AR)	$\frac{D_{\max}}{D_{\min}}$	These parameters use the length and area method to characterize the shape of the aggregate; however, it loses a large amount of shape detail information	[19]
Rectangular degree (RD)	$\frac{A}{A_{M E R}}$

—the total number of CT images,

i

—the

i

th CT image.

p_{c o n v e x}

—perimeter of an equivalent ellipse;

p_{e l l i p s e}

—perimeter of the bounding polygon.

A_{p a r t i c l e}

—area of the real projection of aggregate;

A_{c o n v e x}

—area of convex aggregate’s projection.

P_{i}

—mass percentage composition of each grade of particle size;

I_{a i}

—angularity coefficient of each grade of particle size.

Liu et al. [47] reconstructed the microstructure model of the asphalt mixes used for uniaxial compression tests based on X-ray CT scans, keeping the original morphology of the aggregates. The 3D finite element model with different aggregate angles was established and simulated for uniaxial compression tests by artificially reducing the aggregate angles by keeping the original morphology unchanged. The effect of the aggregate angle on the mechanical response of the asphalt mixture was quantitatively analyzed. In order to obtain the probability distribution of shape features, GA values were proposed in the previous study to quantify the 2D morphology of agglomerate particles by AIMS tests. The GA values of the aggregates are shown in Figure 3.

Figure 3. Distribution of GA values for aggregates.

Zhu et al. [48] collected the aggregate characteristics by a 3D blue-ray scanner, and evaluated the effect of coarse aggregate size and angles on the volumetric properties of asphalt mixtures. Huang et al. [49] prepared mixes with similar aggregate gradation using five different angular levels of coarse-grained aggregates and three different asphalt binders and evaluated the coarse aggregate angles (CAA) by compacted voids in coarse aggregates (VCA) and triaxial shear tests.

In summary, researchers have proposed various evaluation indicators to further improve the road performance of road aggregates from different angles by studying the angularity of aggregates. Currently, many shape indices have been proposed by mathematical methods. However, most of the studies of these indices are based on graphical and statistical laws. Their correlation with the mechanical properties of aggregates lacks research. At the same time, the applicability of these indices for the prediction of the structural properties of the mixture still lacks a uniform classification standard.

2.3. Texture

The aggregate texture is mainly distinguished by the smoothness and roughness of aggregate particles. Miller et al. [50] proposed an analytical method to describe the surface texture of asphalt pavements using laser profiles to estimate friction characteristics. Al-Assi et al. [51] used Close-Range Photogrammetry (CRP) to measure the changes in micro texture on asphalt pavements. Khasawneh et al. [52] proposed an equation based on nonlinear regression to calculate the texture of asphalt pavements and the steady state’s texture features. Du et al. [53] utilized a 2D wavelet method to describe the texture obtained by a 3D laser scanner on asphalt pavement. Hu et al. [54] conducted an analysis of the influence of texture on the skid resistance of asphalt pavements by using the collected pavement surface texture data. The texture index of aggregate with calculation equations and main characteristics is shown in Table 3.

Table 3. Texture index of aggregate with calculation equations and main characteristics.

Texture Indexes	Equations *	Characteristics	Literatures
Texture factor (TF)	$\sum_{p = 1}^{N} \sum_{q = 1}^{N} [{(\frac{a (p, q)}{a_{0}})}^{2} + {(\frac{b (p, q)}{a_{0}})}^{2}] - A F$	Fast measurement of aggregated morphological features based on 2D fast Fourier transform.	[55]
Texture index (TI)	$\frac{1}{3 N} {\sum_{i = 1}^{3} \sum_{j = 1}^{N} (D_{i, j} (x, y))}^{2}$	Quantization of textures using wavelet methods.	[55]
Erosion dilation area ratio (EDR)	$\frac{A_{1} - A_{2}}{A_{1}} \times 100 %$ $\frac{A}{A_{M E R}}$	[20]
$\frac{A_{1} - A_{2}}{A_{1}} \times 100 %$	$\frac{A_{1} - A_{2}}{A_{1}} \times 100 %$	Changes in area after erosion and dilatation cycles are directly related to aggregate texture.	[56]	Axial coefficient (AC)	$\frac{L_{\max}}{L_{\min}}$ Higher values of GA indicate more angular aggregate particles.
Comprehensive erosion dilation area ratio (CEDR)	$\frac{\sum_{i}^{n} (a_{i} \times E D R_{i})}{\sum_{i}^{n} a_{i}}$ $\frac{L_{\max}}{L_{\min}}$	[	43]
$\frac{\sum_{i}^{n} (a_{i} \times E D R_{i})}{\sum_{i}^{n} a_{i}}$	$\frac{\sum_{i}^{n} (a_{i} \times E D R_{i})}{\sum_{i}^{n} a_{i}}$	[21]
$\frac{\sum_{i}^{n} (a_{i} \times E D R_{i})}{\sum_{i}^{n} a_{i}}$	$\frac{\sum_{i}^{n} (a_{i} \times E D R_{i})}{\sum_{i}^{n} a_{i}}$	Aggregate gradation is considered.	[57	Eccentricity (E)	$\frac{c}{a}$	Eccentricity could characterize the shape of the aggregate particles and reflect the elliptical oblate degree.	[19]
]	Angularity index (AI) (radius)	$\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$
Surface texture	$\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$ $\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ $\sum_{θ = 0}^{θ = 360 - Δ θ} \frac{\|R_{p θ} - R_{E E θ}\|}{R_{E E θ}} θ$	As the defining equation of angularity.	$\sum_{m = n_{2} + 1}^{m = n_{3}} \sqrt{a_{m}^{2} + b_{m}^{2}}$ [27]
The shape, angularity and texture of the aggregates are considered separately.	[	38	]	Roundness	$\frac{4 π A}{L^{2}}$	Roundness was the inverse of the form factor. In early imaging techniques, it was used to calculate 2D features of particle shape.	3D angularity (3DA)	$\frac{1}{3} \frac{\sum_{i = 1}^{n_{t}} \frac{P_{t i} \cdot A_{t i}}{P_{t e}}}{\sum_{i = 1}^{n_{t}} A_{t i}} + \frac{\sum_{i = 1}^{n_{r}} \frac{P_{r i} \cdot A_{r i}}{P_{r e}}}{\sum_{i = 1}^{n_{r}} A_{r i}} + \frac{\sum_{i = 1}^{n_{f}} \frac{P_{f i} \cdot A_{f i}}{P_{f e}}}{\sum_{i = 1}^{n_{f}} A_{f i}}$
Fourier series analysis method	$\frac{1}{3} \frac{\sum_{i = 1}^{n_{t}} \frac{P_{t i} \cdot A_{t i}}{P_{t e}}}{\sum_{i = 1}^{n_{t}} A_{t i}} + \frac{\sum_{i = 1}^{n_{r}} \frac{P_{r i} \cdot A_{r i}}{P_{r e}}}{\sum_{i = 1}^{n_{r}} A_{r i}} + \frac{\sum_{i = 1}^{n_{f}} \frac{P_{f i} \cdot A_{f i}}{P_{f e}}}{\sum_{i = 1}^{n_{f}} A_{f i}}$	$\frac{1}{2} \sum_{m = n_{2} + 1}^{m = \infty} [{(\frac{a_{m}}{a_{0}})}^{2} + {(\frac{b_{m}}{a_{o}})}^{2}]$ $\frac{1}{3} \frac{\sum_{i = 1}^{n_{t}} \frac{P_{t i} \cdot A_{t i}}{P_{t e}}}{\sum_{i = 1}^{n_{t}} A_{t i}} + \frac{\sum_{i = 1}^{n_{r}} \frac{P_{r i} \cdot A_{r i}}{P_{r e}}}{\sum_{i = 1}^{n_{r}} A_{r i}} + \frac{\sum_{i = 1}^{n_{f}} \frac{P_{f i} \cdot A_{f i}}{P_{f e}}}{\sum_{i = 1}^{n_{f}} A_{f i}}$ $\frac{1}{2} \sum_{m = n_{2} + 1}^{m = \infty} [{(\frac{a_{m}}{a_{0}})}^{2} + {(\frac{b_{m}}{a_{o}})}^{2}]$ $\frac{1}{3} \frac{\sum_{i = 1}^{n_{t}} \frac{P_{t i} \cdot A_{t i}}{P_{t e}}}{\sum_{i = 1}^{n_{t}} A_{t i}} + \frac{\sum_{i = 1}^{n_{r}} \frac{P_{r i} \cdot A_{r i}}{P_{r e}}}{\sum_{i = 1}^{n_{r}} A_{r i}} + \frac{\sum_{i = 1}^{n_{f}} \frac{P_{f i} \cdot A_{f i}}{P_{f e}}}{\sum_{i = 1}^{n_{f}} A_{f i}}$	[22]
$\frac{1}{3} \frac{\sum_{i = 1}^{n_{t}} \frac{P_{t i} \cdot A_{t i}}{P_{t e}}}{\sum_{i = 1}^{n_{t}} A_{t i}} + \frac{\sum_{i = 1}^{n_{r}} \frac{P_{r i} \cdot A_{r i}}{P_{r e}}}{\sum_{i = 1}^{n_{r}} A_{r i}} + \frac{\sum_{i = 1}^{n_{f}} \frac{P_{f i} \cdot A_{f i}}{P_{f e}}}{\sum_{i = 1}^{n_{f}} A_{f i}}$	3D parameters for more accurate characterization of aggregate morphology.	[	43]	Form index using Fourier series (FRFORM)
Calculation of the surface fabric of aggregates used in PIAS.	]	$R (θ) = a_{0} + \sum_{m = 1}^{\infty} [a_{m} \cos (m θ) + b_{m} \sin (m θ)]$	Angularity parameter $R (θ) = a_{0} + \sum_{m = 1}^{\infty} [a_{m} \cos (m θ) + b_{m} \sin (m θ)]$	${(\frac{p_{c o n v e x}}{p_{e l l i p s e}})}^{2}$
Direct measurement of aggregate dimensions	${(\frac{P}{P_{c o n v e x}})}^{2}$ ${(\frac{p_{c o n v e x}}{p_{e l l i p s e}})}^{2}$ ${(\frac{P}{P_{c o n v e x}})}^{2}$ ${(\frac{p_{c o n v e x}}{p_{e l l i p s e}})}^{2}$ ${(\frac{P}{P_{c o n v e x}})}^{2}$ ${(\frac{p_{c o n v e x}}{p_{e l l i p s e}})}^{2}$ ${(\frac{P}{P_{c o n v e x}})}^{2}$ ${(\frac{p_{c o n v e x}}{p_{e l l i p s e}})}^{2}$ ${(\frac{P}{P_{c o n v e x}})}^{2}$ ${(\frac{p_{c o n v e x}}{p_{e l l i p s e}})}^{2}$ ${(\frac{P}{P_{c o n v e x}})}^{2}$ ${(\frac{p_{c o n v e x}}{p_{e l l i p s e}})}^{2}$ ${(\frac{P}{P_{c o n v e x}})}^{2}$ ${(\frac{p_{c o n v e x}}{p_{e l l i p s e}})}^{2}$ ${(\frac{P}{P_{c o n v e x}})}^{2}$ $R (θ) = a_{0} + \sum_{m = 1}^{\infty} [a_{m} \cos (m θ) + b_{m} \sin (m θ)]$	Fourier series can be used to analyze the form, angularity, and texture of aggregate shape.	[23]
For evaluating textures, the texture parameters are related to the perimeter.	[		44	]	Form (shape) index(Fourier series)	$\frac{1}{2} \sum_{m = 1}^{m = n_{1}} [{(\frac{a_{m}}{a_{0}})}^{2} + {(\frac{b_{m}}{a_{0}})}^{2}]$
Three-dimensional Texture ( ${MT}_{3 D} SHF$ )	$\frac{1}{2} \sum_{m = 1}^{m = n_{1}} [{(\frac{a_{m}}{a_{0}})}^{2} + {(\frac{b_{m}}{a_{0}})}^{2}]$	[23]
Ratio of Breadth to Width	$\frac{\min (x_{c})}{\max (x_{F e})}$	Breadth to width ratio can be used to describe the form of aggregate particles.	[24]
Symmetry	$\frac{1}{2} [1 + \min (\frac{r_{1}}{r_{2}})]$	Symmetry is a term that some imaging systems use to describe aggregate form.	[24]
Shape Factor (SF)	$\frac{D_{S}}{\sqrt{D_{L} \times D_{i}}}$	Shape Factor was a typical index in a system that represented in terms of three dimensions (longest, intermediate, and shortest dimensions).	[25]
Flat and elongated ratio (FER)	$\frac{D_{L}}{D_{S}}$	FER represents the ratio between the longest dimension and the shortest dimension of a particle.	[26]
Sphericity ( $ψ$ )	$\sqrt[3]{\frac{D_{S} \times D_{i}}{D_{L}^{2}}}$	Sphericity is a 3D parameter, which can characterize the aggregate morphology more accurately.	[27]

θ

—directional angle;

Δ θ

—incremental difference in the angles;

R_{θ}

—radius of the particle at angle

θ

;

R_{θ + Δ θ}

—radius of the particle at angle

θ + Δ θ

D_{m a x}

—the maximum length of the minimum external rectangle of the aggregate profile;

D_{m i n}

—minimum size of the minimum external rectangle.

A

—Area of the aggregate profile;

A_{M E R}

—area of minimum circumscribed rectangle.

L_{m a x}

—major axis length of equivalent ellipse;

L_{m i n}

—minor axis length of equivalent ellipse.

c

—focal length of ellipse with same second moment;

a

—major axis of the ellipse.

L

—Circumference of the aggregate profile.

m

—frequency;

a_{0}

—average radius;

a_{m}, b_{m}

—Fourier coefficients

n_{1}

—threshold frequencies separating shape and angularity.

x_{c}

—maximum chord;

x_{F e}

—the Feret diameter.

r_{1}

r_{2}

—the distances of the center of gravity to the edge in a given direction.

D_{L}

—length, longest dimension;

D_{i}

—width, intermediate dimension, the maximum dimension in the plane perpendicular to the length;

D_{S}

—thickness, shortest dimension.

Mora et al. [28] proposed a method for measuring the sphericity, shape factor and convexity of concrete aggregates, which estimated the thickness and volume of the particles, measured the thickness-dependent shape parameters, and evaluated the weighted average of individual particle shape parameters in aggregate specimens. Al-Rousan et al. [29] evaluated the accuracy of the analytical methods used in currently available imaging systems. Komba et al. [30] reconstructed a three-dimensional model of aggregate particles after obtaining data from three-dimensional laser scans of various types of aggregates. This laser scanning technique can better quantify the morphological characteristics of aggregates by calculating the sphericity from surface area and volume, the sphericity from three orthogonal dimensions, and the flat-length ratio from the longest and shortest dimensions of the aggregate particles.

The brittleness index is an important shape parameter in aggregates widely used in infrastructure construction, such as roads, and Anochie-Boateng et al. [31] developed a new equation for the brittleness index and verified its reliability. Laser scanning technology provided more accurate evaluation results and improved the selection of construction aggregates compared to traditional evaluation methods. Pan [32] and Tutumluer [33] evaluated the surface characteristics of crushed and uncrushed aggregates in asphalt mixtures using the 3D laser scanning technique, which had a good correlation with the experimental results. Ge et al. [34] proposed the flatness elongation index, sphericity index and specific surface area to characterize various aspects of the morphological characteristics of coarse aggregates and derived that small-size aggregates have a larger proportion of elongated particle number.

From the above research, it can be seen that the 3D scanning technique can well-characterize the aggregate’s shape. However, previous studies using 3D scanning technology have mainly focused on the shape, volume, and surface properties of aggregates from one or different quarries, with less research on the morphological decay pattern of aggregates during application.

2.2. Angularity

Aggregate angularity essentially describes the degree of sharpness of the angles. Corner roundness was proposed by Wadell [35] and others to describe the sharpness of the combination of corners and edges of a particle and was generally considered to be related to the radius of curvature or geometric fractal of the particle profile at the sharp end, while for quantitative evaluation, most researches have used the method of fitting similar polygons to obtain the angularity parameter of an aggregate. Fractals were first created by Mandelbrot [36] and have been widely used in various disciplines through continuous development. Fractals are used as a theoretical tool to describe and evaluate some characteristics such as hierarchical continuity, irregularity or randomness [37], and asphalt mixtures exhibit statistical self-similarity in volume and properties, which is suitable for quantitative analysis with fractal theory.

Wang et al. [38] proposed a unified Fourier morphological analysis method to quantify the shape characteristics of aggregates. Zhang et al. [39] proposed an index to characterize the combined effect of rib angle and surface texture of coarse aggregates, i.e., AT, and investigated the statistical distribution pattern of the AT index for different grain sizes of basalt aggregates and their composite aggregates. Rajan et al. [40] studied the comparison of rib angles of different grain sizes of aggregates under different crusher production. The angularity index of aggregate with calculation equations and main characteristics is shown in Table 2.

Table 2. Angularity index of aggregate with calculation equations and main characteristics.

Angularity Indexes	Equations *	Characteristics	Literatures
Angularity	$\frac{1}{2} \sum_{m = n_{1} + 1}^{m = n_{2}} [{(\frac{a_{m}}{a_{0}})}^{2} + {(\frac{b_{m}}{a_{0}})}^{2}]$	Fourier series analysis was applied to measure the angularity of aggregate in PIAS.	[38]
Angularity factor (AF)	$A F = \sum_{p = 1}^{N_{A}} \sum_{q = 1}^{N_{A}} [{(\frac{a (p, q)}{a_{0}})}^{2} + {(\frac{b (p, q)}{a_{0}})}^{2}]$	The fast Fourier transform method could be applied to calculate and analyze the angularity factor (AF) of aggregate in FTI.	[41]
Surface erosion-dilation method	$\frac{A_{1} - A_{2}}{A_{1}} \times 100 %$	As a popular image processing technique, erosion dilation was utilized to analyze the angularity of aggregate particles in the UIAIA.	[42]
Gradient angularity (GA)	$\frac{1}{\frac{n}{3} - 1} \sum_{i = 1}^{n - 3} \|θ_{i} - θ_{i + 3}\|$
All of them optimize the profile of the particles of the aggregate but cannot retain the original morphological characteristics of the aggregate.
	[	44	]
	Convexity	$\sqrt{\frac{A_{p a r t i c l e}}{A_{c o n v e x}}}$	[45]
Average angularitycoefficient	$I = \frac{\sum_{i = 1}^{n} P_{i} \times I_{a i}}{100}$	[46]

n_{1}

n_{2}

—threshold frequencies separating shape, angularity, and texture.

N_{A}

—threshold value;

N

—size of

z (x, y)

matrix;

a, b

—the real and imaginary roots of FFT2 coefficients;

a_{0}

—the average value of

z (x, y)

A_{1}

—initial area (in pixel) of 2D projection of particle in the image;

A_{2}

—area (in pixel) of the particle after performing a sequence of ‘n’ cycles of erosion—dilation.

θ

—directional angle;

i

—the

i

th point on the boundary of aggregate;

n

—total number of the points on the boundary of aggregate.

R_{p θ}

—radius of the particle at a directional angle

θ

;

R_{E E θ}

—radius of an equivalent ellipse at the same

θ

P_{t i}, P_{r i}, P_{f i}

—the parameters of grain particles;

P_{t e}, P_{r e}, P_{f e}

—the equivalent ellipse perimeters of grain particles;

A_{t i}, A_{r i}, A_{f i}

—the areas of grain particles;

n_{t}, n_{r}, n_{f}

[
23
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$	${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$ ${MT}_{3 D N} = \frac{M T_{3 D}}{a_{00}}$	3D textures can be calculated.Can differentiate between aggregates with similar sizes.Used to compare aggregates of different sizes and volumes.	[58]

N

—the level of decomposition;

i

—a value of 1, 2, or 3 for the three detailed images of texture;

j

—the wavelet coefficient index.

A_{1}

—the 2D profile area of the aggregate;

A_{2}

—the 2D profile area of the aggregate after the erosion and expansion.

a_{i}

—the proportion of the

i

th grade aggregate in the gradation;

E D R_{i}

—erosion expansion area ratio of the

i

th grade of aggregate in the sample.

P

—the perimeter of a particle measured on a black and white image;

P_{c o n v e x}

—the perimeter of a bounding polygon.

M T_{3 D}

—3D microtexture.

An asphalt mixture is a typical non-homogeneous composite material consisting of irregularly shaped and distributed aggregates, asphalt binder and voids [59,60]. The aggregate morphological characteristics significantly influence the properties of asphalt mixtures [61,62], which can be found in existing studies. Improving the performance of asphalt mixes based on morphological characteristics of aggregates accounting for 90% of the total mass of asphalt mixes is an important research direction. Traditional methods of testing aggregate morphological characteristics are time-consuming and rely heavily on subjective judgments, resulting in imprecise experimental results. In most of the tests, it is difficult for the researchers to control the aggregate morphological characteristics variables [63]. Aragão et al. [64] mixed aggregates from different sources and tested them using AIMS to obtain aggregates with different angles and textures; nevertheless, it was impossible to distinguish which morphological characteristics were responsible for the effect. Wang et al. [65] employed a modified Los Angeles (LA) wear test to obtain aggregates with different angles, calculated the fractal dimension and changed the size and surface texture of the aggregates. Puzzo et al. [66] proposed a 3D model evaluation method for pavement texture based on digital image processing techniques, which used photographs to generate a Digital Surface Model (DSM) to calculate the digital Mean Texture Depth (MTD) and calculate and analyze other texture parameters for contour extraction.

Cui [55] evaluated the effect of aggregate morphology variation on the performance of asphalt pavements and investigated the asphalt coverage ratio and the mechanical properties of the mixture for different types of aggregates to quantify the aggregate morphological characteristics by AIMS. An accurate definition of the aggregate morphological characteristics is important for the initial design phase of pavements, extending the service life of roads and reducing the use of natural aggregates and the generation of solid waste.

Inadequate adhesion between asphalt and aggregate can cause severe distress in pavement structure and service life [67], and the conventional approach, which does not consider active adhesion mechanisms, is hardly convincing. Cui et al. [68] investigated and designed an active adhesion-based test method to compensate for the adhesion deficiency and used digital image processing techniques to quantify the adhesive properties. The results showed that digital image processing could measure the asphalt coating ratio accurately and effectively.

In summary, the texture is a micro-morphological characteristic of the aggregate. The greater the complexity of the surface texture, the greater the internal friction angle that can be generated when the aggregates are embedded and locked with each other after crushing. Aggregate texture has fewer means of characterization relative to shape and angularity. At present, there is no excellent way to judge the state of the textural distribution of asphalt pavement surface texture. Further development of evaluation metrics is needed to better predict pavement friction performance.

[
23
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$
${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$	${MT}_{3 D} = \sum_{n = 26}^{n_{\max}} \sum_{m = - 1}^{n} \|a_{n m}\|$ ${MT}_{3 D N} = \frac{M T_{3 D}}{a_{00}}$	3D textures can be calculated.Can differentiate between aggregates with similar sizes.Used to compare aggregates of different sizes and volumes.	[58]