Advances in Building Technologies and Construction MaterialsView this Special Issue
Concrete Image Segmentation Based on Multiscale Mathematic Morphology Operators and Otsu Method
The aim of the current study lies in the development of a reformative technique of image segmentation for Computed Tomography (CT) concrete images with the strength grades of C30 and C40. The results, through the comparison of the traditional threshold algorithms, indicate that three threshold algorithms and five edge detectors fail to meet the demand of segmentation for Computed Tomography concrete images. The paper proposes a new segmentation method, by combining multiscale noise suppression morphology edge detector with Otsu method, which is more appropriate for the segmentation of Computed Tomography concrete images with low contrast. This method cannot only locate the boundaries between objects and background with high accuracy, but also obtain a complete edge and eliminate noise.
The concrete can be considered as a multiphase composite material system, which consists of mortar matrix, aggregate, and interfacial transition zone (ITZ), with various size pores distributed inside the concrete. Pores play a crucial role in strength and durability performance of concrete. For example, porosity is directly related to strength [1–4]. Mean pore size and spacing factor have affected significantly the frost resistance of concrete [5, 6]. Increasing pore connectivity factor will lead to the permeability resistance decreasing of concrete [7, 8]. In recent years, tremendous interests have been aroused by people in pore characteristics, which are crucial for the profound understanding of concrete deterioration mechanism and improving concrete performance.
Due to the diversity pore sizes and shapes distributed in the microstructure of concrete, the analysis of pores becomes difficult. Usually, the analysis of pore characteristics is through the means of indirect methods like mercury intrusion method (MIP) [2–4, 7–9]. However, the hypothesis of MIP method that pore is cylindrical is not always the actual situation and the predrying process of samples before testing could result in the irreversible deformation of pore structure. Some scholars [10–12] argue that MIP method for characterizing pore structure is not appropriate. With the nondestructive analysis technique developing, the Computed Tomography (CT) scanner is being introduced into the characterization of microstructure in materials science fields [13–15].
The computer-aided image processing method may be used to extract a target object in CT images. The image processing method that originated in the 1920s becomes an increasingly powerful tool to solve the hot topic in the civil engineering fields. Başyiğit et al.  evaluated the relationship between microstructure and compressive strength of concrete based on the image analysis software (Image J); Marinoni et al.  employed the threshold segmentation algorithm and the filtering technique to investigate the mortar morphology. Soroushian et al.  carried on the quantification analysis on the microcrack and void in the microstructure of concrete using the automatic thresholding method. One can find that image segmentation methods can mainly be classified into three categories: threshold segmentation algorithm [19–25], edge-based algorithm [26–28], and interdisciplinary application of algorithms [29–32] (e.g., mathematic morphology). However, there is no single accepted method of achieving the segmentation for all images. So finding or matching a new effective segmentation method is indispensable to obtain better segmentation results for CT concrete image.
The paper is organized as follows. First, the preparation of pavement concrete samples and CT testing method are presented in Section 2. Section 3 presents the image methods used in the segmentation and the comparative analysis of these methods. It has been shown that the single traditional method fail to perform segmentation for CT images with low-contrast edge effectively. So a new method, a combination of mathematic morphology and Otsu method, is developed in Section 4. Section 5 draws the conclusions.
2. Experimental Preparation and CT Testing Method
2.1. Sample Preparation
Pavement concrete with the workability slump of 20–40 mm and the strength grade of C30 and C40 is prepared by mixing these raw materials with different proportion. The basic materials are ordinary portland cement 42.5R, river sand, crushed stone (size range 5–26.5 mm), and superplasticizer (0.6–1% recommended dosage and 30% water reduction ratio). Six cube specimens of dimensions 100 × 100 × 100 mm consist of five cube specimens for testing 28-day compressive strength and one for CT scanning. The mixing proportions of raw materials and the compressive strength of specimen are presented in Table 1.
2.2. CT Testing Method
The X-ray image analysis of concrete specimens is carried out nondestructively in the micro-CT scanner, shown in Figures 1(a) and 1(b). As X-rays conduct through the concrete specimens, they are attenuated due to the absorption of objects with different material density. The two-dimensional or three-dimensional image can be reconstructed based on X-ray attenuation information, which is measured by an X-ray muhidetector. Concrete specimens are scanned through micro-CT scanner with a 225 kv voltages, 0.6 mA current, cone beam scan mode, 0.5 micron detection accuracy, and 0.12–0.13 millimeter scanning interval.
3. Traditional Methods for Image Segmentation
The aim of the segmentation for concrete specimens is to separate objects from background image. However, the segmentation of CT images with low contrast is not an easy task. Traditional image segmentation methods are based on the similarities and discontinuities in the gray images [22–28]. In this section, the comparative analysis of threshold algorithms and edge detecting techniques is carried out.
3.1. Threshold Segmentation Methods
The threshold approach based on the gray histogram is a simple practical image segmentation technique, which only requires a gray threshold value, but if the threshold value is not appropriate, the segmentation performance of the method would be directly affected. Common approaches founded on threshold method include two-mode algorithm, iterative algorithm, and Otsu algorithm.
3.1.1. Two-Mode Threshold Algorithm
Prewitt  first proposes the two-mode method, which is a typical global threshold approach. When two peaks appear in the histogram of image, the threshold is generally located at the valley of the gray histogram. So the approach strongly depends on the image operator’s experience and is only applicable to such kind of images that the obvious gray difference exists between background and objects.
3.1.2. Iterative Threshold Algorithm
Iterative algorithm [20, 23] is an adaptive threshold method. The threshold is obtained by optimizing an objective function. The main steps are described as follows: step one, the initial threshold is determined according to formula (1); step two, the image is divided into two regions based on the threshold , and then, the average gray for each region is calculated again; step three, the threshold is updated by a new threshold according to formula (2); repeat step two and step three until and are approximately equal:where is the minimum gray value of image. is the maximum gray value. and are the average gray values of 1st region and 2nd region. is the frequency of occurrences of gray .
3.1.3. Otsu Method
Otsu method also is known as the maximal variance between-class method, which was put forward in 1979. In most cases, Otsu method [21, 22, 24, 25] can obtain better segmentation results. Suppose that the gray is a given image range from 1 to ; the threshold was determined by maximizing the variance between classes according to formula (3). An image may be divided into two classes (objects and background) by the threshold : where is the number of pixels at gray level , is the total number of pixels, consists of pixels with gray levels and consists of pixels with gray levels , and are the cumulative probabilities, and are the variances of the classes and , is the average gray level, and is the variance between classes.
3.1.4. Comparisons of Three Threshold Methods
The above three threshold methods are used to implement segmentation of three C30 (a–c) and three C40 (d–f) CT concrete images. The processed images are shown in Figure 2. Some conclusions can be drawn: using the two-mode algorithm, the threshold value needs to be manually adjusted to obtain better segmentation performance, so it is time-consuming; applying the iterative algorithm, the whole gray level in the image is classified into the background since the final iterative threshold is 0; although the ideal segmentation performance, Otsu method, can be obtained in many applications, Figure 2. shows that the CT concrete images with low-contrast subject to inaccurate segmentation.
3.2. Traditional Edge Detection Algorithms
An edge is the boundary between background and objects, in which the materials density and the image intensity show abrupt changes at edges. The popular edge detection techniques (e.g., Roberts, Prewitt, Sobel, Canny, and LOG operators) have been widely applied to detect discontinuities in gray level. The corresponding algorithms can be expressed by the following formulae, respectively:
Figure 3 shows the results of edge detection using the different edge detectors. It can be detected that the traditional edge operators have good performance in edge detection and localization, but the question reported in the literatures [26–28] remains: these detectors may lead to the loss of boundary information, so the severed edges need to be connected by applying other theories; removing noise ability of the operators like Roberts, Prewitt, and Sobel operators is poor; although the Canny and LOG operators have stronger ability to eliminate noise, part of boundaries are still incomplete; and the threshold of segmentation also needs to be adjusted manually.
(a) Original CT image
(b) Sobel operator
(c) Roberts operator
(d) Prewitt operator
(e) Canny operator
(f) LOG operator
(g) Original CT image
(h) Sobel operator
(i) Roberts operator
(j) Prewitt operator
(k) Canny operator
(l) LOG operator
4. Image Segmentation Based on Mathematical Morphology Method and Otsu Method
Mathematical morphology is an interdisciplinary theory for the analysis and processing of digital images, based on set theory. It was first introduced into image processing fields by Matheron and Serra in 1964. As a powerful image processing tool, now it has been widely used in many domains like medical image processing , remote sensing image analysis , industrial inspection , materials science , and so forth. Morphological image processing consists of a set of operators (e.g., erosion, dilation, opening, and closing) that transform images according to the different geometrical structures characterizations.
4.1. Mathematical Morphology Method
In morphology method, erosion and dilation operators are the basic structuring element sets. Many improved image processing operators are formed by combining the basic operators for obtaining many image processing tasks. For a gray image, let and be the original image (the domain ) and the structuring element (the domain ) respectively, and denotes the translation parameters. The general framework of erosion and dilation is given by
The erosion and dilation of the gray image by the structuring element mean subtracting and adding from and to center , respectively. Let us choose 3 × 3 window as the structuring element ; the erosion and dilation processes of by are described individually in Figures 4 and 5. Suppose that “4” is the center; the erosion of by can be understood as the center of is moved to the point “4,” and then let us subtract every element in window from 4 to find the minimum in the new window . Finally, the point “4” is replaced by the minimum within the new window while the dilation of by means that the point “4” is replaced by returning the maximum value within the new window . The final processed images are obtained by implementing the same process for every point in the gray .
4.2. Morphological Filters
The morphological filter is a nonlinear filter approach based on a series of morphological operators. Compared with the smoothing filter or the sharpening filter, the morphological filter can better filter noise and protect the image details [36–38]. The opening operator and the closing operator are the basic ones in the morphological filter method. A random filter approach is also generated by combinations of the basic operators. The opening operator and the closing operator may be expressed by the formula (10). The opening operator can often filter the peak noise, while the closing operator can filter the valley noise. So one can obtain a new algorithm by combining these basic operators for removing various noise in an image:
To verify the performance of the morphological filter, in this paper we considered the following morphological filter operators including the single opening filtering, the single closing filtering, and the mixed filtering. Experiment results are shown in Figures 6(a)–6(f), in which (a) denotes the original image, (b) is an artificial image added the salt-and-pepper noise, and (c)–(f) are the filtered images. Some conclusions can be drawn: the opening operator can only filter the white point and the closing operator can only filter the black point in the image, while the combinations of opening and closing can remove both the white noise and the black noise.
(a) Original CT image
(b) Image with adding noise
(c) Single opening filter
(d) Single closing filter
(e) Opening followed by closing
(f) Closing followed by opening
4.3. Proposed Segmentation Method by Combining Multiscale Noise Suppression Morphology Edge Detector with Otsu Method (MNSMO)
Many morphology edge detection algorithms [39–41] have been developed for performing a special image processing task. The new operators are generated by the combinations of the basic morphologic operators. For example, you can get the internal boundaries of objects in an image by subtracting the image by erosion from the original image, while the external boundaries of objects can be obtained by subtracting the original image from the image by dilation. Besides, the structuring element dimensions may be adjusted to obtain more subtle results. In the paper, the designed multiscale noise suppression morphology edge detection operator can not only be used for detecting the boundaries between the background and variable size objects in a concrete image, but also can be used for removing noise.
4.3.1. Designing Idea of MNSMO Algorithm
Combinations of the opening operator and the closing operator are used to detect the edge and remove noise in the original . The structuring elements with variable dimensions are used to detect the subtle boundaries information. Through comparing the testing results, the structuring element with the “disk” shape (shown in Figure 7) is more suitable for the pore edge detection. Finally, Otsu algorithm is applied to carry out the binary processing. The designed multiscale noise suppression morphology edge detection operator is shown in the formula below:
4.3.2. Description of MNSMO Algorithm
The flowchart for the MNSMO algorithm is shown in Figure 8.
The main image processing steps are summarized as follows.
Step 1. The original image is filtered by using the morphological operator of the opening followed by the closing.
Step 2. The disk-shaped multiscale structuring elements are designated as 1, 3, 5, and 7.
Step 3. It is the multiscale edge detection through using the erosion and dilation operators.
Step 4. The Otsu method is used to implement the binarization for the above processed image.
4.4. Application of MNSMO Method in the Image Segmentation
The segmentation of six concrete images including three C30 (A1–A3) specimens and three C40 (A1–A3) specimens is implemented by using the MNSMO method in the matlab7.0 platform. Figure 9 illustrates the segmentation results. For each C30 (G1–G3) or C40 (G1–G3) image, the pore edge and the pore shape can be automatically and accurately located. Compared with other image segmentation methods, it can be seen that the MNSMO method has several distinct advantages: noise in the image can be effectively eliminated; both the internal boundaries and the external boundaries of the pores have been detected accurately and no discontinuity of the edge appears in the six images; and by comparison, it is more simple for implementing the segmentation for an image using the MNSMO method that only the structuring element shape and dimensions need be adjusted manually; besides it is not sensitive to the low contrast of the image.
Three threshold methods (e.g., two-mode algorithm, iterative algorithm, and Otsu algorithm) are inaccurate and time-consuming for the segmentation of a CT concrete image. Although traditional edge detection algorithms can better locate the boundaries between objects and background, the edge discontinuity needs to be further solved. And it is poor to restrain noise for some algorithms such as Sobel, Roberts, and Prewitt.
In this paper, a new segmentation method is proposed based on mathematical morphology theory, named “MNSMO.” The MNSMO method is a combination operator of noise suppression operator, edge detection operator, and automatic threshold segmentation algorithm. First, the operator of opening followed by closing is used to filter noise in an image; then the edge is located by designing the eroding and dilating operators by using the multiscale dimensions and the “disk” shape structuring element. Finally, Otsu method is applied to implement the binarization of image.
Compared with the segmentation methods reported in the literatures, “MNSMO” not only can accurately locate the boundaries between objects and background, but is an automatic image processing method, so it avoids subjecting to human error.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This study was financially supported by the Fundamental Research Funds for the Central Universities (nos. 2013G5210010 and 2013G2313001) and the National Natural Science Foundation of China (no. 51278059).
S. Chatterji, “A discussion of the paper ‘Mercury porosimetry—an inappropriate method for the measurement of pore size distributions in cement-based materials’ by S. Diamond,” Cement and Concrete Research, vol. 31, no. 11, pp. 1657–1658, 2001.View at: Google Scholar
M. A. Snyder, D. G. Vlachos, and V. Nikolakis, “Quantitative analysis of membrane morphology, microstructure, and polycrystallinity via laser scanning confocal microscopy: application to NaX zeolite membranes,” Journal of Membrane Science, vol. 290, no. 1-2, pp. 1–18, 2007.View at: Publisher Site | Google Scholar
B. Wei and Z. Liu, “The effective particles edge detection method based on laplace,” Procedia Engineering, vol. 29, pp. 3096–3099, 2012.View at: Google Scholar