Next Article in Journal
Test-Retest Reliability of PODOSmart® Gait Analysis Insoles
Previous Article in Journal
Design of an Automated Multiposition Dynamic Wheelchair
Previous Article in Special Issue
Characterization of Temperature Gradients According to Height in a Baroque Church by Means of Wireless Sensors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Image-Based Automated Width Measurement of Surface Cracking

by
Miguel Carrasco
1,
Gerardo Araya-Letelier
2,
Ramiro Velázquez
3 and
Paolo Visconti
3,4,*
1
Facultad de Ingeniería y Ciencias, Universidad Adolfo Ibáñez, Av. Diagonal las Torres 2640, Santiago 7941169, Chile
2
Escuela de Construcción Civil, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago 7820436, Chile
3
Facultad de Ingeniería, Universidad Panamericana, Av. Josemaría Escrivá de Balaguer 101, Aguascalientes 20296, Mexico
4
Department of Innovation Engineering, University of Salento, Via per Monteroni, 73100 Lecce, Italy
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(22), 7534; https://doi.org/10.3390/s21227534
Submission received: 30 September 2021 / Revised: 6 November 2021 / Accepted: 10 November 2021 / Published: 12 November 2021
(This article belongs to the Special Issue Sensors and Data Processing Techniques for Cultural Heritage)

Abstract

:
The detection of cracks is an important monitoring task in civil engineering infrastructure devoted to ensuring durability, structural safety, and integrity. It has been traditionally performed by visual inspection, and the measurement of crack width has been manually obtained with a crack-width comparator gauge (CWCG). Unfortunately, this technique is time-consuming, suffers from subjective judgement, and is error-prone due to the difficulty of ensuring a correct spatial measurement as the CWCG may not be correctly positioned in accordance with the crack orientation. Although algorithms for automatic crack detection have been developed, most of them have specifically focused on solving the segmentation problem through Deep Learning techniques failing to address the underlying problem: crack width evaluation, which is critical for the assessment of civil structures. This paper proposes a novel automated method for surface cracking width measurement based on digital image processing techniques. Our proposal consists of three stages: anisotropic smoothing, segmentation, and stabilized central points by k-means adjustment and allows the characterization of both crack width and curvature-related orientation. The method is validated by assessing the surface cracking of fiber-reinforced earthen construction materials. The preliminary results show that the proposal is robust, efficient, and highly accurate at estimating crack width in digital images. The method effectively discards false cracks and detects real ones as small as 0.15 mm width regardless of the lighting conditions.

1. Introduction

Over the last two decades, several techniques have been proposed to detect crack formation in materials and structures together with estimating its dimensions (width, length, and depth). The first image-based algorithms devoted to crack detection explored segmentation techniques targeting different types of civil engineering materials and infrastructures such as concrete, walls, and load testbeds. Recently, numerous advances have been made to develop techniques that consider the time evolution of loading with Digital Image Correlation (DIC) techniques [1,2] or hybrid methods that use a set of techniques and their combinations to improve the segmentation process [3,4,5]. In this way, it has been possible to perform experiments on different types of samples, which in some cases have allowed the evaluation of the use of different fiber types (e.g., synthetic, vegetal or animal) on different construction materials (e.g., mortar, concrete, and adobe) [6,7,8], as well as the study of the loading conditions over time on columns [9], beams and slabs [3], and walls [4] and the design of analytical models for estimating crack width [10] and surface analysis using machine learning techniques [11,12,13,14].
One of the main challenges for extracting and isolating the crack from the rest of the image is the implementation of segmentation algorithms. Research in this field has focused on two paradigms: (1) algorithms based on heuristic image processing and (2) algorithms based on Deep Learning (DL) tools.
The former includes techniques that evaluate the changes in the structure as a function of load [1,15,16,17], techniques that use one or more cameras for homography [18] or photogrammetry [15,19,20,21,22], and those that simulate a loading process and then analyze the material under pressure [19,23]. In general, the proposed solutions can be classified into techniques based on region growth [16,22], combinations of mathematical morphology and optimal threshold selection [5,24,25], techniques based on the Hessian matrix [26,27], and methods based on a sequence of partial filters that perform noise reduction as segmentation goes on [4,5,13]. The main advantage of techniques based on heuristic processing is that they can be directly applied to the samples without requiring any previous training. Nevertheless, they strongly depend on the type of material, implying that they might not always work properly.
The latter techniques have recently raised interest among researchers mainly because DL allows the implementation of robust solutions that can be applied in a wide variety of materials/structural elements such as paving, walls, and columns. In addition, it allows the use of different illumination conditions, orientation, distance, and objects [27,28,29,30]. DL-based algorithms rely on a learning process, which involves a previous labeling of images, so that the network can make a prediction through the learning process. Several types of network architectures have been explored in the literature [14,30,31]. The main difference across them is the generated output: images composed of blocks classified as failures [28,30], pixel-level binary segmentation [29,32,33,34], or a combination of both [34,35].
Both heuristic and DL paradigms exhibit different types of responses and require specific adjustments according to the exploited approach [31,36]. Surprisingly, despite the great advances in terms of segmentation, the literature presents limited progress with respect to crack width analysis [26]. Most related work found in the literature addresses crack width and length estimation without considering that width might change along with the material [27,29,33]. In these works, the middle section of the crack is indeed analyzed. However, no technique to measure the width as it progresses along the crack’s path is observed. An accurate image-based estimation of the crack width remains an open problem [9,24,31]. It is important to mention that different crack widths can be related to very different mechanical/durability effects within civil engineering materials [37]. Thus, an accurate crack width estimation is crucial since very fine cracks might be harmless and do not require further attention while wider cracks might require immediate attention and action. One of the main challenges in this topic is how to correctly measure the width while taking into consideration that the crack might exhibit a highly variable path, and therefore its width must be measured orthogonally to the path. Once this problem is solved, it becomes then relevant to consider the number of points to evaluate along the path.
Currently, the most popular method for crack measurement and characterization is by visual inspection using a crack-width comparator gauge (CWCG) [38], where the width of the crack is recorded manually over a distance interval using this CWCG (Figure 1). Although this technique allows for the rapid measurement of cracks in materials and structures, it is prone to serious errors: (i) the overlap between the measuring rule and the crack is not accurate as it depends on the inspector’s experience; (ii) depending on the orientation of the crack curvature, the crack measurement angle may not be accurate as it might not be tangentially measured; (iii) the measurement interval between two points might not be regular and may under or over represent the average crack width as cracks do not necessarily have a normal distribution [9] and (iv) the measurement is made by visual observation; therefore, it is subject to the inspector’s subjective measurement [39]. The above problems lead to the fact that the crack measurement can be underestimated or overestimated.
To avoid these errors, it is necessary to perform a non-contact measurement of the sample in such a way that there is no intervention on the crack to be measured. This requires the acquisition of images in which cracks are clearly exposed. However, there are some limitations inherent to the acquisition; for example, the incident angle of the light could generate shadows on the surface (Figure 2) that could lead to errors in crack detection. Depending on the material under analysis, the surface may generate false crack regions [31]. In this context, it is necessary to establish the correct position of the light source for the object to be analyzed. In general, the light should be positioned perpendicularly to the object to eliminate shadows and avoid false crack regions. This situation has been analyzed with different approaches, all of which are strongly dependent on the type of sample to be used. However, techniques using DL-based segmentation have proven to be more versatile for handling samples of different types of material, angles, and distances [30].
This paper proposes a novel image-based method for estimating crack width with a three-stage algorithm that effectively considers width variations along the crack’s path. By measuring the orientation angle that is normal to the crack curvature, the method allows for a more accurate estimation of the cracks present in the material or structure. Much of the work found in the literature has not focused on this problem. On the contrary, research has focused on solving the segmentation problem. Although segmentation is a highly relevant and complex problem, so is width estimation. Its correct measurement is essential to determine the actual condition of civil structures. Automating this process can be done after addressing the segmentation problem.
The rest of the paper is organized as follows: Section 2 describes the method for crack width estimation. Section 3 presents the results obtained with different sample materials. Finally, Section 4 concludes by summarizing the main concepts and results and giving perspectives on future work.

