Double Low-Rank and Sparse Decomposition for Surface Defect Segmentation of Steel Sheet

Featured Application: The proposed DLRSD-based segmentation method can be applied for other industrial products, such as glass, fabric, LCD and AMOLED. Abstract: Surface defect segmentation supports real-time surface defect detection system of steel sheet by reducing redundant information and highlighting the critical defect regions for high-level image understanding. Existing defect segmentation methods usually lack adaptiveness to different shape, size and scale of the defect object. Based on the observation that the defective area can be regarded as the salient part of image, a saliency detection model using double low-rank and sparse decomposition (DLRSD) is proposed for surface defect segmentation. The proposed method adopts a low-rank assumption which characterizes the defective sub-regions and defect-free background sub-regions respectively. In addition, DLRSD model uses sparse constrains for background sub-regions so as to improve the robustness to noise and uneven illumination simultaneously. Then the Laplacian regularization among spatially adjacent sub-regions is incorporated into the DLRSD model in order to uniformly highlight the defect object. Our proposed DLRSD-based segmentation method consists of three steps: ﬁrstly, using DLRSD model to obtain the defect foreground image; then, enhancing the foreground image to establish the good foundation for segmentation; ﬁnally, the Otsu’s method is used to choose an optimal threshold automatically for segmentation. Experimental results demonstrate that the proposed method outperforms state-of-the-art approaches in terms of both subjective and objective tests. Meanwhile, the proposed method is applicable to industrial detection with limited computational resources.


Introduction
Surface defect detection plays an important role in quality enhancement in industrial product manufacturing. However, traditional defect detection is performed by human eyes, which yields low efficiency and high missing rate. Currently, vision-based automated defect detection has drawn much attention, which has important theoretical and practical value [1][2][3][4]. In automatic surface inspection of steel sheet, segmentation of surface defect is a significant step, which generates a binary map to identify defects. In the past two decades, commonly-used segmentation methods can be classified into three categories: statistical-based methods, filter-based methods and model-based methods. Statistical-based methods, such as Otsu's method [5], gray level co-occurrence matrix, local binary pattern, maximum entropy, region growing and morphological watersheds, are used to evaluate the spatial distribution of pixel intensities for segmentation. Filter-based methods, such as discrete Fourier discrete Fourier transform [6], discrete Gabor transform [7] and discrete wavelet transform [8,9], apply a bank of filters to the image, in which the energies of the filters response are utilized as features to segment the defects. Model-based approaches obtain certain models with specific feature distributions or other attributes using diverse descriptors, for instance, level set, fuzzy theory, partial differential equations and texture patterns.
Most recently, with the development of saliency detection technology, segmentation methods that use saliency map are gradually rising in the industrial defect inspection field. This method constructs a saliency map that highlights the defect regions standing out from the rest of the image, which provide the good foundation for segmentation. Guan et al. [10] proposed saliency map construction method using Gaussian pyramid decomposition. Then segmentation is conducted with the saliency map. This model exhibits good performance for strip steel defect detection. Li et al. [11] devised a low-rank representation-based saliency detection model for textile fabric defect detection. Zhao et al. [12] also presented a novel saliency detection model, which obviously improve the accuracy of automated defect segmentation.
These methods achieve good results on defect segmentation for a certain and homogeneous texture, but remain a challenging issue for segmentation with miscellaneous textures due to random disturbance. Specially, as the surface defect image of steel sheet has a low signal-to-noise ratio, low contrast between defect object and background, heterogeneous and scattered defect, cluttered and complicated background, these methods still lack of accuracy and suffer from limited adaptability and robustness in industrial practice.
Usually, a defect-free surface in industrial products has consistent texture. The emergence of defects can be regarded as the foreground object superposed in the regular-texture background. As shown in Figure 1, a surface defect image of steel sheet is decomposed into two parts: relatively homogeneous background image and a defect foreground image that is the desired image for the following segmentation. Inspired by the above analysis, an easy-to-implement method based on double low-rank and sparse decomposition (DLRSD) is proposed in this paper for surface defect segmentation. Considering double low-rank and sparse characteristics of surface defect image, combined with a local consistency constrain among spatially adjacent sub-regions by imposing Laplacian regularization, the feature matrix that form by can be adaptively decomposed into foreground feature matrix that form by defect foreground image and background feature matrix that form by background image in a certain feature space, respectively. Specifically, the foreground image is served as the source image for segmentation, which can better cope with the intra-class variations and background clutters, leading to a higher performance. Theoretical analysis and experimental results demonstrate the feasibility and effectiveness of the proposed DLRSD-based segmentation method for the surface defect of steel sheet. At the same time, it provides an interesting perspective for the industrial product's surface defect segmentation.
The rest of this paper is organized as follows. In Section 2, we review some existing saliency detection methods, especially the structural matrix decomposition-based methods. In Section 3 we introduce the proposed DLRSD model, including formulation and optimization. Section 4 presents the DLRSD-based defect segmentation method. Also, we give more detail on enhancing the original defect foreground image. Section 5 describes experimental results between our proposed method and some state-of-the-art methods. Finally, conclusions are given in Section 6. Inspired by the above analysis, an easy-to-implement method based on double low-rank and sparse decomposition (DLRSD) is proposed in this paper for surface defect segmentation. Considering double low-rank and sparse characteristics of surface defect image, combined with a local consistency constrain among spatially adjacent sub-regions by imposing Laplacian regularization, the feature matrix that form by I can be adaptively decomposed into foreground feature matrix that form by defect foreground image F and background feature matrix that form by background image B in a certain feature space, respectively. Specifically, the foreground image F is served as the source image for segmentation, which can better cope with the intra-class variations and background clutters, leading to a higher performance. Theoretical analysis and experimental results demonstrate the feasibility and effectiveness of the proposed DLRSD-based segmentation method for the surface defect of steel sheet. At the same time, it provides an interesting perspective for the industrial product's surface defect segmentation.
The rest of this paper is organized as follows. In Section 2, we review some existing saliency detection methods, especially the structural matrix decomposition-based methods. In Section 3 we introduce the proposed DLRSD model, including formulation and optimization. Section 4 presents the DLRSD-based defect segmentation method. Also, we give more detail on enhancing the original defect foreground image. Section 5 describes experimental results between our proposed method and some state-of-the-art methods. Finally, conclusions are given in Section 6.

