# Automatic Measurement Method of Slump Flow of Cement Paste Based on Dynamic Image Processing Technology

Authors: Shiqin He; Pengfei Gao; Chunyue Wang; Ding Nie; Hui Wang
Journal: Advances in Materials Science and Engineering
 (2022)

---

## Abstract

The slump flow and flow time T200 are the main parameters to evaluate the rheological properties of cement paste. Generally, the manual measurement method is used and it is difficult to get T200 accurately. In this paper, an automatic measurement method of slump flow based on dynamic image processing technology is studied. Firstly, the cement paste flowability test video is divided into frames every 100 ms; secondly, based on the principle of perspective transformation, an automatic point selection method is proposed to correct the frame image to solve the problem of shooting angle; thirdly, the watershed algorithm based on distance transformation is used to segment the boundary of the cement paste area, and the slump flow is calculated according to the polar coordinates of the boundary; lastly, the relationship curve between slump flow and time is automatically output, and the slump flow at any time is obtained. The method is verified by nine groups of cement paste flowability tests. The results show that this method can fastly and accurately obtain the slump flow and flow time T200. Compared with the traditional manual measurement method, the measurement error of slump flow is less than 2%, and the measurement error of flow time T200 is less than 8%. The method enables real-time dynamic measurement of the slump flow.


---

## Body

## 1. Introduction

It is generally accepted that fresh cement-based materials such as concrete, cement mortar, and cement paste are Bingham fluids. Yield stress and plastic viscosity are the two basic physical parameters characterizing their rheological behavior [1]. The rheology of cement paste is an important parameter for evaluating the workability of self-compacting concrete. The two rheological parameters of cement paste can be calculated from the slump flow and the flow time T200, which are obtained from the cement paste flowability test [2, 3]. The slump flow is usually measured manually, while T200 is difficult to be obtained directly by the manual measurement and needs to read the cement paste flowability test video, which is a tedious operation. With the continuous upgrading of information technology and artificial intelligence industries, the automation of image processing has gradually become the future development trend. In image processing, orthogonal polynomial transform plays an important role, which is widely used in face recognition, target classification, template matching, and image edge detection [4, 5]. At present, there are few studies on the rheological process of cement paste by using image processing technology.Tregger et al. [6] used MATLAB software to take the position with the maximum image brightness gradient as the boundary of the cement paste and calculate the slump flow by the boundary coordinates information. This method strictly requires the camera imaging surface to be parallel to the glass plate, and T200 cannot be obtained accurately. Shen [7] took photos with a smartphone and estimated the slump flow by the proportion of concrete and a predefined sign in the image. This method requires artificially specifying concrete boundary points on the mobile phone, and the object needs to be photographed squarely. Li et al. [8] first calculated the spatial location coordinates of the identification image using the Vuforia Technology and took the distance between the two endpoints as the slump flow. This method solves the problem of shooting angle, but it still requires artificially selected boundary points. Based on the Vuforia package, Nie [9] first corrected the image according to the coordinate information of the identification image and took the position with the maximum image brightness gradient vector as the boundary of the cement paste. This method can accurately obtain the slump flow, but there is an obvious error in the measurement of T200. The improvement is difficult because the internal parameters of Vuforia cannot be adjusted.In fact, the measurement of the flow time T200 and the slump flow boils down to finding the diameter of the cement paste at a specific moment. In this paper, based on Python language and OpenCV library [10], an automated detection technology of slump flow is studied. Through the framing processing, image correction, boundary segmentation, and other steps of the cement paste flowability test video, the relationship curve between slump flow and time can be output automatically, and the slump flow corresponding to any moment can be obtained; the automatic detection process is achieved. Moreover, the measurement accuracy is greatly improved. The image processing technology is applied to the field of cement-based material workability testing to improve efficiency and lay the foundation for smart production of concrete.


## 2. Processing Method for Images of Cement Paste

### 2.1. Method of Image Correction Automatically

It is difficult to achieve image alignment when shooting the cement paste flowability test video, so the images need to be corrected first. In this paper, the perspective transformation [11] method is used to correct the glass plate, as shown in(1)x′,y′,w′=u,v,wa11a12a13a21a22a23a31a32a33,where x′,y′,w′ is the transformed coordinate, u,v,w is the coordinate before the transformation, a11a12a13a21a22a23a31a32a33 is the transformation matrix, and a33 is a constant with the value 1. Since the image being processed is two-dimensional, w is constant at 1. As shown in the following equation,(2)x=x′w′,y=y′w′.The two-dimensional image transformation equation can be obtained. As shown in the following equation,(3)x=a11u+a21v+a31a13u+a23v+1,y=a12u+a22v+a32a13u+a23v+1.Therefore, the transformation matrix can be solved only by giving the coordinates of the original image’s four corner pointsu,v and the transformed image’s target coordinates. The critical problem is how to determine the coordinates of the four corner points of the image before the transformation. In this paper, two methods are proposed to determine the coordinates of four corner points before transformation: manual point selection and automatic point selection methods.


#### 2.1.1. Manual Point Selection Method

This method calls the mouse event in OpenCV, using the function setMouseCallback(). By tapping the four corner points of the glass plate in the image, their coordinate values can be obtained [12]. The order of the selected points is recommended to be A, B, C, and D. The correction results of the manual point selection method are shown in Figure 1.Figure 1
Correction results of manual point selection method: (a) original image; (b) corrected image.
(a)(b)


#### 2.1.2. Automatic Point Selection Method

Based on the analysis of the characteristics of the corner coordinates of the original image, a method for automatically obtaining the four corner coordinates of the glass plate before transformation is proposed. As shown in Figure2, M and K are points with the smallest and largest X coordinate values; N and L are points with the smallest and largest Y coordinate values, respectively. The NumPy library can obtain these four corner coordinates in Python [13]. Just rotate the glass plate counterclockwise at a certain angle to make its corner coordinates comply with Figure 2. The result of automatic correction of the glass plate using the perspective transformation function is shown in Figure 3. The result of the automatic correction of the glass plate using the perspective transformation function is shown in Figure 3.Figure 2
Diagram of four corner points on the coordinate axis.Figure 3
Correction results of automatic point selection method: (a) original image; (b) corrected image.
(a)(b)


### 2.2. Image Boundary Segmentation

After image correction, the cement paste area needs to be segmented from the glass plate. The boundary of the cement paste can be found to calculate the slump flow, that is, the maximum diameter of the cement paste area. In this paper, three algorithms are compared as follows.


#### 2.2.1. Watershed Algorithm Based on Distance Transformation

The foreground and background are labeled with different colours (or grayscale values), and the areas of uncertainty are labeled with 0. Then, the watershed algorithm [14] is implemented to update the labels, and the boundary information of the target is obtained when the foreground and background meet [15]. In this paper, the area of cement paste is the foreground, and the area other than the cement paste is the background. The main steps are as follows:(1)
Binarization [16], preset threshold Q, and pixel points with grayscale values less than Q are set as black, and those greater than or equal to Q are set as white.(2)
Morphological operations [17] remove all white noise in the image.(3)
Distance transformation, generates a distance image that describes the distance between pixel points in the grayscale image and the nearest background pixel.(4)
Label creation, the areas of uncertainty are labeled with 0, and the background is labeled with a positive integer above 1.(5)
Watershed operation, the foreground and background are distinguished by updating the label.(6)
Foreground segmentation, the pixel positions labeled with 1 (background area) are synchronized to the original image.The results of the implementation of each step are shown in Figure4.Figure 4
The segmentation results of the watershed algorithm.


#### 2.2.2. HSV Colour Segmentation Algorithm

Compared with the RGB model, the HSV model is easier to track objects of a specific colour [18]. The H-value indicates the hue, the S-value indicates the saturation, and the V-value indicates the luminance. The colour of the cement paste area is relatively uniform and obviously different from other areas. Therefore, a range of pixel values can be obtained by extracting the HSV value of the cement paste area, the pixel value located in this range is set to 255 (white), and the pixel value below or above this range is set to 0 (black). The results are shown in Figure 5.Figure 5
The segmentation results of HSV.


#### 2.2.3. GrabCut Foreground Segmentation Algorithm

The GrabCut algorithm is based on colour features and edge features of the image [19]. The premise of image segmentation using the GrabCut algorithm is that there is a significant colour difference between foreground and background. The implementation process is to select the foreground area (cement paste) with the mouse; through the selected foreground information, obtain the remaining foreground and background; get the segmentation result after several iterations of calculation. The segmentation results may be different for different ways of foreground selection. As shown in Figure 6, the selection of foreground 2 cannot get a satisfactory segmentation result.Figure 6
The segmentation results of GrabCut algorithm: (a) foreground 1, (b) the segmentation result of foreground 1, (c) foreground 2, and (d) the segmentation result of foreground 2.
(a)(b)(c)(d)


### 2.3. Calculation of Slump Flow

In this paper, the watershed algorithm based on distance transformation is selected to segment the boundary of the corrected area of cement paste. The extreme points of the boundary are obtained by using the NumPy library [20]. As shown in Figure 7, the maximum diameter of the cement paste area can be calculated according to(4)D=P2−P1+P4−P32,where (P2 − P1) is the distance between points P1 and P2, and (P4 − P3) is the distance between points P3 and P4. After the video is frame-separated, the above operation is performed for all images in turn. The visualization tool Matplotlib is used to draw the relationship curve between slump flow and time, as shown in Figure 8.Figure 7
The boundary and extremum of cement paste.Figure 8
The SF-T curve.


## 2.1. Method of Image Correction Automatically

It is difficult to achieve image alignment when shooting the cement paste flowability test video, so the images need to be corrected first. In this paper, the perspective transformation [11] method is used to correct the glass plate, as shown in(1)x′,y′,w′=u,v,wa11a12a13a21a22a23a31a32a33,where x′,y′,w′ is the transformed coordinate, u,v,w is the coordinate before the transformation, a11a12a13a21a22a23a31a32a33 is the transformation matrix, and a33 is a constant with the value 1. Since the image being processed is two-dimensional, w is constant at 1. As shown in the following equation,(2)x=x′w′,y=y′w′.The two-dimensional image transformation equation can be obtained. As shown in the following equation,(3)x=a11u+a21v+a31a13u+a23v+1,y=a12u+a22v+a32a13u+a23v+1.Therefore, the transformation matrix can be solved only by giving the coordinates of the original image’s four corner pointsu,v and the transformed image’s target coordinates. The critical problem is how to determine the coordinates of the four corner points of the image before the transformation. In this paper, two methods are proposed to determine the coordinates of four corner points before transformation: manual point selection and automatic point selection methods.


### 2.1.1. Manual Point Selection Method

This method calls the mouse event in OpenCV, using the function setMouseCallback(). By tapping the four corner points of the glass plate in the image, their coordinate values can be obtained [12]. The order of the selected points is recommended to be A, B, C, and D. The correction results of the manual point selection method are shown in Figure 1.Figure 1
Correction results of manual point selection method: (a) original image; (b) corrected image.
(a)(b)


### 2.1.2. Automatic Point Selection Method

Based on the analysis of the characteristics of the corner coordinates of the original image, a method for automatically obtaining the four corner coordinates of the glass plate before transformation is proposed. As shown in Figure2, M and K are points with the smallest and largest X coordinate values; N and L are points with the smallest and largest Y coordinate values, respectively. The NumPy library can obtain these four corner coordinates in Python [13]. Just rotate the glass plate counterclockwise at a certain angle to make its corner coordinates comply with Figure 2. The result of automatic correction of the glass plate using the perspective transformation function is shown in Figure 3. The result of the automatic correction of the glass plate using the perspective transformation function is shown in Figure 3.Figure 2
Diagram of four corner points on the coordinate axis.Figure 3
Correction results of automatic point selection method: (a) original image; (b) corrected image.
(a)(b)


## 2.1.1. Manual Point Selection Method

This method calls the mouse event in OpenCV, using the function setMouseCallback(). By tapping the four corner points of the glass plate in the image, their coordinate values can be obtained [12]. The order of the selected points is recommended to be A, B, C, and D. The correction results of the manual point selection method are shown in Figure 1.Figure 1
Correction results of manual point selection method: (a) original image; (b) corrected image.
(a)(b)


## 2.1.2. Automatic Point Selection Method

Based on the analysis of the characteristics of the corner coordinates of the original image, a method for automatically obtaining the four corner coordinates of the glass plate before transformation is proposed. As shown in Figure2, M and K are points with the smallest and largest X coordinate values; N and L are points with the smallest and largest Y coordinate values, respectively. The NumPy library can obtain these four corner coordinates in Python [13]. Just rotate the glass plate counterclockwise at a certain angle to make its corner coordinates comply with Figure 2. The result of automatic correction of the glass plate using the perspective transformation function is shown in Figure 3. The result of the automatic correction of the glass plate using the perspective transformation function is shown in Figure 3.Figure 2
Diagram of four corner points on the coordinate axis.Figure 3
Correction results of automatic point selection method: (a) original image; (b) corrected image.
(a)(b)


## 2.2. Image Boundary Segmentation

After image correction, the cement paste area needs to be segmented from the glass plate. The boundary of the cement paste can be found to calculate the slump flow, that is, the maximum diameter of the cement paste area. In this paper, three algorithms are compared as follows.


### 2.2.1. Watershed Algorithm Based on Distance Transformation

The foreground and background are labeled with different colours (or grayscale values), and the areas of uncertainty are labeled with 0. Then, the watershed algorithm [14] is implemented to update the labels, and the boundary information of the target is obtained when the foreground and background meet [15]. In this paper, the area of cement paste is the foreground, and the area other than the cement paste is the background. The main steps are as follows:(1)
Binarization [16], preset threshold Q, and pixel points with grayscale values less than Q are set as black, and those greater than or equal to Q are set as white.(2)
Morphological operations [17] remove all white noise in the image.(3)
Distance transformation, generates a distance image that describes the distance between pixel points in the grayscale image and the nearest background pixel.(4)
Label creation, the areas of uncertainty are labeled with 0, and the background is labeled with a positive integer above 1.(5)
Watershed operation, the foreground and background are distinguished by updating the label.(6)
Foreground segmentation, the pixel positions labeled with 1 (background area) are synchronized to the original image.The results of the implementation of each step are shown in Figure4.Figure 4
The segmentation results of the watershed algorithm.


### 2.2.2. HSV Colour Segmentation Algorithm

Compared with the RGB model, the HSV model is easier to track objects of a specific colour [18]. The H-value indicates the hue, the S-value indicates the saturation, and the V-value indicates the luminance. The colour of the cement paste area is relatively uniform and obviously different from other areas. Therefore, a range of pixel values can be obtained by extracting the HSV value of the cement paste area, the pixel value located in this range is set to 255 (white), and the pixel value below or above this range is set to 0 (black). The results are shown in Figure 5.Figure 5
The segmentation results of HSV.


### 2.2.3. GrabCut Foreground Segmentation Algorithm

The GrabCut algorithm is based on colour features and edge features of the image [19]. The premise of image segmentation using the GrabCut algorithm is that there is a significant colour difference between foreground and background. The implementation process is to select the foreground area (cement paste) with the mouse; through the selected foreground information, obtain the remaining foreground and background; get the segmentation result after several iterations of calculation. The segmentation results may be different for different ways of foreground selection. As shown in Figure 6, the selection of foreground 2 cannot get a satisfactory segmentation result.Figure 6
The segmentation results of GrabCut algorithm: (a) foreground 1, (b) the segmentation result of foreground 1, (c) foreground 2, and (d) the segmentation result of foreground 2.
(a)(b)(c)(d)


## 2.2.1. Watershed Algorithm Based on Distance Transformation

The foreground and background are labeled with different colours (or grayscale values), and the areas of uncertainty are labeled with 0. Then, the watershed algorithm [14] is implemented to update the labels, and the boundary information of the target is obtained when the foreground and background meet [15]. In this paper, the area of cement paste is the foreground, and the area other than the cement paste is the background. The main steps are as follows:(1)
Binarization [16], preset threshold Q, and pixel points with grayscale values less than Q are set as black, and those greater than or equal to Q are set as white.(2)
Morphological operations [17] remove all white noise in the image.(3)
Distance transformation, generates a distance image that describes the distance between pixel points in the grayscale image and the nearest background pixel.(4)
Label creation, the areas of uncertainty are labeled with 0, and the background is labeled with a positive integer above 1.(5)
Watershed operation, the foreground and background are distinguished by updating the label.(6)
Foreground segmentation, the pixel positions labeled with 1 (background area) are synchronized to the original image.The results of the implementation of each step are shown in Figure4.Figure 4
The segmentation results of the watershed algorithm.


## 2.2.2. HSV Colour Segmentation Algorithm

Compared with the RGB model, the HSV model is easier to track objects of a specific colour [18]. The H-value indicates the hue, the S-value indicates the saturation, and the V-value indicates the luminance. The colour of the cement paste area is relatively uniform and obviously different from other areas. Therefore, a range of pixel values can be obtained by extracting the HSV value of the cement paste area, the pixel value located in this range is set to 255 (white), and the pixel value below or above this range is set to 0 (black). The results are shown in Figure 5.Figure 5
The segmentation results of HSV.


## 2.2.3. GrabCut Foreground Segmentation Algorithm

The GrabCut algorithm is based on colour features and edge features of the image [19]. The premise of image segmentation using the GrabCut algorithm is that there is a significant colour difference between foreground and background. The implementation process is to select the foreground area (cement paste) with the mouse; through the selected foreground information, obtain the remaining foreground and background; get the segmentation result after several iterations of calculation. The segmentation results may be different for different ways of foreground selection. As shown in Figure 6, the selection of foreground 2 cannot get a satisfactory segmentation result.Figure 6
The segmentation results of GrabCut algorithm: (a) foreground 1, (b) the segmentation result of foreground 1, (c) foreground 2, and (d) the segmentation result of foreground 2.
(a)(b)(c)(d)


## 2.3. Calculation of Slump Flow

In this paper, the watershed algorithm based on distance transformation is selected to segment the boundary of the corrected area of cement paste. The extreme points of the boundary are obtained by using the NumPy library [20]. As shown in Figure 7, the maximum diameter of the cement paste area can be calculated according to(4)D=P2−P1+P4−P32,where (P2 − P1) is the distance between points P1 and P2, and (P4 − P3) is the distance between points P3 and P4. After the video is frame-separated, the above operation is performed for all images in turn. The visualization tool Matplotlib is used to draw the relationship curve between slump flow and time, as shown in Figure 8.Figure 7
The boundary and extremum of cement paste.Figure 8
The SF-T curve.


## 3. Automatic Test Method for Flowability of Cement Paste

The mapped desktop colour is not conducive to identifying the boundary of the glass plate, so the bottom of the glass plate can be padded with white paper or sprayed with white paint, and the outer ring wrapped with a yellow rubber band (size in 5 cm–8 cm), as shown in Figure9(a). In addition, the laptop needs to be connected to an external USB camera, and the whole procedure is executed by the laptop, as shown in Figure 9(b).Figure 9
The layout of test scenario: (a) test scenario; (b) sketch map.
(a)(b)The process of automatic measurement of slump flow of cement paste is as follows:(1)
Start-up procedure while the tester lifts the slump cone, and shoot the video of the cement paste flowability test automatically.(2)
Frame the video according to the time series of 100 ms, and store the picture set in the local folder.(3)
The perspective transformation function is used to correct the image based on the automatic point selection method.(4)
For the corrected image, the watershed algorithm is used to identify the boundary of cement paste and determine the extreme points.(5)
The slump flow is calculated, and the SF-T curve is drawn.The flowchart is shown in Figure10.Figure 10
The flowchart.In order to facilitate the conversion between the image pixel distance and the actual size, the size of the corrected glass plate can be set as the actual size, and the distance between the pixel value of the minimum point and the maximum point of the horizontal coordinate is the actual distance.


## 4. The Test Results and Analysis for Slump Flow Test of Cement Paste

Based on the above method, nine groups of flowability tests of cement paste were conducted with the water powder ratio being 0.9, 1.0, and 1.1 and the admixture content was 0.5%, 0.6%, and 0.7%. The results of slump flow by using a steel ruler and the automatic measurement results are shown in Table1. It can be seen that the relative errors of the measurement results are minor and are below 2%. The method can record the slump flow every 100 ms, and the SF-T curve shows that the flow rate decreases gradually with time, which shows a high agreement with the manual measurement and can replace the manual measurement.Table 1
Comparison of SF measurement results of cement paste flowability test.
Water powder ratioAdmixture content (%)Results of steel ruler measurement (mm)Results of automatic measurement (mm)Relative error (%)0.90.5260255.71.650.90.6285280.01.750.90.7300297.70.771.00.5290295.41.861.00.6300298.00.671.00.7310307.50.811.10.5310312.60.841.10.6320320.180.061.10.7330328.70.39To compare with the results of automatic measurement, T200 can be obtained by observing the pictures. For example, if the diameter of the cement paste area in thenth picture after the framing process is exactly 200 mm, then T200 can be calculated with(5)T200=100ms⋅n.Concurrently, this paper uses the algorithm in reference [9] to measure T200. The comparison results of the three methods are shown in Table 2. It can be seen that the relative error between the accurate value and the calculation result of the algorithm proposed in this paper is less than 8%. Compared to the algorithm in reference [9], there is a great improvement in the measurement accuracy. It can be used in practical experiments.Table 2
Comparison of T200 measurement results of cement paste flowability test
Water powder ratioAdmixture content (%)① Results of observing the pictures (mm)② Results of the algorithm in reference [9] (mm)③ Results of the algorithm in this paper (mm)Relative error of ① and ② (%)Relative error of ① and ③ (%)0.90.54600433747005.722.170.90.63900375038003.852.560.90.73300297134009.973.031.00.531002664290014.066.451.00.62500249923000.048.001.00.7230002200100.004.351.10.520001655200017.250.001.10.6130001400100.007.691.10.78000900100.000.13


## 5. Conclusion

In this paper, based on dynamic image processing technology, the automatic detection method for the flowability test of cement paste is studied. This method is verified by nine groups of flowability tests of cement paste. The conclusions are as follows.The corner coordinates of the glass plate before transformation can be obtained fastly and accurately by using the automatic point selection method; the relationship curve between slump flow and time can be output automatically to get the slump flow at any moment; the errors of measurement results about slump flow are below 2%; the errors of measurement results about T200 are below 8%; this method is efficient and straightforward, which can replace the manual measurement method, and provides effective support for the mix proportion optimization system of self-compacting concrete.This work has the following limitations: the method has certain requirements for the brightness of the measurement environment, in addition to the interference of water damage and human shadows, which needs further improvement. In the future, this technology can be further applied to test the workability of concrete and develop a feedback mechanism based on the test results to adjust the mix proportion in real-time and realize smart production of concrete.



---

*Source: 9950524-2022-06-13.xml*