2. Method

To determine the width of the cracks on an image, the proposed algorithm exploits three stages: (i) preliminary filtering; (ii) adaptive segmentation; and (iii) profile estimation (Figure 3). This section describes these stages. The following subsections describe each of these stages together with their corresponding intermediate steps.

2.1. Step I: Preliminary Filtering

One of the main challenges in the crack detection process is to reduce the false cracks, i.e., the regions that have visual characteristics similar to those of a real crack, such as shadows or elements embedded in the material, such as fibers, and digital noise. Therefore, the objective of the preliminary filtering stage is to reduce the number of false cracks through two filtering sub-stages: (i) color space change; and (ii) Perona–Malik filter.
The first sub-stage involves a color change using the well-known L*a*b color space [40]. This technique maximizes the luminance channel while allowing the attenuation of the effects due to shadows or light levels on the object. For this purpose, the algorithm performs an L*a*b channel transformation on the original RGB image. Only the L-channel of this transformation is retained.
The second sub-stage uses the filtering technique proposed by Perona and Malik [41,42,43]. This filter achieves noise reduction while preserving the fundamental structures of the image such as the edges. Let us consider a grayscale image L of dimensions (n × m) represented by gray values f x , y 2 in positions (x, y). The proposed filter is defined by Equation (1):
t L = d i v c L   L ,   w i t h   t = 0 L = L
where div(∙) represents the divergence operator, ∇ is the gradient operator, and c(∙) is a non-increasing smooth diffusion operator with c(0) = 1 and c(s)→0 with s→∞, which depends on the location of the point in the image.
According to variable t, this operation is iterative. To achieve a result, it requires a number of previously user-defined iterations. In our proposal, the number of iterations is a fixed parameter. Concerning the diffusion operator, different models can be found in the literature [41]. In this work, we implement the linear model proposed by Weickert [44] that uses a technique known as coherence filter [45]. Figure 4 shows a performance example of the filtering stage.

2.2. Step II: Binary Segmentation

Let us assume that the result of the previous stage stopped at time t. In addition, let us consider W = t L as the image resulting from the coherence filter. Applying operator · over matrix W generates a vector w = W of dimensions (nm × 1) (see Appendix A). The next procedure involves the search of threshold δ, which relates each value of vector w to a binary value, in other words, the process known as segmentation. To that end, let us further define vector w as a monotonically increasing function with values ordered through an order function w + = s o r t w . The threshold’s value can be determined by Equation (2):
δ = a r g m a x i w + i n m i i = 1 , , n m
Image segmentation can be defined as the process of relating all pixels in the image to a binary value f x , y > δ Β x , y = 1 , otherwise Β x , y = 0 , where Β x , y is a matrix of dimensions (n × m) (Figure 5).
Once the binary edges have been determined, our proposal proceeds to determine the central position of the potential cracks. For this purpose, the method employs the skeleton algorithm based on the distance transformation process [46], which allows the simplification of the image to determine the centroids of each crack.
The distance transform is a measure of the distance between a pixel and its edge. The farther the pixel from the edge, the higher its value [24,47]. Subsequently, the method applies the top-hat filter in such a way that only the pixels with the maximum value remain stable. This procedure eases the application of the topological filter to find the skeleton of the structure (Figure 6). Even though it is possible to directly apply this last filter to the image shown in Figure 6a, the resulting skeleton would be incorrect since other regions that are not part of the central structure of the crack could appear given that the crack exhibits a greater width in other regions.
Let S be the binary matrix containing the crack skeleton (Figure 6d) and M a coordinate matrix satisfying Equation (3):
M i = x 1 , . . , m y 1 , . . , n     x , y   |   S x , y = 1
Depending on the number of coordinates having the logical value 1, the number of indices may vary for each p-tuples. Since the pixels of the skeleton are neighbors to each other, one way to reduce the spacing is to determine a minimum distance dmin corresponding to the minimum distance between a point and its neighbor.
The proposed procedure iterates the number of clusters of the k-means algorithm [48] until a certain distance between the centers is reached. Assume that C = k m e a n s M ,   k is a subset of tuples of matrix M, i.e., C M for k centroids. The purpose is to find distance dmin according to k given in Equation (4):
d m i n > a r g m i n k   C i C j   i , j   |   i j k +
Figure 7 illustrates this concept.