Related Work
During the past few years, there are many methods attempting to segment the salient object from the saliency map of an input image [13][14][15][16][17][18][19]. The quality and effectiveness of segmentation are decided by the quality of saliency map [20]. Based on the milestone work, some structural matrix decomposition-based methods transform a saliency detection problem into a feature subspace decomposition problem, which can improve detection results in terms of both speed and accuracy. Particularly, many studies conclude that low-rank matrix decomposition-based methods can obtain better saliency detection performance. These methods assume that an image can be represented as a combination of a highly redundant part (e.g., visually consistent background regions) and a sparse part (e.g., salient object foreground regions). Therefore, given the feature matrix of an input image, it can be decomposed into a low-rank matrix corresponding to the non-salient background and a sparse matrix corresponding to the salient foreground objects. Yan et al. [21] employed sparse coding as a feature representation vector of image. Zou et al. [22] designed multi-scale superpixel segmentation to construct the feature matrix and prior matrix. Although Shen et al. [23] adopted learnt linear transformation of the feature space to integrate low-level features and high-level prior knowledge, the learnt transform matrix is correlated for training data set. Unfortunately, the sparsity assumption of the salient objects can't be guaranteed universally, especially when the salient objects with big size occupy most of the image, and then suffer from limited adaptability. Therefore, Peng et al. [24] developed tree-structured sparsity-inducing regularization and Laplacian regularization to disentangle the salient objects and background precisely, and then obtained competitive results. But it may be difficult to suppress some small background regions with distinctive appearances because of the constructed index-tree is not precise enough. Subsequently, Sun et al. [25] presented diversity-induced regularization based on Hilbert-Schmidt independence criterion, which make the background much cleaner in the saliency map and boost the saliency detection performance. But, they don't consider the low-rank characteristic for the foreground regions and background regions simultaneously, and ignore the spatial and pattern relations of image regions, which may lead to very noisy saliency map and influences on the final segmentation performance.
To solve the problems mentioned above, the proposed DLRSD model considers the correlation between defective regions and defect-free regions, which is different from existing methods in essence. Besides, it uses the nuclear norm to depict the low-rank property of defect object rather than consider it as the sparse noises, which can produce more accurate and reliable saliency map that represents the defect foreground image.

Double Low-Rank and Sparse Decomposition Model
In this section, we will introduce the proposed DLRSD model and optimization procedure in details.

Problem Formulation
Let {R 1 , R 2 , · · · , R K } be a set of K non-overlapping sub-regions of a surface defect image I, all the feature vectors of sub-regions can construct the feature matrix D. The proposed DLRSD model is to design an effective model to decompose the feature matrix D into a feature matrix S that represents a defect foreground image F and a feature matrix L that represents a background image B: In order to separate defect regions and background regions accurately, some constrains are needed for characterizing two feature matrices S and L. According to the surface defect image I that is pre-processed by superpixel segmentation, both defect foreground and background contain multiple homogeneous and highly similar sub-regions, for each defective sub-region, the corresponding locations in saliency map has high probability in larger brightness, indicating that it has higher saliency value. Besides, different defective sub-regions are highly correlated and the corresponding feature vectors lie in a low-dimensional subspace. Therefore, the feature matrix S is expected to be low-rank. Meanwhile, most of background sub-regions tend to have lower saliency value. They are strongly correlated and lie in a low-dimensional feature subspace that is independent of the defect foreground subspace. The strong correlations among the background sub-regions suggest that feature matrix L may have the low-rank property. What is more, in order to reduce the influence of noises and enhance the robustness to uneven illumination, we assume that the background lies in a sparse feature subspace and can be characterized by a sparse matrix.
Based on above analysis, the structured matrix decomposition model can be constructed as follows: min where rank(·) denotes the rank of matrix; ||·|| 0 denotes l 0 norm of matrix, which equals the number of non-zero element of matrix; Θ(S, L) denotes the regularization to enlarge the margin and reduce the coherence between the feature subspaces induced by S and L; D ∈ R d×K represents the feature matrix; α > 0, β > 0 and γ > 0 are regularization parameters.
To separate the defect object from the background easily, spatially adjacent sub-regions with smaller spatial distance and more similar feature vector should be assigned to similar and higher weight values, the local invariance assumption [26] based Laplacian regularization Θ(S, L) [24] can be defined as follows: where tr(·) denotes the trace of a matrix; s i denotes the i-th column of matrix S; the element w ij of affinity matrix W ∈ R K×K denotes the weight that represents the feature similarity between sub-regions R i and R j ; M ∈ R K×K is a Laplacian matrix. According to the undirected graph model from a surface defect image, each sub-region is represented by a node, the affinity matrix W is where p i ∈ R 2 and p j ∈ R 2 denote the central coordinate of R i and R j ; represents spatial connectivity between R i and R j , which represents the spatial contiguity; exp gives the feature similarity between R i and R j ; σ p and σ f are two scalars. The Laplacian matrix M is In particular, the Laplacian regularization Θ(S, L) can preserve the local consistency and invariance among the spatially adjacent sub-regions with similar saliency values in saliency maps. More specifically, the defect foreground is more uniformly highlighted and the background noise is also better suppressed, and eventually separates the defect from the background as much as possible.

Optimization
As rank(·) and ||·|| 0 are not convex, Equation (2) is NP-hard problem. A common heuristic criterion is to replace rank(·) and ||·|| 0 are replaced by nuclear norm ||·|| * and l 1 norm ||·|| 1 respectively. It has been shown that nuclear norm-based models can obtain the optimal low-rank solution in many kinds of applications [27,28]. Then, Equation (2) can be converted to the following convex surrogate optimization problem: where ||·|| * equals the sum of singular values of matrix; ||·|| 1 equals the sum of the absolute values of To solve Equation (6) efficiently, the alternating direction method (ADM) algorithm [29] can be adopted. By introducing the auxiliary variables H and J, the augmented Lagrange function is given as follows: where ||·|| 2 F denotes the Frobenius norm of matrix, which is defined as the sum of squares of each element of matrix; Y 1 , Y 2 and Y 3 are Lagrange multipliers; µ > 0 is a penalty parameter.
Therefore, Equation (7) can be converted to the following equivalent optimization problem: The above optimization problem can be solved by alternately updating one variable while others fixed. The detailed ADM algorithm for proposed DLRSD model is summarized in Algorithm 1.
(1) Updating H In order to solve H, the optimal solution can be obtained by Equation (9): Differentiating it with respect to H, and let it to be zero, therefore The close-form solution can be obtained as follows: (2) Updating J In order to solve J, the optimal solution can be obtained by Equation (12): The solution is where Ψ γ µ (·) denotes soft-thresholding shrinkage operator, which is defined as where T denotes a matrix, T ij denotes the (i, j)-th element of T, sgn(T) is the matrix whose entries are the signs of those of T.
(3) Updating L In order to solve L, the optimal solution can be obtained by Equation (15): It can be rewritten as follows: Its solution is , svd(·) denotes singular value decomposition operator.
(4) Updating S In order to solve S, the optimal solution can be obtained by Equation (18): It can be rewritten as follows: Its solution is where (U, Σ, V) = svd 1 2 where 0 < ρ < 1.

DLRSD-Based Surface Defect Segmentation
In this section, we describe how to apply the proposed DLRSD model to surface defect segmentation. The segmentation method has three stages. In first stage, we use DLRSD model to obtain the defect foreground image F. While in second stage, we utilize regression optimization to enhance F. At last, the segmentation is finished by Otsu's method. The framework of DLRSD-based segmentation method is shown in Figure 2, the detailed procedure is summarized in Algorithm 2.
achieve the better structural information about defect image, we conduct superpixel segmentation for image by adaptive simple linear iterative clustering (ASLIC) algorithm [30]. Each compact, edge-aware and perceptually homogeneous sub-region , ,…, can be represented by feature vector ∈ ℝ , where represents the mean feature vector of all pixels that belong to , where denotes the number of sub-regions. By arranging into a matrix, the feature matrix = , , ⋯ , ∈ ℝ × of image is obtained. Step

Matrix Decomposition
According to Algorithm 1, the input feature matrix is decomposed into structured components and . According to the obtained * = , , ⋯ , ∈ ℝ × and * = , , ⋯ , ∈ ℝ × , each column of these two matrixes represents the feature vector of corresponding sub-region, respectively. Then, we transfer * and * from the feature domain to the spatial domain for constructing saliency map. The saliency value of each sub-region in foreground image and background image are max and max , respectively, where ∈ ℝ × and ∈ ℝ × denotes the -th column of * and * , max • denotes the maximum component of the vector, = 1, 2, … , . After allocating the saliency value to corresponding pixels and normalizing, the defect foreground image and background image can be obtained.

Enhancement
As shown in Figure 2, the original foreground image can be enhanced in consistency, completeness of defect objects and suppression of background noise. In the paper, the regression optimization method is adopted by combining foreground image and background image . The optimization problem can be formulated as follows:

Feature Matrix Construction
According to [23,24], for each pixel {I i } i=1,2,...,N of a surface defect image I, where N denotes the number of pixels, different types of low-level visual features, including gray-scale, Gabor filters and steerable pyramids, are extracted.
(1) Gray-scale The pixel value of each pixel in defect image I is extracted for gray-scale feature, which is normalized by subtracting its mean value over the entire image.  All those 25 features are then stacked vertically to construct a 25-dimension feature vector {f i } i=1,2,...,N ∈ R d for each pixel. Then, in order to improve the efficiency of defect detection and achieve the better structural information about defect image, we conduct superpixel segmentation for image I by adaptive simple linear iterative clustering (ASLIC) algorithm [30]. Each compact, edge-aware and perceptually homogeneous sub-region R j j=1,2,...,K can be represented by feature vector f j ∈ R d , where f j represents the mean feature vector of all pixels that belong to R j , where K denotes the number of sub-regions. By arranging f j into a matrix, the feature matrix D = f 1 , f 2 , · · · , f K ∈ R d×K of image I is obtained.

Matrix Decomposition
According to Algorithm 1, the input feature matrix D is decomposed into structured components S and L. According to the obtained S * = (s 1 , s 2 , · · · , s K ) ∈ R d×K and L * = (l 1 , l 2 , · · · , l K ) ∈ R d×K , each column of these two matrixes represents the feature vector of corresponding sub-region, respectively. Then, we transfer S * and L * from the feature domain to the spatial domain for constructing saliency map. The saliency value of each sub-region in foreground image F and background image B are max s j and max l j , respectively, where s j ∈ R d×1 and l j ∈ R d×1 denotes the j-th column of S * and L * , max(·) denotes the maximum component of the vector, j = 1, 2, . . . , K. After allocating the saliency value to corresponding pixels and normalizing, the defect foreground image F and background image B can be obtained.

Enhancement
As shown in Figure 2, the original foreground image F can be enhanced in consistency, completeness of defect objects and suppression of background noise. In the paper, the regression optimization method is adopted by combining foreground image F and background image B. The optimization problem can be formulated as follows: where w f i denotes saliency value of sub-region in foreground image F, w f i = Val P j ; w b i denotes saliency value of sub-region in background image B, w b i = Val Q j ; s i denotes the optimized saliency value of sub-region in foreground image F. (23) can be reformulated as follows: where 1 ∈ R K×1 denotes a one vector, M ∈ R K×K denotes the same Laplacian matrix in Equation (5).
Differentiating it with respect to s, and let it to be zero, therefore The solution is Through Equation (26), the sub-regions within the same class (foreground or background) have more similar saliency values while the sub-regions from different classes (foreground and background) have different saliency values. The saliency value of defect sub-region in foreground image is bigger, while the saliency value of background sub-region is smaller, so that the surface defect object can be highlighted further.

Segmentation
After obtaining the enhanced foreground image F, the high-quality binary image can be obtained through a simple Otsu's method. In binary image of surface defect, white pixel represents surface defect regions, and black pixel represents background regions.

Algorithm 2 DLRSD-based defect segmentation.
Input: Surface defect image I Output: Binary segmentation image 1: Construct the feature matrix D of I 2: Run Algorithm 1 to get the defect foreground feature matrix S 3: Enhance defect foreground image F 4: Segment enhanced F by Otsu's method

Experiment
In this section, several experiments are conducted to verify the superiority of our proposed method. We first introduce the experimental setups, which include parameters settings and evaluation metrics. Then, computational complexity, convergence, noise immunity and segmentation results are discussed. At last, the qualitative and quantitative comparisons are presented.

Experimental Setup
In order to verify and evaluate the effectiveness and robustness of the proposed method, we have adopted the NEU surface defect database established by Kechen Song [12] in our experiments. The size of each surface defect image is 200 × 200 and the number of image is 300 per class. Two typical surface defect images, such as Patch and Scratch, are selected in the experiments. Our proposed method is compared with eight representative saliency detection methods quantitatively and qualitatively, such as RPCA [28], IS [13], ULR [23], RBD [31], SBD [32], DSR [33], RS [16] and SMF [24], where RPCA, IS, ULR, RBD, SBD, DSR, RS and SMF represent the method of robust principal component analysis, image signature, unified low rank matrix recovery, robust background detection, spaces of background-based distribution, dense and sparse reconstruction, ranking saliency and structured matrix decomposition, respectively. Only a few examples are shown in the paper, the whole segmentation results are uploaded in Baidu Disk (https://pan.baidu.com/s/1QkwFfWsUE9hKL86prlL4nw, Code: iydw).

Parameters Settings
In Equation (6), α represents the redundancy of defect foreground, β represents the uniformity of defect foreground, γ represents the sparsity of background. We conduct some experiments to study the detection performance variation with respect to different α, β and γ, which shows that the detection performance can achieve a high level at α ∈ (0.2, 0.4), β ∈ (0.9, 1.3) and γ ∈ (0.05, 0.25). In order to achieve the better segmentation results, α, β and γ are set to 0.35, 1.2 and 0.1, respectively. For other methods in our comparison, we use the source codes provided by the authors with default parameters.

Evaluation Metrics
The qualitative evaluation metrics refers to evaluate the detection performance based on human subjective feeling. For example, the boundary of surface defect is clear, and the contrast between defect object and background is obvious.
There are five quantitative evaluation metrics, including precision-recall (P-R) curve, receiver operating characteristic (ROC) curve, average F-Measure (F ζ ), area under ROC (AUC) and mean square error (MAE). They are defined as follows: where a pixel that belonging to defect is defined as a positive example, and a pixel that belonging to background is defined as a negative example; true positive (TP) indicates that the positive pixel is judged correctly, true negative (TN) indicates that the negative pixel is judged correctly, false positive (FP) indicates that the positive pixel is judged as the negative pixel mistakenly, false negative (FN) indicates that the negative pixel is judged as the positive pixel mistakenly; precision = TP/(TP + FP), recall = TP/(TP + FN); N represents the number of surface defect image samples of the same class, H and W denotes the height and width of surface defect image, respectively; precision is defined as the percentage of defect pixels correctly assigned, while recall is the ratio of correctly detected defect pixels to all true defect pixels. F ζ represents the weighted harmonic mean of precision and recall. Besides, P-R curve is obtained by binarizing the saliency map using a number of thresholds ranging from 0 to 255; TPR represents true positive rate, FPR represents false positive rate; MAE measures the dissimilarity between the saliency map S and the ground truth G.

Analysis of Computational Complexity
According to Algorithm 1, the main computational load is singular value decomposition operation in updating matrix S and L. As the size of matrix D is d × K, the computational complexity is reduced from dK 2 to (drK) by the low-rank constraint, where r denotes the rank of matrix. In our experiments, d = 25, K = 100, so the computational complexity is low.

Analysis of Convergence
According to Algorithm 1 and ADM algorithm, when penalty parameter sequence {µ k } is increasing monotonically and bounded, the Lagrange multipliers Y 1 , Y 2 and Y 3 can converge to the optimal solution linearly; when {µ k } is increasing monotonically and unbounded, Y 1 , Y 2 and Y 3 can converge to the optimal solution super-linearly. As shown in Figure 3, the x-axis denotes the iteration number, and the y-axis is the value of objective function. We can see that the objective function value converges in a very fast manner, usually within 40 iterations, which also proves the fast convergence property of the proposed DLRSD model.

Analysis of Segmentation Results
From enhanced defect foreground image shown in Figure 4c, it has achieved the goal of "highlight the foreground and suppressing the background". It can accurately extract the entire defect object and assigns nearly uniform saliency values to all sub-regions within the defect objects. Figure 4d shows that the segmentation images are similar to ground truth, the whole defect object can be uniformly highlighted, and boundary of defect object is well-defined. Therefore, we locate the defects accurately.

Analysis of Segmentation Results
From enhanced defect foreground image shown in Figure 4c, it has achieved the goal of "highlight the foreground and suppressing the background". It can accurately extract the entire defect object and assigns nearly uniform saliency values to all sub-regions within the defect objects. Figure 4d shows that the segmentation images are similar to ground truth, the whole defect object can be uniformly highlighted, and boundary of defect object is well-defined. Therefore, we locate the defects accurately.

Analysis of Robustness to Noise
Considering the surface defect image is polluted by Gaussian noise with SNR, including 22 dB, 18 dB, 14 dB and 10 dB, the same experiments are conducted to verify the robustness of the proposed DLRSD model. According to Table 1, when SNR decreases gradually, the AUC and MAE can remain a high level, especially when SNR = 18 dB, AUC can remain around 0.8. It's shown that the proposed DLRSD model is robust to noise and can lead to better saliency detection result, which establishes the good foundation for segmentation. The experimental results also indicate that adding sparse constraint for background can reduce the influence from noises, which is a reasonable strategy for surface defect detection.

Qualitative Comparison
The qualitative comparison results by the proposed method and other eight methods are shown in Figure 5. It's shown that most saliency detection methods can handle well simple images with relatively homogenous background (e.g., row 4, 5, 7 and 8). They can uniformly highlight the whole defect object and generate high-quality saliency map and segmentation image. However, for some complex defect images containing multiple objects (e.g., row 5, 6, 10 and 11), having a cluttered background (e.g., row 6), and showing there are similarities between the defect objects and background (e.g., row 2 and 9), the whole defect objects could not be uniformly highlighted, and parts of the background being falsely taken as the defect objects. It can be seen that the contrast of saliency maps obtained by RPCA, DSR and RS is low and ambiguous, especially for Patch defects (e.g., row 5 and 6), which is difficult to define a proper threshold to segment the defects. The saliency maps obtained by RBD and SBD miss detecting parts of the defect objects, while some incorrectly include background regions into detection results. Hence, there are some missing defects and fake defects in their final segmentation image. Differently, although IS, ULR and SMF produce the good saliency map, there are many pixels that belonging to the background are misjudged by defect, and some background regions also stand out with the defect regions. By contrast, our proposed method separates the defect objects from the background successfully and locates various defects precisely. It more efficiently highlights the complete defect object with well-defined boundaries and effectively suppresses the backgrounds than the other saliency detection methods. These results illustrate our proposed method not only enhances the contrast between surface defect and background effectively but also improves the robustness to the different illumination conditions, various shapes, scales, directions and locations of surface defect.   Figure 6 shows the quantitative results of the proposed DLRSD model against eight state-of-the-art methods. It is known that it perform competitively and is both better than the other methods in terms of the P-R curve, ROC curve and F-Measure curve. Especially, the precision can remain above 90% within a large threshold range.  Figure 6 shows the quantitative results of the proposed DLRSD model against eight state-of-theart methods. It is known that it perform competitively and is both better than the other methods in terms of the P-R curve, ROC curve and F-Measure curve. Especially, the precision can remain above 90% within a large threshold range.  Table 2 summarizes the quantitative results of all the eight methods. We can see that the proposed DLRSD model has achieved the best performance in AUC, Fβ and MAE. Compared with SMF, it increased by 8.52% and 4.05% in AUC and F ζ , respectively, decreased by 5.01% in MAE. All experiments are run in Matlab 2018a on a PC with an Intel Core i7-4790@2.90GHz CPU and 8GB RAM, the running time of the proposed DLRSD model is slightly slower than RS but much faster than ULR and SMF. Based on the above qualitative and quantitative analyses, it confirms that our proposed method consistently outperforms some state-of-the-art methods and verifies the effectiveness of the proposed structural constraints in separating the low-rank and sparse subspaces.

Conclusions
Based on the salient characteristics of the defects in the surface defect image of steel sheet, we formulate the defect segmentation as a problem of saliency detection. We design a double low-rank and sparse decomposition model to obtain high-quality defect foreground image directly, which provides a robust way to segment the surface defect. We experimentally compare our proposed method with some state-of-the-art methods on surface defect images. The experimental results prove  Table 2 summarizes the quantitative results of all the eight methods. We can see that the proposed DLRSD model has achieved the best performance in AUC, F β and MAE. Compared with SMF, it increased by 8.52% and 4.05% in AUC and F ζ , respectively, decreased by 5.01% in MAE. All experiments are run in Matlab 2018a on a PC with an Intel Core i7-4790@2.90GHz CPU and 8GB RAM, the running time of the proposed DLRSD model is slightly slower than RS but much faster than ULR and SMF.  [13] 0.7140 0.2814 0.2485 0.0032 ULR [23] 0.7843 0.4780 0.2976 4.5504 RBD [31] 0.7125 0.4607 0.2090 0.0331 SBD [32] 0.6907 0.5038 0.2390 0.6619 DSR [33] 0.7786 0.6264 0.1626 1.1797 RS [16] 0.7469 0.6454 0.1758 0.1281 SMF [24] 0 Based on the above qualitative and quantitative analyses, it confirms that our proposed method consistently outperforms some state-of-the-art methods and verifies the effectiveness of the proposed structural constraints in separating the low-rank and sparse subspaces.

Conclusions
Based on the salient characteristics of the defects in the surface defect image of steel sheet, we formulate the defect segmentation as a problem of saliency detection. We design a double low-rank and sparse decomposition model to obtain high-quality defect foreground image directly, which provides a robust way to segment the surface defect. We experimentally compare our proposed method with some state-of-the-art methods on surface defect images. The experimental results prove that the proposed method performs efficiently and competitively for the surface defect segmentation task and has a strong adaptive ability for the complex and varying surface defects of steel sheet. Our proposed method is an unsupervised framework, which skips the training process and therefore enjoys more flexibility. In the future, we will focus on combining our proposed method with convolutional auto-encoder and expanding the method to other industrial products' defect detection.