Structure Label Matrix Completion for PolSAR Image Classification

Terrain classification is a hot topic in polarimetric synthetic aperture radar (PolSAR) image interpretation that aims at assigning a label to every pixel and forms a label matrix for a PolSAR image. From the perspective of human interpretation, classification is not in terms of pixels, but decomposed perceptual groups and structures. Therefore, a new perspective of label matrix completion is proposed to treat the classification task as a label matrix inversion process. Firstly, a matrix completion framework is built to avoid processing the large-scale PolSAR image data for traditional feature and classifier learning. Secondly, a light network is used to obtain the known labels for the complete task by uniform down-sampling of the entire image that aims to keep the shape of regions in a PolSAR image and reduce the computational complexity. Finally, the zeroth- and first-order label information is proposed as the prior distribution for label matrix completion to keep the structure and texture. The experiments are tested on real PolSAR images, which demonstrates that the proposed method can realize excellent classification results with less computation time and labeled pixels. Classification results with different down-sampling rates in the light network also prove the robustness of this method.


Introduction
Polarimetric synthetic aperture radar (PolSAR) image analysis has attracted increasingly more attention for rich multi-polarization information compared with single-polarization synthetic aperture radar (SAR) [1]. PolSAR image terrain classification is a widespread application and is of great importance in city planning, sea monitoring, geological exploration, and plant-growth status assessment [2]. The classification task aims to tag every pixel a label corresponding to one definite terrain. Thus, the classification result is a label matrix with every element corresponding to the label of one pixel in a PolSAR image. In the last few years, numerous methods have been proposed to obtain the label matrix in supervised, unsupervised, and semi-supervised ways [3][4][5][6][7]. These methods usually pay attention to two aspects, namely, feature extraction and classification for the supervised way or clustering for the unsupervised way.

Related Works
In the case of feature acquisition, two alternative directions are followed according to researcher interest. First, it is well known that a PolSAR image is obtained by electromagnetic imaging. Such a

Motivation and Contributions
However, the above-listed PolSAR image classification methods must extract features and classify all the pixels in one image to obtain the label matrix. Supervised classification methods need labeled pixels to learn the models, and a better result is usually obtained with more labeled pixels. It is a tedious task for PolSAR images to label a significant number of pixels [21,43]. Furthermore, unsupervised and semi-supervised classification methods suffer from a huge amount of computation, especially for large-scale PolSAR images [44,45]. In addition, it may be difficult for the existing methods to meet the processing demand of massive data collected by the PolSAR systems. Thus, many classification methods were developed by considering the contextual information of the label matrix to alleviate the labor-intensive labeling. The representative way is an additional prior distribution with a Markov random field (MRF), and better classification results are obtained [11,21,44,46]. However, traditional methods with contextual information still need heavy computation for a large number of pixels in one PolSAR image [21,46].
From the perspective of human interpretation, classification is not in terms of pixels, but decomposed perceptual groups and structures. Because pixels in a consistent region belong to the same class with a high probability, feature and classifier learning for all the pixels in this region are redundant. From this point of view, in this paper realizing the PolSAR classification task as label matrix completion is proposed. Matrix completion is an inverse process to fill the missing terms with a few of the elements known in the matrix, which can be completed accurately by the prior distribution of all the elements. Then, it turns the classification task into matrix completion with a solid mathematical foundation method in theory as in [47]. The most important point in label matrix completion is how to obtain the known labels and the prior distribution of the entire label matrix for one PolSAR image. Regarding the known labels, it is well known that the matrix can be completed with as many elements known as possible. In this paper, we obtain the known labels from a uniform down-sampling of PolSAR images. In this way, the traditional classification approach of feature and classifier learning is performed on fewer pixels. Meanwhile, the uniform sampling also keeps the region of the original PolSAR image. Moreover, in order to retain the texture and structure information, the labels of the down-sampled image are transferred by zeroth-and first-order label prior distributions to the entire label matrix. The zeroth-order prior term retains the label in the corresponding location of the entire image. The first-order prior term describes the structure information of the region consistency and sharp boundary between different regions.
As a novel approach, in this paper, a label matrix completion method for PolSAR image classification is proposed, and the main contributions are the following.

•
First, label matrix completion is introduced for the PolSAR image classification for the first time to the best of our knowledge. Label matrix completion solves a matrix inverse problem to obtain the classification result instead of the map from PolSAR data to the label field. • Second, a uniform down-sampled PolSAR image is used to obtain the known labels for the entire label matrix. In this way, the heavy computational burden in semi-supervised and unsupervised classification methods is relieved. • Third, the zeroth-and first-order label prior distributions are proposed to complete the entire label matrix with the known labels of the sampled PolSAR image. Thus, the final label matrix can be obtained with fewer labels known than traditional classification methods by this label prior to distribution. This results in state-of-the-art classification performance on both accuracy value and computation time.
The rest of this paper is organized as follows. In Section 2, the proposed structure label matrix completion method for PolSAR image classification tasks is introduced in detail. In Section 3, the optimization process and analysis of the model and computation complexity are illustrated to demonstrate the practicability of the proposed method. The analysis of parameters and experimental results on PolSAR images are presented in Section 4 to validate the classification effectiveness. Conclusions are drawn and planned future work described in Section 5.

Proposed Structure Label Matrix Completion Method for PolSAR Image Classification
To describe the proposed method, the PolSAR classification task is derived into a label matrix completion problem as a matrix inverse equation first. Then, a uniform under-sampled PolSAR image is introduced to obtain the known labels. Finally, the entire label matrix is completed by the zeroth-and first-order label prior distributions with the known labels to retain the region and structure information.

From Classification to Label Matrix Completion
Generally, let X ∈ R h×w×d be a PolSAR image with h, w and d denoting its height, width, and dimension of every pixel, respectively. Above all, labeled pixels and their corresponding labels are selected randomly to form the subset {X L , Y L }, while the remaining pixels are unlabeled and denoted {X U , Y U }. The traditional PolSAR image classification task aims to assign each unlabeled pixel with a proper label to produce an estimated classification map matrix Y U .
The traditional classification method consists of PolSAR image feature transformation or learning, as well as classifier learning. With the learned determinate model, the classification map of unlabeled pixels will be predicted. This process is illustrated as follows: where Φ and Θ are the designed feature and classifier models, respectively. This process must extract features and classification for all the pixels in one PolSAR image with a large computational burden. Furthermore, classification performance for different PolSAR images may vary widely. In fact, it is not necessary to extract features and classification for every pixel in a homogeneous region because pixels in one region usually belong to the same terrain. In addition, redundant operations cost much computation time. Therefore, a new perspective is considered inspired by label matrix completion to achieve less computational overhead for PolSAR image classification. In the following, basic knowledge about matrix completion will be introduced first.
, the following is a standard: formulation for matrix completion [47].
where p(Y) is any information known for the recovered matrix. In addition, e > 0 and Dis(P Ω obs (Ŷ), P Ω obs (Y) is the distance loss used to measure the fitness of the recovered matrix Y to the observed matrixŶ. Note that, for the cause of noisy labels existing in the measurement matrixŶ, the above constraint allows for a slight discrepancy between the recovered and observed matrices.
For PolSAR image classification, observed matrixŶ is replaced by known labels of pixels filling its corresponding position. It can be noted that three terms should be designed to perform classification with the label matrix completion framework. Known pixels in the observed matrixŶ are obtained from a light classification network for a down-sampled PolSAR image. Then, the prior term of the label matrix is designed with texture and structure information. Finally, the square loss function is designed for distance measure Dis, resulting in the following restraint expression form: where F is the Frobenius norm ||a|| F = ∑ i a 2 i (square root of the sum of squares for all the elements in the matrix). In the following, the detailed description of the learning process of label prior p(Y) and observed label matrixŶ are introduced.
In this paper, the sampling method is introduced into the input domain of a classification model to alleviate the labor-intensive labeling and heavy computational burden for highly redundant pixels. Specifically, a smaller imageX ∈ Rĥ ×ŵ×d of the entire PolSAR image X is obtained by the uniform down-sampling method as shown in Figure 1. To reduce the computational burden, the smaller image is used for classification to obtain a label matrixŶ. Thus, the overall pixels are much fewer in number than in the traditional classification method that extracts features and learns parameters from the entire image. Then, the classification result of this down-sampled image is transferred to the label matrix Y for the entire image. On the one hand, the down-sampled image needs fewer labeled pixels and less computation time to train the model. On the other hand, the region structure information will be retained in the sampled image of the entire image. Thus, it is desirable to design the loss so as to transfer all label information without missing any important ones from the known labels.

Completion of Label Matrix by Matrix Completion
Label matrix completion requires two important functions: obtaining the known labels and understanding the prior distribution of the entire label matrix. As shown in the entire framework depicted in Figure 2, the light network processes a uniform down-sampled PolSAR image as in Section 2.1 to obtain the known labels. Then, the prediction result is used to complete the entire label matrix with prior distribution from the zeroth-and first-order perspectives to keep the texture and structure information. The specific form for the zeroth-and first-order information is introduced in the following. As shown in the left-hand part of Figure 3, the squares in green color are the pixels classified with a light network for the uniform down-sampled matrix from the entire image, while the white squares are pixels to be completed by the known ones represented by the green squares. Because the green squares are transferred from the sampled image to the entire image in the corresponding index, it is named the zeroth-order prior.  Regarding the first-order prior distribution, the structure and contextual information of an entire PolSAR image is considered. In this paper, a two-structure prior is designed with the fact that pixels in the same region belong to the same class and the boundary between different regions is obvious. On the one hand, to make pixels in the same region have the same label, a spatially smooth term is improved by MRF priority. On the other hand, the boundary between different regions means the labels on the two sides of this boundary are different. Thus, a sharp edge is modeled as a Laplace prior as shown in the right-hand part of Figure 3. Specifically, we formulate the objective function for our label prior p(Y) learning as where Y d is the predicted label matrix for the uniform down-sampled PolSAR image with the ith vector as one-hot for pixel i. Z indicates the zeroth-order information about the priority of the entire image. This term aims to fill the entire label matrix Y with Y d in the corresponding location index of the entire image. The other elements are zero in the rest of the locations, which results in the label matrixŶ as in Equation (3). F is the first-order information transferred from the sampling image to the entire label matrix, which contains the spatial and contextual information for the label matrix. Specifically, we where N i indicates the eight-neighborhood index of pixel x and Y is the logits output of the corresponding network. The first term ensures the spatial consistency that elements in the neighbor location have the same label, while L(Y) is the Laplace variance that clarifies the edge. Here, this term is defined as a quadratic function of Y; therefore, a differentiable function can be easily integrated into the learning process. It is given by [48] L(Y) = 1 The light network in Figure 2 is operated on the uniform down-sampled PolSAR image, which is obtained from the entire image and results in less computational burden. Here, the semi-supervised method is used for combining the labeled and unlabeled pixels for a better classification result in the light network. With excellent classification performance, the method proposed for the imbalanced PolSAR images is used as shown in our preliminary work [24]. The cost-sensitive latent space learning network is built based on the parametric feature and classifier learning framework for imbalanced PolSAR images, where classifier and latent space learning are defined as optimizing the posterior and likelihood function for labeled pixels, respectively. On the one hand, the down-sampled PolSAR images have fewer samples in minor classes, which results in more imbalanced data distribution than in the entire image. On the other hand, the proposed classification method can obtain high accuracy and robust results for balanced and imbalanced PolSAR images.

Model Analysis
The label matrix completion task in Equation (3) consists of two crucial aspects: the known labels learning and matrix completion with label prior distribution. To keep the region and structure of the entire PolSAR image, the known labels are obtained by classifying a uniform down-sampled PolSAR image. The sampled image has the same structure as the entire PolSAR image, but there are fewer pixels to be processed. Notably, the classification method for this down-sampled image can be any supervised, unsupervised, or semi-supervised method with good results because it is just a way to obtain the known labels to complete the entire label matrix. The labels obtained from the uniform down-sampled image are the zeroth-order information to complete the entire label matrix. With first-order prior distribution information, a smooth and edge-shaped label matrix is obtained.
The proposed method in this paper is a new framework incorporating the innovative perspective of label matrix completion. Instead of designing a map function from the PolSAR data field to the label field [6,21,36,49], a label matrix inverse recovery method is used for the PolSAR image classification task to avoid the tedious learning processes of feature and classifier. The uniform down-sampling of the entire image may retain the structure and shape information for terrains in a PolSAR image compared with random sampling. Meanwhile, the information of the PolSAR data is used in this down-sampled image to consider the physical characteristic of PolSAR data. Most important, the classification model adopted for light lines can be any method that works effectively. This characteristic increases the generalization of the proposed model for different PolSAR images and makes the proposed method a new baseline for PolSAR image classification methods from the perspective of label matrix completion.

Optimization and Complexity Analysis
In this section, a detailed optimization scheme is presented to solve the proposed classification model. First, a form of the maximum a posteriori (MAP) conversion of the probability Equation (4) is obtained to update variables and parameters in this method. Then, the computational complexities of the proposed model for each step are discussed (Algorithm 1).

Optimization Process
To solve Equation (4), a two-stage process can be obtained by minimizing an error function for label prior distribution learning and matrix completion. The final loss function can be obtained by combining Equation (3) to Equation (6) as whereŶ = Z (Y d ) denotes filling an observed label matrix with Y d in the corresponding down-sampling location, α and β are the weighting coefficients that indicate the importance of these two terms in first-order priority. First, Y d is obtained with a light classification network. Then, the label matrix completion task enters into solving Equation (8), which is a differentiable squared function with a closed-form solution for Y. The solution for this pixel i is [21,48]: where The final classification label is an index of the maximum value in Equation (8) for pixel i.

Complex Analysis
Considering the proposed optimization scheme of the proposed model in the preceding section, the entire complexity comprises two aspects corresponding to known labels learning with a light network and label matrix completion, respectively. Specifically, as for the known labels learning task, the computation complexity for every layer is O((ĥŵ) 2 E 2 d l d l−1 ), where E is the size of the convolution kernel and d the number of latent units. The computational complexity relies heavily on the size of the down-sampled PolSAR image, which is uniformly sampled from the entire image. The complexity is lower when the down-sampling rate is smaller. Regarding label matrix completion, the computational cost is illustrated in Equation (7) with O((hw) 2 . In theory, it can be observed that the total time computational complexity of the proposed method is acceptable.

Results
In this section, the proposed structure label matrix completion method is evaluated with several experiments. First, the effectiveness of the proposed method on a PolSAR image is investigated with different down-sampling rates. Next, the effectiveness of the structure priority distribution for the entire label matrix is validated. Finally, the effectiveness of the proposed label matrix completion with a structure prior for two different PolSAR datasets is evaluated to establish the superiority of our proposed framework. All experiments are conducted on a desktop PC equipped with an i5 CPU@3.2 GHz and 16 GB of memory in the TensorFlow environment. In the following, the total accuracy is the number of well-classified pixels divided by the number of all pixels. Meanwhile, the classification accuracy and kappa coefficient are given in Tables 3 and 5 on two different PolSAR datasets, respectively. The overall accuracy and kappa coefficient are calculated by the following equation [21,26]: where Nc denotes the number of pixels with the correct label compared with the "Total" number of the pixels N. n kk is the diagonal element of the confusion matrices and n k,: , n :,k are the k-th row and k-th column of the confusion matrices, respectively. They denote the number of classes k divided into the k-th class and k-th class separated into others, respectively.

Parameter Analysis
Information for the PolSAR database tested and detailed here is shown in Figure 4. Figure   Parameters play an important role in the classification performance, which is usually chosen within a range. The parameters to be tested are the coefficients α and β before the MRF term and Laplace term, the aims of which are to retain smoothness in one region and the edge between two regions, respectively. In our experience, we vary the parameter from 10 −3 to 10 for α and 10 −10 to 10 −1 with the interval of 10 −1 for β. Classification accuracy versus different parameter setting values is shown in Figure 5. On the one hand, the classification-accuracy values are almost unchanged as the change of parameter β, which reveals that the Laplace term is robust in the proposed method. On the other hand, classification accuracy increases as parameter α varies from 10 −3 to 1, while decreasing from 1 to 10 and thus reaching the highest value when α is 1. As a whole, it can be seen that the classification-accuracy values exhibit little difference for various α. Thus, in the following experiments, the chosen parameter values are 1 and 10 −4 for α and β, respectively.

Classification Performance with Different Down-Sampling Rates
To validate the robustness of the label matrix completion method, the down-sampled PolSAR images with different down-sampling rates were tested on dataset Flevoland, as shown in Figure 4. In Table 1, the number of pixels in the down-sampled images are given to illustrate the computational cost compared with the entire image. First, the classification results are shown in Figure 6a-d with down-sampling rates of 1 4 , 1 9 , 1 16 , and 1 25 , respectively. Then, a priority matrix for the entire image is obtained by adding the labels into the corresponding locations in the large PolSAR image, as shown in   In Table 1, it can be seen that many fewer samples are to be classified with the down-sampling operation, which results in low computational and storage costs. In this way, a fast classification process will be realized. The computation times used by different down-sampled images and the entire image are given in Table 2, which shows that less time is needed for down-sampled images than that for the entire image. From the classification accuracies in Table 2, it can be concluded that satisfactory values are obtained when the sampling rate is 1 4 , 1 9 , and 1 16 . The accuracy values of these label matrix completion results are 99%, 97%, and 96%, respectively. Even though the accuracy of the training pixels used is 0.6%, the classification accuracy value is high enough. More importantly, the classification accuracy values are good enough for almost all categories when the down-sampling rate is 1 25 . The classification accuracy is poor because the down-sampled image would significantly destroy the structure information of the entire image. It is relatively obvious that pixels in the minor class Building could not be classified correctly. It also results in low accuracy values for down-sampled images with 1 16 and 1 9 sampling rates. However, the reduction of accuracy value is relatively low compared with the computation time, as shown in Table 1.
The runtimes in Table 1 were tested on the down-sampled images and the entire images with the same light classification network setting as in [50]. The classification results correspond to Figure 6 and Table 2. It can be concluded that the computation times of different down-sampled images are proportional to the sampling rates. The time for the 1 25 down-sampling-rate image is only 13.765 s, which is less than the 344 s the whole image needs. It costs only 1 4 time for the whole image compared to that for the 1 4 down-sampling data. Similar results are obtained with the other down-sampled PolSAR images. Therefore, with a suitable light network, a huge PolSAR image can reach a satisfactory classification result with much lower computation time and storage. This situation makes the proposed method able to be widely applied in all the excellent classification methods of the light network to process the down-sampled PolSAR images. Regarding the visual classification results with the proposed label matrix completion method, experiments were conducted on images with different down-sampling rates. It is obvious that the visual results of the down-sampled images are satisfactory, as shown in Figure 6a-d, which have the correct structure and smooth regions. Almost every terrain can be classified correctly in these images. In the completed label prior map images in Figure 6e-h, the structure of regions are retained for the special down-sampling strategy with an adaptive prior learning method. With different sampling rates, the completed label prior image varies greatly when more information is obtained with a high sampling rate, as in Figure 6e. Thus, the final completed label matrix of the image with high sampling rate is better than that with low sampling rate. Figure 6i illustrates better performance in keeping the structure of regions, while the label matrix in Figure 6l with a 1 25 sampling rate exhibits obvious over-fitting where regions in the unconnected part are continuous, such as the terrain Grass. The results of images with 1 9 and 1 16 sampling rates are satisfactory in keeping region structure and smoothness. In general, the performance of label matrix completion is excellent for PolSAR image terrain classification on accuracy values, visual results, and computation time.
To show the robustness of the proposed method with different labeled data and sampling rates, the box plots of different sampling rates are given in Figure 7. From the plot results, it can be seen that the classification-accuracy value increases with decreasing sampling rate, because the higher classification-accuracy value is obtained with more labeled data. With the same sampling rate, the classification-accuracy value is stable. The variance of the classification-accuracy value decreases with different sampling rates from 1 25 to 1 4 . The classification-accuracy result of the down-sampled image is the lowest with a 1 25 sampling rate. Furthermore, the completion results will vary greatly for fewer known labels in the inversion of the matrix completion problem.  Figure 7. Classification accuracy vs different sampling rates for five runs.

Structure Prior for Label Matrix Completion
The effectiveness of structure prior information for label matrix completion was validated. All of the experiments were carried out on the PolSAR image Flevoland database from NASA/JPL AIRSAR with a 1 16 down-sampling rate. The classification visual results and accuracy values are given in Figure 8. Figure 8a is the result of the down-sampled image, where the accuracy value is 94.94%. Filling the down-sampled image into the entire image as shown in Figure 8b, the overall accuracy value is 5%. In addition, the structure of regions can be retained by the uniform sampling method as in Figure 8b. Thus, the classification result with the MRF prior is good in terms of both the visual result and accuracy value in Figure 8c. However, the boundary between different regions is not smooth in the result without the Laplace priority, and the classification visual result in Figure 8d is the best for every region being homogeneous and the boundary is smooth. This is obvious in the regions that are outlined by red boxes as shown in Figure 8c,d.

Classification Performance on PolSAR Data Sets
It has been proved that the proposed label matrix completion method for PolSAR image classification can obtain satisfactory results with low computation times under different down-sampling rates. The classification performance of visual results and accuracy values was compared with several state-of-the-art methods, the results of which were taken from several cited papers.

Dataset Flevoland from NASA/JPL AIRSAR
The methods compared here with the Flevoland dataset include the former works variational mixture Wishart (VMW) [26] and robust-semi-supervised (RS) [21] methods, a traditional machine-learning classifier SVM [29], CNN-based methods, a real-valued CNN (RV-CNN) [32], a complex-valued CNN (CV-CNN) [51], and a graph semi-supervised (GSS) method [44]. The method called CLSL is a former work-and cost-sensitive latent space learning model for imbalanced PolSAR images [50]. The proposed SLMC is classifying a down-sampling PolSAR image with the CLSL method following label matrix completion. All of the compared methods are for PolSAR semi-supervised classification consisting of pixel-and region-based methods. The classification results for these methods are shown in Figure 9, and the accuracy values are given in Table 3.
The image in Figure 9a is the classification result of VMW, which is a pixel-based method heavily influenced by the speckle noise in the PolSAR data. A similar situation also occurs in the classification result of SVM, as shown in Figure 9b. Figure 9c of the RS method has a better classification result with the feature-learning process for the original data. The results in Figure 9d,e are images classified by the RV-CNN and CV-CNN methods, respectively. Both of these latter methods can obtain smooth regions for spatial feature learning for pixels with its neighboring pixels. By computing a graph among the entire image, graph-based methods in Figure 9e can also obtain excellent classification results. However, graph-based methods are heavily limited by the computation of the graph for the entire image, which increases the computational and storage burden. This situation occurs in the deep-network-based methods, such as the RV-CNN and CV-CNN methods, as shown in Figure 9d,e, which also exhibit good classification results. In Figure 9f, the classification result of the proposed CLSL method is good for an imbalanced PolSAR image. The last image in Figure 9g is the result of the proposed label matrix completion method with a 1 4 down-sampling rate. This method can realize excellent classification results with less time and storage.  [5]. (b) WBRBM classifier [9]. (c) WDBN classifier [9]. (d) WDSN classifier [10]. (e) CNN classifier [11]. (f) MLRsubMLL classifier [58]. (g) Wishart classifier [1]. (h) SSA classifier [27]. (i) SNC classifier [29]. over-smoothed and misclassified after iterations of optimization. For example, the road area [marked by red rectangles in Fig. 5(j) and (k)] which should be long and continuous is cut to several disconnected segments for our proposed methods. This is because these small object areas are mixed with pixels of different terrain classes, or the edges of them are weak and unclear. Therefore, the small objects shrink or even disappear after rounds of successive MRF smoothing. The time complexity of our proposed model depends on the time complexities of the two tasks conducted in the optimization process, i.e., CNN learning and label propagation. The time complexity of CNN learning is determined by the number of convolutional layers, the filter number in each layer, and the input and output feature channel numbers of each layer in the CNN, which can be referred to [59]. The belief propagation algorithm is employed to solve the label propagation problem in this paper, which has a linear time complexity [54]. The total running time of the KWGraphCNN method for the Flevoland area data is 14 min 35 s, including the network training (13 min 38 s) and label propagation (57 s).
Regarding the classification accuracy in Table 3, a similar conclusion can be obtained. First, the graph-based and deep-network-based methods obtain higher overall classification accuracy from 33% to 51% than the pixel-based VMW method. The number of training pixels used in the VMW method is unified as 300 for every class, which has a robust result for the number of training pixels. The proposed method obtains the same accuracy of 99% with half the training pixels used in the GSS method and 1 4 used in the CLSL method. Compared with deep-network-based classification methods, e.g., RV-CNN, CV-CNN, and SRDNN, the proposed method can obtain an accuracy value higher than 5%. Furthermore, it can be seen that 1 4 computation time is needed compared with the other methods for the down-sampling process of the entire image, as shown in Table 4. Semi-supervised methods, e.g., the VMW, RS, RV-CNN, and CLSL methods, need feature learning for all the pixels by a parameterized network. The SRDNN method costs more time than the above methods for an extra constraint of super-pixel information for the network. There is an extra relation matrix to be computed for all the pixels of graph-based semi-supervised methods, e.g., the SAG and GSS methods. Thus, these methods cost more time than the others. The proposed method only processes a down-sampling of pixels, resulting in less computation time and space. The overall time is proportional to the number of pixels to be processed. Thus, it is believed that the computational burden will be much less with a lower sampling rate.

Dataset Oberpfaffenhofen from ESAR Airborne Platform
The Pauli RGB pseudo-colored images database tested is shown in Figure 10. Figure 10a is the Pauli RGB image of an L-band multi-look PolSAR database called Oberpfaffenhofen of size 1300 × 1200 obtained by the ESAR airborne platform provided by the German Aerospace Center. Figure 10b,c are the ground-truth image and all terrains in this dataset, respectively. There are three different terrains in this dataset, including built-up area, woodland, and open area, as shown in Figure 10j. Terrains in this PolSAR dataset are balanced and homogeneous in a large area, which makes it suitable for the proposed label matrix completion method with a large down-sampling rate. With a few pixels labeled, the final classification result is obtained with satisfactory performance. The classification visual results and accuracy values are given in Figure 10c-i and Table 5, respectively. The methods compared on this dataset are the SVM [29], RS [21], and CNN-based RV-CNN [32] and CV-CNN [51] methods, and the SRDNN method [52].     (j) Figure 10. Pauli RGB pseudo-colored images and their terrain information from PolSAR database Oberpfaffenhofen from ESAR Airborne Platform: (a) Pauli RGB pseudo-colored image, (b) ground-truth image, (c-i) classification results on dataset Flevoland with different methods, including SVM [29], RS [21], RV-CNN [32], CV-CNN [51], SRDNN [52] methods, and our former CLSL [50] model and proposed SLMC method, respectively, and (j) terrains in the image in (a).
From the classification visual and quantified performance, it can be seen that a different method could achieve a similar label map. The largest accuracy value of the proposed method is 15% higher than the smallest one. While only a 1% accuracy value is higher than that of most of compared methods, it validates that all the methods have excellent performance on this dataset. However, the training pixels used in different methods vary greatly. The proposed method uses a 1 64 down-sampling image with only 0.15% training pixels and can realize an excellent classification result. CNN-based methods realize this with 1% pixels, which is 10 times the number used in the proposed method. Not only are the computation time and storage reduced, but the human cost in label tagging is also relieved. The classification visual result in Figure 10d is produced by the RS method, which aims to make the method robust to noisy PolSAR data and labels. This method falsely classifies several built-up areas as woodland. In addition, there are some significantly noisy blocks in the visual results of the RV-CNN and CV-CNN methods. The visual result in Figure 10g of the SRDNN method has some misclassified pixels in the top right-hand corner of the terrain buildup. The visual result of the proposed method is relatively satisfactory in the entire image.

Conclusions and Future Work
In this paper, the PolSAR image classification task is first tackled from a label matrix completion perspective. The known labels are learned by a light network that processes a uniform down-sampled PolSAR image to relieve the heavy computational burden and retain the regional information of the entire PolSAR image. Then, the label transfer framework realizes the final label matrix completion for the entire PolSAR image. Compared with other state-of-the-art classification methods, the proposed label matrix completion model can achieve good classification results. It can be concluded from the experiments that the proposed method realizes excellent classification results at a low computation-time cost. In addition, from the experiments described in Section 4, it can also be found that PolSAR images with continuous large areas are more suitable for the proposed method than those with small and imbalanced regions.
In fact, the proposed structure label matrix completion method for PolSAR classification task is a label prior transfer process. This task could be realized to transfer information from a down-sampling PolSAR image as in this paper and from the different data sources such as infrared or optics. Thus, the future work will pay attention to the transfer learning for label prior of the label matrix completion task in PolSAR classification.