2.3. Step III. Profile Width

Once the central point and neighboring spacing have been determined in the image, the next step involves the calculation of an optimal orientation angle that is normal to the crack curvature at each point. For this purpose, we use a 40 × 40-pixel mask centered on each point described above, and over this point, we extract a binary sub-region, previously defined as matrix B. The central angle can be estimated as the angle of the major axis of an ellipse that fits the boundary of the binary region (details of this algorithm can be found in [49]). This ensures that the measurement of the ellipse’s width is correct and in accordance with the path and curvature of the crack (Figure 8a). Once both the angle and the position have been defined, the method proceeds to determine the width by superimposing a profile line greater than the width of the crack (Figure 8b).
As the cracks have a significantly different intensity level with respect to the background, it is possible to determine a cluster that identifies two groups of pixels {background, crack} through the k-means algorithm. Let Pi be the i-th vector comprising the intensity levels of length t. The width of the crack can be calculated as (Equation (5)):
w i d t h i = ( k m e a n s P i 1 ,   t / 2 = = c l a s s _ c r a c k + ( k m e a n s P i t 2 + 1 , t = = c l a s s _ c r a c k
where class crack is previously defined as the lower intensity value of the two groups, and the values of the vector indices range from 1 to t.

3. Results

The proposed method was evaluated on a set of earthen construction material (ECM) samples since this type of construction material is very susceptible to cracking [50,51] and, thus, different types of fibers (i.e., vegetal, synthetic, and animal) are used to reduce cracking and to compare this cracking reduction to an ECM sample without fibers (plain mixture). As previously explained, the algorithm searches for the optimal angle of orientation on the crack’s curvature according to a user-defined spacing for estimating its width and measurement angle, in the same way as a human operator would do with a crack monitor (Figure 9). Since each crack is locally analyzed, the computing time varies according to the number of cracks present in the sample.
It is important to highlight that the number of width samples can be modified by the user implying that the resulting crack width will vary according to the number of samples (Figure 10). Due to sample variation, it is possible to estimate the error and the standard deviation. Therefore, it is possible to obtain a confidence interval for the measurement of the samples and to estimate the optimal number for such measurements. This fact is relevant since it allows obtaining a more accurate measurement compared to the current methods, which are limited in the number of crack width samples.
The dataset used to validate the proposed method was generated from an ongoing study that analyzed the influence of the width of different fiber types (i.e., vegetal, synthetic, animal) in terms of the cracking reduction of ECMs that were also compared to a plain ECM (no fiber) (Figure 11). Depending on factors such as fiber type and dosages, it is possible to observe different cracking patterns, which generate a variable number of cracks, angles, and widths [51,52]. Additionally, to evaluate the method’s performance under different scenarios, samples were exposed to both natural and artificial lighting.
Despite the different lighting conditions in the image dataset, the proposed method exhibits a stable behavior due to the extraction of the luminance channel from the L*a*b transformation, thus reducing the effect produced by the type of light source in the sample (see Figure 12).
On the other hand, in noisy areas generated by non-regular surfaces, the coherence filter allows a reduction in most of the structures that do not correspond to real cracks, thus maintaining their continuity and reducing the number of false alarms (Figure 13c). This effect is more noticeable in some samples of vegetal fibers (V/2 and V/3 fiber), for which surfaces exhibit irregular zones that could be erroneously segmented as cracks (Figure 13c). However, a possible undesired effect is that some measurement points can be filtered out when the crack width is very small (less than 0.15 mm).
In general, the results indicate that the crack samples, which exhibited different types of fibers and lighting, were correctly analyzed by the proposed method and that the width of the cracks was determined with high accuracy (Figure 13b). Furthermore, for those areas in which there is no relevant change between the surface and the background, the method discards the possible candidate regions (Figure 13c). This is achieved upon the analysis of the distance between the crack’s background and surface. The above concept is represented by purple (fuchsia) lines without the blue markers.
To statistically validate the proposed method, a manual width estimation process was performed on the images of the dataset. A total of 330 crack width measurements were extracted. For each sample, 30 crack width measurements were obtained taking into account the crack curvature. Each of the manual points was randomly selected with respect to the dispersion of the cracks in the samples. Both results, manual and automatic, were compared by means of a t-test with non-equivalent variance (Welch Test [53]).
Table 1 summarizes the results obtained. Note that each sample comes from the same distribution according to the values obtained in Z-score and p-value for each of the samples (manual versus automatic). The most significant differences were found in those images with a low point extraction. In such images, a higher number of samples is necessary to make a valid comparison (Figure 11, sample V/4).
According to the conditions defined in Table 1, the crack width distribution was evaluated to analyze its behavior (Figure 14). Recall that the method performs a point spacing between the set of measurements, thus ensuring that each measurement is completely unbiased in a parameterized way (parameter k). The differences in crack width across the samples are associated with a mechanical stress present in the interaction between the fiber and the mixture [50,52]. Still, it is possible to visualize zero, one, or more underlying distributions depending on the type of crack present in the sample.
We discuss the possible types of distributions below.
(1)
No clear distribution: Samples V/3 and V/4 show this case. In these samples, it was possible to extract only a limited number of measurement points due to the small number of cracks present in the sample. Thus, it is not possible to associate the data with a specific distribution. This limitation of the method is a particular point for future improvement.
(2)
Normal distribution: This phenomenon is present in samples with a regular cracking pattern. It is clearly observed in samples A/1, A/2, and A/3. Additionally, in sample A/2, the number of measurements was increased from 75 to 606 through parameter k (Figure 10). By increasing the number of measurements, it is possible to observe a normal behavior and a lower error in relation to the manual process (p-value 0.906 versus p-value 0.308).
(3)
Bimodal distribution: This phenomenon can be clearly observed in samples No Fiber/1, A/4, and A/5 (Figure 13). These samples exhibit two types of cracks (coarse and fine), which have a vertical and/or horizontal cracking behavior. In some cases, a greater width was found in the horizontal cracks. However, it is worth mentioning that all the cracks have internal angles that can only be appreciated in images with a higher magnification (Figure 13a).
The previous analysis is complemented with a boxplot (Figure 15). Overall, it can be observed that the cracks exhibit a variation of less than five pixels (0.75 mm) in the IQR range for each sample. The above does not apply to samples A/4 and A/5, in which the cracking pattern has a larger variation; wide open and closed crack widths can be found in the same sample with an approximate variation of 15 pixels (2.25 mm). It is important to note that the proposed method was evaluated on cracks with a width greater than one pixel (1 pixel = 0.15 mm); however, this relationship depends on the image resolution, which can be improved if the resolution of the sample images is increased. In addition, it can be observed that the median and mean in most of the analyzed samples are very close to each other, thus ensuring the numerical stability of the method by greatly reducing the number of false alarms (see Figure 14, median: green dotted line, mean: orange line). As previously discussed, the main variations can be found in samples V/3 and V/4 due to an insufficient number of measurements to define a distribution. In the rest of the distributions, the behavior is normal, either with univariate or with bivariate distributions due to the type of cracks present in the samples.

4. Conclusions

This paper has presented an efficient and highly accurate method for estimating crack width in digital images and this method can be easily applied to the cracking assessment of construction materials, such as mortar, concrete and adobe, which is an important task to evaluate in civil engineering materials and structures. The proposed method considers the orientation of the crack curvature and ensures that the measurement is unbiased and parameterized (i.e., at regular intervals) and was evaluated on a dataset of earthen construction material (ECM) samples that present samples reinforced with different fiber types (i.e., vegetal, synthetic, animal) as well as a plain sample (no fiber).
The dataset presents different crack widths that are within the range of 0.15 to 5 mm. Images were digitized using two types of light sources (natural and artificial) to evaluate the method’s performance in real situations (laboratory and external environment). For external environments, two additional lighting conditions were used: sunlight and afternoon sun to check the effect of shadows on the cracks.
The results show that there is no statistical difference between the manual and automated measurement. This represents a great advantage since the measurement time of the proposed method is at least 20 times shorter than that of the manual process. In addition, it was possible to obtain a robust confidence interval due to the larger number of measurement points compared to the manual method, which is tedious, time-consuming, prone to human error, and potentially biased. The method is independent to the type of light used, which facilitates the analysis of samples in outdoor environments.
The main difference between our method and DL-based techniques is that the latter focuses exclusively on the segmentation stage. It is true that segmentation is fundamental for crack detection, and our method acknowledges the importance of involving a segmentation stage. However, determining how best to estimate crack width remains an open problem with just a few works in the literature addressing this topic. Most of the methods already reported in the literature perform crack analysis by considering groups of pixels without any continuity. However, cracks exhibit a tortuous behavior, and it becomes essential to preserve their path along the material. Threshold-based algorithms directly applied on the image tend to break this continuity thus limiting any further analysis on the crack width. This limitation justifies the use of the Perona–Malik technique in our proposal which allows us to highlight the path’s continuity while reducing noise (Table 2, noise removal column).
Another major advantage of our method lies in the estimation of the actual crack width. Some works assume that width can be measured either horizontally or vertically and do not consider the crack’s motion or angle. A limited number of studies report on this analysis but do so without considering the crack as an object of variable width. This last concept is fundamental as, at the pixel level, the crack can be visualized as a group of pixels with tortuous motion, with width and angle variations according to the type of stress applied to the sample (Table 2, route tracing column).
Finally, width estimations must be performed on the whole sample at regular intervals, so that the measurement result is unbiased. Most of the previous studies do not take this point into consideration; they only analyze crack bodies, or they only perform an analysis at the pixel level and do not measure the real width of the cracks as suggested by the international building codes (Table 2, spatial width distribution).
The major limitation of the proposed method can be appreciated in samples with limited crack width, which considerably reduces the number of measurement samples. In these cases, in future work, an analysis of the type of crack should be defined prior to the definition of the number of measurements on each sample.

Author Contributions

Conceptualization, M.C. and G.A.-L.; methodology, M.C. and R.V.; software, M.C. and R.V.; validation, G.A.-L., M.C. and R.V.; formal analysis, P.V.; investigation, M.C., G.A.-L. and R.V.; resources, P.V.; data curation, R.V.; writing—original draft preparation, R.V and G.A.-L.; writing—review and editing, R.V. and G.A.-L.; visualization, M.C.; supervision, P.V.; project administration, P.V.; funding acquisition, P.V. All authors have read and agreed to the published version of the manuscript.

Funding

G.A. wants to thank the financial support provided by Pontificia Universidad Católica de Chile grants PIA-3581-050-81 and Proyecto de investigación interdisciplina II190041.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

https://github.com/mlacarrasco/crack_width (accessed on 30 September 2021).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Let X be a matrix of dimensions (n × m) and let ei be the i-th canonical basis vector of dimensions (m × 1), i.e., ei = [0, ⋯, 0, 1, 0, ⋯, 0] with a value of 1 at the i-th position and the rest with zeros. Operator · can be defined as Equation (A1):
X = i = 1 m e i I X e i
where ⊗ represents the Kronecker product, I is the identity matrix of dimensions (n × n). The result of this operation for every matrix X of dimensions (n × m) is transformed into a vector of length (nm × 1). The inverse process, i.e., transforming a vector of length (nm × 1) into a matrix of dimensions (nm × 1) can be defined by Equation (A2):
( X ) = i = 1 m e i I ( X ) e i T = X
We will write · as the vector-matrix transformation process, it being consistent with the number of existing elements in the vector and resulting in a matrix of dimensions (n × m).

References

  1. Zhao, P.; Zsaki, A.M.; Nokken, M.R. Using digital image correlation to evaluate plastic shrinkage cracking in cement-based materials. Constr. Build. Mater. 2018, 182, 108–117. [Google Scholar] [CrossRef]
  2. Bertelsen, I.M.G.; Kragh, C.; Cardinaud, G.; Ottosen, L.M.; Fischer, G. Quantification of plastic shrinkage cracking in mortars using digital image correlation. Cem. Concr. Res. 2019, 123, 105761. [Google Scholar] [CrossRef]
  3. Davoudi, R.; Miller, G.R.; Kutz, J.N. structural load estimation using machine vision and surface crack patterns for shear-critical RC beams and slabs. J. Comput. Civ. Eng. 2018, 32, 04018024. [Google Scholar] [CrossRef]
  4. Hoang, N.-D. Image processing-based recognition of wall defects using machine learning approaches and steerable filters. Comput. Intel. Neurosci. 2018, 2018, 7913952. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Han, H.; Deng, H.; Dong, Q.; Gu, X.; Zhang, T.; Wang, Y. An advanced Otsu method integrated with edge detection and decision tree for crack detection in highway transportation infrastructure. Adv. Mater. Sci. Eng. 2021, 2021, 9205509. [Google Scholar] [CrossRef]
  6. Araya-Letelier, G.; Maturana, P.; Carrasco, M.; Antico, F.C.; Gómez, M.S. Mechanical-damage behavior of mortars reinforced with recycled polypropylene fibers. Sustainability 2019, 11, 2200. [Google Scholar] [CrossRef] [Green Version]
  7. Concha-Riedel, J.; Araya-Letelier, G.; Antico, F.C.; Reidel, U.; Glade, A. Influence of jute fibers to improve flexural toughness, impact resistance and drying shrinkage cracking in adobe mixes. In Earthen Dwellings and Structures; Reddy, B.V.V., Mani, M., Walker, P., Eds.; Springer: Singapore, 2019; ISBN 978-981-13-5883-8. [Google Scholar]
  8. Araya-Letelier, G.; Gonzalez-Calderon, H.; Kunze, S.; Burbano-Garcia, C.; Reidel, U.; Sandoval, C.; Bas, F. Waste-based natural fiber reinforcement of adobe mixtures: Physical, mechanical, damage and durability performance assessment. J. Clean. Prod. 2020, 273, 122806. [Google Scholar] [CrossRef]
  9. Bertelsen, I.M.G.; Ottosen, L.M.; Fischer, G. Quantitative analysis of the influence of synthetic fibres on plastic shrinkage cracking using digital image correlation. Constr. Build. Mater. 2019, 199, 124–137. [Google Scholar] [CrossRef]
  10. Löfgren, I. Calculation of crack width and crack spacing. In Proceedings of the Nordic Mini-Seminar: “Fibre Reinforced Concrete”, Trondheim, Norway, 15 November 2007. [Google Scholar]
  11. Tang, P.; Huber, D.; Akinci, B. Characterization of laser scanners and algorithms for detecting flatness defects on concrete surfaces. J. Comput. Civ. Eng. 2011, 25, 31–42. [Google Scholar] [CrossRef]
  12. Pragalath, H.; Seshathiri, S.; Rathod, H.; Esakki, B.; Gupta, R. Deterioration assessment of infrastructure using fuzzy logic and image processing algorithm. J. Perform. Constr. Facil. 2018, 32, 04018009. [Google Scholar] [CrossRef]
  13. Nishikawa, T.; Yoshida, J.; Sugiyama, T.; Fujino, Y. Concrete crack detection by multiple sequential image filtering: Concrete crack detection by image processing. Comput. Aided Civil Infrastruct. Eng. 2012, 27, 29–47. [Google Scholar] [CrossRef]
  14. Hsieh, Y.-A.; Tsai, Y.J. Machine learning for crack detection: Review and model performance comparison. J. Comput. Civ. Eng. 2020, 34, 04020038. [Google Scholar] [CrossRef]
  15. Christen, R.; Bergamini, A.; Motavalli, M. High precision measurement of surface cracks using an optical system. Meas. Sci. Technol. 2009, 20, 20–077001. [Google Scholar] [CrossRef]
  16. Yamaguchi, T.; Nakamura, S.; Saegusa, R.; Hashimoto, S. Image-based crack detection for real concrete surfaces. IEEJ Trans. Electron. Eng. 2008, 3, 128–135. [Google Scholar] [CrossRef]
  17. Dare, P.; Hanley, H.; Fraser, C.; Riedel, B.; Niemeier, W. An operational application of automatic feature extraction: The measurement of cracks in concrete structures. Photogramm. Rec. 2002, 17, 453–464. [Google Scholar] [CrossRef]
  18. Hartley, R.; Zisserman, A. Multiple View Geometry in Computer Vision; Cambridge University Press: New York, NY, USA, 2003. [Google Scholar]
  19. Detchev, I.; Habib, A.; El-Badry, M.; Bahreh, V.M. Detection of cracks in a concrete beam-column joint using target gridding. In Proceedings of the American Society of Photogrammetry and Remote Sensing Annual Conference, Baltimore, MD, USA, 24–28 March 2013; p. 9. [Google Scholar]
  20. Sohn, H.-G.; Lim, Y.-M.; Yun, K.-H.; Kim, G.-H. Monitoring crack changes in concrete structures. Comput. Aided Civil Eng. 2005, 20, 52–61. [Google Scholar] [CrossRef]
  21. Habib, A.; Lichti, D.; El-Badry, M.; Detchev, I. Optical remote sensing systems for structural deflection measurement and crack characterization. In Proceedings of the 9th International Conference on Short and Medium Span Bridges, Calgary, AB, Canada, 15–18 July 2014. [Google Scholar] [CrossRef]
  22. Yamaguchi, T.; Hashimoto, S. Fast crack detection method for large-size concrete surface images using percolation-based image processing. Mach. Vis. Appl. 2010, 21, 797–809. [Google Scholar] [CrossRef]
  23. Chen, L.-C.; Shao, Y.-C.; Jan, H.-H.; Huang, C.-W.; Tien, Y.-M. measuring system for cracks in concrete using multitemporal images. J. Surv. Eng. 2006, 132, 77–82. [Google Scholar] [CrossRef]
  24. Berrocal, C.G.; Löfgren, I.; Lundgren, K.; Görander, N.; Halldén, C. Characterisation of bending cracks in R/FRC using image analysis. Cem. Concr. Res. 2016, 90, 104–116. [Google Scholar] [CrossRef] [Green Version]
  25. Hoang, N.-D. Detection of surface crack in building structures using image processing technique with an improved otsu method for image thresholding. Adv. Civil Eng. 2018, 2018, 3924120. [Google Scholar] [CrossRef] [Green Version]
  26. Cho, H.; Yoon, H.-J.; Jung, J.-Y. Image-based crack detection using crack width transform (CWT) algorithm. IEEE Access 2018, 6, 60100–60114. [Google Scholar] [CrossRef]
  27. Chen, Q.; Huang, Y.; Sun, H.; Huang, W. Pavement crack detection using hessian structure propagation. Adv. Eng. Inf. 2021, 49, 101303. [Google Scholar] [CrossRef]
  28. Alipour, M.; Harris, D.K. Increasing the robustness of material-specific deep learning models for crack detection across different materials. Eng. Struct. 2020, 206, 110157. [Google Scholar] [CrossRef]
  29. Kang, D.; Benipal, S.S.; Gopal, D.L.; Cha, Y.-J. Hybrid pixel-level concrete crack segmentation and quantification across complex backgrounds using deep learning. Autom. Constr. 2020, 118, 103291. [Google Scholar] [CrossRef]
  30. Rao, A.S.; Nguyen, T.; Palaniswami, M.; Ngo, T. Vision-Based Automated Crack Detection Using Convolutional Neural Networks for Condition Assessment of Infrastructure. Struct. Health Monit. 2021, 20, 2124–2142. [Google Scholar] [CrossRef]
  31. Munawar, H.S.; Hammad, A.W.A.; Haddad, A.; Soares, C.A.P.; Waller, S.T. Image-based crack detection methods: A review. Infrastructures 2021, 6, 115. [Google Scholar] [CrossRef]
  32. Yamane, T.; Chun, P. Crack detection from a concrete surface image based on semantic segmentation using deep learning. ACT 2020, 18, 493–504. [Google Scholar] [CrossRef]
  33. Park, S.E.; Eem, S.-H.; Jeon, H. Concrete crack detection and quantification using deep learning and structured light. Constr. Build. Mater. 2020, 252, 119096. [Google Scholar] [CrossRef]
  34. Liu, J.; Yang, X.; Lau, S.; Wang, X.; Luo, S.; Lee, V.C.; Ding, L. Automated pavement crack detection and segmentation based on two-step convolutional neural network. Comput. Aided Civil Infrastr. Eng. 2020, 35, 1291–1305. [Google Scholar] [CrossRef]
  35. Abdellatif, M.; Peel, H.; Cohn, A.G.; Fuentes, R. Combining block-based and pixel-based approaches to improve crack detection and localisation. Autom. Constr. 2021, 122, 103492. [Google Scholar] [CrossRef]
  36. Mohan, A.; Poobal, S. Crack detection using image processing: A critical review and analysis. Alex. Eng. J. 2018, 57, 787–798. [Google Scholar] [CrossRef]
  37. Alexander, M.G.; Bentur, A.; Mindess, S. Durability of Concrete: Design and Construction; CRC Press: Boca Raton, FL, USA, 2017; ISBN 978-1-4822-3725-2. [Google Scholar]
  38. ACI Committee 224 American Concrete Institute. Control of Cracking in Concrete Structures ACI 224 R-01; American Concrete Institute: Farmington Hills, OK, USA, 2001; ISBN 978-0-87031-056-0. [Google Scholar]
  39. Newman, T.S.; Jain, A.K. A survey of automated visual inspection. Comput. Vis. Image Underst. 1995, 61, 231–262. [Google Scholar] [CrossRef]
  40. Chen, H.-C.; Chien, W.-J.; Wang, S.-J. Contrast-based color image segmentation. IEEE Signal Process. Lett. 2004, 11, 641–644. [Google Scholar] [CrossRef]
  41. Perona, P.; Malik, J. Scale-space and edge detection using anisotropic diffusion. IEEE Trans. Pattern Anal. Mach. Intel. 1990, 12, 629–639. [Google Scholar] [CrossRef] [Green Version]
  42. Kamalaveni, V.; Rajalakshmi, R.A.; Narayanankutty, K.A. Image denoising using variations of perona-malik model with different edge stopping functions. Procedia Comput. Sci. 2015, 58, 673–682. [Google Scholar] [CrossRef] [Green Version]
  43. Liu, K.; Tan, J.; Ai, L. Hybrid regularizers-based adaptive anisotropic diffusion for image denoising. SpringerPlus 2016, 5, 404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Weickert, J. Anisotropic Diffusion in Image Processing; B.G. Teubner: Stuttgart, Germany, 1998. [Google Scholar]
  45. Barbu, T. Robust anisotropic diffusion scheme for image noise removal. Procedia Comput. Sci. 2014, 35, 522–530. [Google Scholar] [CrossRef] [Green Version]
  46. Telea, A.; van Wijk, J.J. An Augmented Fast Marching Method for Computing Skeletons and Centerlines. In Proceedings of the IEEE TCVG Symposium on Visualization, Barcelona, Spain, 27–29 May 2002. [Google Scholar]
  47. Gonzalez, R.; Woods, R.E. Digital Image Processing; Addision Wesley: San Francisco, CA, USA, 1992. [Google Scholar]
  48. Dhanachandra, N.; Manglem, K.; Chanu, Y.J. Image Segmentation using K-means clustering algorithm and subtractive clustering algorithm. Procedia Comput. Sci. 2015, 54, 764–771. [Google Scholar] [CrossRef] [Green Version]
  49. Fitzgibbon, A.W.; Pilu, M.; Fisher, R.B. Direct least squares fitting of ellipses. In Proceedings of the 13th International Conference on Pattern Recognition, Vienna, Austria, 25–29 August 1996; Volume 1, pp. 253–257. [Google Scholar]
  50. Araya-Letelier, G.; Antico, F.C.; Concha-Riedel, J.; Glade, A.; Wiener, M.J. Effectiveness of polypropylene fibers on impact and shrinkage cracking behavior of adobe mixes. In Earthen Dwellings and Structures Effectiveness; Reddy, B.V.V., Mani, M., Walker, P., Eds.; Springer: Singapore, 2019; ISBN 978-981-13-5883-8. [Google Scholar]
  51. Araya-Letelier, G.; Concha-Riedel, J.; Antico, F.C.; Valdés, C.; Cáceres, G. Influence of natural fiber dosage and length on adobe mixes damage-mechanical behavior. Constr. Build. Mater. 2018, 174, 645–655. [Google Scholar] [CrossRef]
  52. Araya-Letelier, G.; Antico, F.C.; Burbano-Garcia, C.; Concha-Riedel, J.; Norambuena-Contreras, J.; Concha, J.; Saavedra Flores, E.I. Experimental evaluation of adobe mixtures reinforced with jute fibers. Constr. Build. Mater. 2021, 276, 122127. [Google Scholar] [CrossRef]
  53. Ruxton, G.D. The unequal variance T-test is an underused alternative to student’s t-test and the Mann–Whitney U test. Behav. Ecol. 2006, 17, 688–690. [Google Scholar] [CrossRef]
Figure 1. Crack-width comparator gauge (CWCG) for measuring the crack width; image resolution 14.7 pixels/mm.
Figure 1. Crack-width comparator gauge (CWCG) for measuring the crack width; image resolution 14.7 pixels/mm.
Sensors 21 07534 g001
Figure 2. Test image with cracks and shadowed regions; original image area (500 mm × 500 mm), sub-image area (82 mm × 82 mm), resolution 11.5 pixel/mm, shadow region area (25 mm × 25 mm): Arrow: detail of a suspected fault region that is actually a shadow.
Figure 2. Test image with cracks and shadowed regions; original image area (500 mm × 500 mm), sub-image area (82 mm × 82 mm), resolution 11.5 pixel/mm, shadow region area (25 mm × 25 mm): Arrow: detail of a suspected fault region that is actually a shadow.
Sensors 21 07534 g002
Figure 3. The proposed three-stage algorithm for crack width estimation. The partial results for each imagen of the algorithm are described in Section 2.1,Section 2.2,Section 2.3.
Figure 3. The proposed three-stage algorithm for crack width estimation. The partial results for each imagen of the algorithm are described in Section 2.1,Section 2.2,Section 2.3.
Sensors 21 07534 g003
Figure 4. Preliminary filtering step: (a) original image captured with a Canon XTI Rebel 6i camera, (b) L-channel from the L*a*b space, and (c) coherence filter applied to the L-channel; sample area 500 mm × 500 mm.
Figure 4. Preliminary filtering step: (a) original image captured with a Canon XTI Rebel 6i camera, (b) L-channel from the L*a*b space, and (c) coherence filter applied to the L-channel; sample area 500 mm × 500 mm.
Sensors 21 07534 g004
Figure 5. (a) Vector w+ in the original image, (b) the maximum index for Equation (2), (c) the coherence filter output, and (d) binary segmentation; sample area 500 mm × 500 mm.
Figure 5. (a) Vector w+ in the original image, (b) the maximum index for Equation (2), (c) the coherence filter output, and (d) binary segmentation; sample area 500 mm × 500 mm.
Sensors 21 07534 g005
Figure 6. (a) A section of the binary image, (b) the Euclidean distance transform on the binary image (in pseudo-color), (c) application of the top-hat filter to the distance transform, and (d) topological skeleton applied to the top-hat image (zoom area 12.8 mm × 12.8 mm).
Figure 6. (a) A section of the binary image, (b) the Euclidean distance transform on the binary image (in pseudo-color), (c) application of the top-hat filter to the distance transform, and (d) topological skeleton applied to the top-hat image (zoom area 12.8 mm × 12.8 mm).
Sensors 21 07534 g006
Figure 7. Central points after the neighbor spacing stabilization process; central point area equal to 28.3 mm × 28.4 mm.
Figure 7. Central points after the neighbor spacing stabilization process; central point area equal to 28.3 mm × 28.4 mm.
Sensors 21 07534 g007
Figure 8. (a) Sub-region and line profile superimposed on the angle normal to the crack, (b) intensity levels of the profile, and (c) width estimation by clustering.
Figure 8. (a) Sub-region and line profile superimposed on the angle normal to the crack, (b) intensity levels of the profile, and (c) width estimation by clustering.
Sensors 21 07534 g008
Figure 9. Plain mixture specimen with optimum calculation of the crack width and optimum angle of orientation according to the curvature.
Figure 9. Plain mixture specimen with optimum calculation of the crack width and optimum angle of orientation according to the curvature.
Sensors 21 07534 g009
Figure 10. Control point spacing according to the k-means parameter k: (a) k = 10, dmin = 8.5, davg = 3.238 mm, (b) k = 15, dmin = 13.5, davg = 3.236, (c) k = 20, dmin = 16.5, davg = 3.364, and (d) k = 40, dmin = 38, davg = 3.158.
Figure 10. Control point spacing according to the k-means parameter k: (a) k = 10, dmin = 8.5, davg = 3.238 mm, (b) k = 15, dmin = 13.5, davg = 3.236, (c) k = 20, dmin = 16.5, davg = 3.364, and (d) k = 40, dmin = 38, davg = 3.158.
Sensors 21 07534 g010
Figure 11. Dataset of samples generated from ECMs reinforced with different fiber types and a plain ECM. Samples were taken with two light sources: natural and artificial (area under analysis 500 mm × 500 mm).
Figure 11. Dataset of samples generated from ECMs reinforced with different fiber types and a plain ECM. Samples were taken with two light sources: natural and artificial (area under analysis 500 mm × 500 mm).
Sensors 21 07534 g011
Figure 12. Results of the proposed method on the dataset of Figure 11 (area under analysis 500 mm × 500 mm).
Figure 12. Results of the proposed method on the dataset of Figure 11 (area under analysis 500 mm × 500 mm).
Sensors 21 07534 g012
Figure 13. Zoom on the analyzed images (a), zoom area: 8.3 mm × 8.4 mm), (b), zoom area: 61.5 mm × 61.5 mm), and (c), zoom area: 51.3 mm × 51.3 mm).
Figure 13. Zoom on the analyzed images (a), zoom area: 8.3 mm × 8.4 mm), (b), zoom area: 61.5 mm × 61.5 mm), and (c), zoom area: 51.3 mm × 51.3 mm).
Sensors 21 07534 g013
Figure 14. Frequency distribution of the measurements performed automatically.
Figure 14. Frequency distribution of the measurements performed automatically.
Sensors 21 07534 g014
Figure 15. Boxplot of the measurements performed automatically.
Figure 15. Boxplot of the measurements performed automatically.
Sensors 21 07534 g015
Table 1. The method’s performance applied automatically and manually for 30 random points.
Table 1. The method’s performance applied automatically and manually for 30 random points.
Automatic MeasurementManual Measurement
(30 Points)
t-Test
Comparison
Image Features
Image CodeCountMean
(Pixel)
std
(Pixels)
Mean
(Pixels)
std
(Pixels)
Z-Scorep-ValueLight TypekFiber Type
V/1305.351.995.261.44−0.250.801Midday sun20Vegetal
V/2309.904.959.552.04−0.430.671Midday sun20Vegetal
V/3106.421.436.562.080.200.842Light (3000 K)10Vegetal
V/454.480.544.321.11−0.500.624Light (3000 K)20Vegetal
S/1397.374.726.214.65−1.030.308Afternoon Sun20Industrial
A/1757.282.407.801.870.880.386Afternoon Sun20Animal
A/26067.613.727.801.870.110.906Afternoon Sun10Animal
A/3958.383.268.833.060.490.622Light (3000 K)20Animal
A/419223.159.3922.718.22−0.310.761Midday sun10Animal
A/523417.867.1517.677.79−0.380.699Light (3000 K)10Animal
NoFiber/19512.436.3612.477.39−0.390.694Afternoon Sun20No Fiber
NoFiber/230814.376.4214.838.140.220.827Light (3000 K)20No Fiber
TOTAL1719 regions
Table 2. Comparison with other techniques and their main features.
Table 2. Comparison with other techniques and their main features.
TechniqueNoise RemovalLength
Estimation
Route Tracing
(Tortuosity)
Light
Source Setup
Spatial Width DistributionRef
Edge distance + Linear fitNone (automatic threshold)NoNoYesNo[15]
Percolation + BinarizationPercolation processingNoNoNoNo[16]
Percolation + Neighbor boundaryPercolation processingNoNoNoNo[22]
Skeletonization between two pointsNone (clumping process)NoYesNoNo[23]
Digital image correlation (DIC)None (automatic threshold)YesNoYesNo[2]
Top-Hat + Otsu BinarizationGaussian function-based spatial filterNoNoNoNo[5]
Feature extraction and SVM algorithmSteerable FilterYesNoNoNo[4]
Genetic AlgorithmMulti-sequential image filterYesYesNoNo[13]
Deep LearningFast-RCNN + TuFFNoNoNoNo[29]
Hessian structure propagationNoneNoYesNoNo[27]
Deep LearningYOLOYesNoYesNo[33]
filtering + edge searching.Frangi filteringYesNoNoNo[26]
Bwdist transform + Arc LengthMorphological operations (Aperture)YesYesYesNo[24]
M2GLDMin-Max Gray Level DiscriminationNoNoNoNo[25]
Proposed techniqueL*a*b + Coherence FilterYesYesNoYes
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Carrasco, M.; Araya-Letelier, G.; Velázquez, R.; Visconti, P. Image-Based Automated Width Measurement of Surface Cracking. Sensors 2021, 21, 7534. https://doi.org/10.3390/s21227534

AMA Style

Carrasco M, Araya-Letelier G, Velázquez R, Visconti P. Image-Based Automated Width Measurement of Surface Cracking. Sensors. 2021; 21(22):7534. https://doi.org/10.3390/s21227534

Chicago/Turabian Style

Carrasco, Miguel, Gerardo Araya-Letelier, Ramiro Velázquez, and Paolo Visconti. 2021. "Image-Based Automated Width Measurement of Surface Cracking" Sensors 21, no. 22: 7534. https://doi.org/10.3390/s21227534

APA Style

Carrasco, M., Araya-Letelier, G., Velázquez, R., & Visconti, P. (2021). Image-Based Automated Width Measurement of Surface Cracking. Sensors, 21(22), 7534. https://doi.org/10.3390/s21227534

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop