Deep Learning Based Lithology Classification Using Dual-Frequency Pol-SAR Data

Lithology classification is a crucial step in the prospecting process, and polarimetric synthetic aperture radar (Pol-SAR) imagery has been extensively used for it. However, despite significant improvements in both information content of Pol-SAR imagery and advanced classification approaches, lithology classification using Pol-SAR data may not provide satisfactory classification accuracy due to high similarity of certain classes. In this paper, a novel Pol-SAR lithology classification method based on a stacked sparse autoencoder (SSAE) is proposed. By using superpixel segmentation, new features can be extracted from dual-frequency Pol-SAR data, which can increase the class separability of the input data. Then, these features and the coherency matrices are incorporated into SSAE to classify the lithology. The classification performance is evaluated on an SIR-C dataset acquired over Xinjiang, China. The experimental result shows that this method is effective for lithology classification and can improve the overall accuracy up to 98.90%.


Introduction
Synthetic aperture radar (SAR) imagery has the capabilities of providing high resolution images for a large area under all-weather and all-day conditions.Therefore, it has broad applications in areas such as geological exploration, environmental monitoring, and natural disaster assessment [1].Moreover, compared with traditional single polarized SAR images, full polarized SAR images (Pol-SAR), which contain four different Pol-SAR channels (e.g., hh, hv, vh, and vv), can provide more abundant information on land cover characteristics.It is able to provide not only the radar scattering cross section (RCS) but also the phase information of the four channels, which can reveal the scattering mechanism of targets from various aspects, such as structure and surface undulations.Hence, the classification of different lithologies based on Pol-SAR images is becoming an important tool.
In the last few decades, a variety of land cover classification algorithms have been proposed [2].They can be broadly categorized into supervised approaches, semi-supervised approaches, and unsupervised approaches according to whether manual annotations are utilized.Kong proposed a maximum-likelihood based classification method for Pol-SAR images based on a complex Gaussian distribution [3].This classifier was applied to single look Pol-SAR images and was modified to multi-look situations by Lee [4].In this work, many features are extracted from Pol-SAR images, which can describe the targets from various aspects.The application of features makes it an important field in which to conduct land cover classification based on polarimetric features.For instance, Cloude and Pottier defined three features including entropy (H), anisotropy (A), and alpha (α) by eigen decomposition.The H and α are widely used for classification [5].Freeman-Durden [6] determined three-component scattering mechanisms based on physical models.In this method, double bounce, volume, and surface scattering correspond with the double-bounce scattering from a dihedral corner reflector, randomly oriented thin cylindrical dipoles, and first-order Bragg scattering, respectively.The Freeman-Durden decomposition is useful for decomposing the scattering mechanism from naturally incoherent scatterers.In Reference [7], Qi synthesized polarization decomposition, and proposed an object-based Pol-SAR image classification based on a decision tree.Through mathematical and physical analysis of Pol-SAR data, Mahdianpair presented a modified coherency matrix.The classification accuracy of wetlands was improved to 92.17% by combining the modified coherency matrix and a random forest algorithm [8].
Recently, deep learning have been proven to be very effective in classification and recognition as it can extract features automatically [9,10].Therefore, a Pol-SAR image classification based on deep learning has attracted the attention of experts [11].Hansch proposed the usage of complex valued neural networks for classification of different land in Pol-SAR images, which confirmed the possibility of deep learning applied to Pol-SAR image classification [12].Lv et al. used a deep belief network as one of the deep learning models to classify different urban land covers in a Pol-SAR image [13].In Reference [14], Zhou introduced a convolutional neural network into Pol-SAR image classification.Experimental results with Pol-SAR data indicated that the deep learning methods could provide competitive results.Xie proved the suitability of a stacked sparse autoencoder (SSAE) in the classification of Pol-SAR images [15].However, this method is pixel-based, without considering the local spatial information.To solve this problem, Zhang et al. took all the neighborhood pixels around the cell under test as the input data, and then the Pol-SAR images were classified by SSAE [16].This method integrated the contextual information of the neighborhood, but it has defects in practical applications especially when the target area is complex.In References [17,18], Hou and Geng proposed improved land cover classification methods based on superpixels to use the local spatial information more properly.In this study, SSAE is selected due to several advantages compared to other methods.First, unlike the convolutional neural network, SSAE is pixel-based, thus, it can classify any pixel including boundary pixels.Second, SSAE can easily handle large and multi-dimension remote sensing datasets, which is especially suitable for this study.Finally, it has a flexible and straightforward structure and has shown good results in various remote sensing applications.
In the processing of classifying different vegetation covers or different urban land covers using Pol-SAR, effective features extraction is the foundation to obtain good results.However, it is difficult to get satisfactory classification features just by using one band Pol-SAR data for lithology, since lithologies are affected by surface weathering over long periods of time, and the surface characters thus become similar among some classes [1].To classify different lithologies, dual-frequency Pol-SAR will play a role because of the different wavelengths and varying penetrability to vegetation cover.L-frequency has the longer wavelength and better penetrability.SSAE can extract some features according to its own rules, but it cannot extract all possible features, especially those with specific physical meanings.Combining the effective handcrafted features with raw data may improve the lithology classification performance of SSAE.In this paper, new polarimetric features were proposed based on dual-frequency Pol-SAR data for the surface fluctuations of different lithologies, and then SSAE was used to extract deep features and classify different lithologies.To the best of our knowledge, this is the first attempt at the classification of lithology by combining SSAE and dual-frequency Pol-SAR data.
The paper is organized as follows: in Section 2, we explain the polarimetric features extraction and SSAE.This section also presents the new lithology classification method and describes the dataset used in the experiments.The proposed method is verified and analyzed in Section 3, at the same time, the proposed method is also compared with other methods in this section.A discussion is found in Section 4. And some useful conclusions are drawn in Section 5.

Pol-SAR Features Extraction
Polarimetric radar measures the complex scattering matrix of a medium with quad polarizations.The scattering matrix in the linear orthogonal polarization basis can be expressed as [19]: where s hv is the scattering element of horizontally polarized transmitting and vertically polarized receiving waves, and the other three elements are defined similarly.In the case of reciprocal backscattering, s vh = s hv .The corresponding vectors k T and k C can be expressed by the following forms [19]: where [ ] T is an ordinary transpose operation.
For the distributed targets, the 3 × 3 coherency matrix T and covariance matrix C can be computed by [19]: where * stands for conjugate, and H denotes conjugate transpose.Both the coherency matrix and covariance matrix contain the scattering characteristics of a distributed target.Pol-SAR data can provide characteristics of the land cover, such as permittivity, surface roughness, material, and so on.To clearly show the difference between different lithologies, some new polarimetric features should be defined and extracted based on the cross polarized ratio, co-polarized correlation coefficient, and the proportion of different scattering mechanisms, etc.

Cross Polarized Ratio
The cross polarized ratio is the echo ratio between the cross channel and co-polarized channels.It can be defined as follows with respect to the coherency matrix [20]: where r xv and r xh are the cross polarized ratios in vv channel and hh channel, respectively.The cross polarized ratio reflects the effect of surface roughness on scattering, which is proportional to the surface roughness, within a certain roughness range.

Co-Polarized Correlation Coefficient
The hh and vv channels may have different echoes due to different scattering mechanisms.The co-polarized correlation coefficient can be used to describe the scattering characters of the illuminating area, and it is defined as Equation ( 8) with respect to the coherency matrix [20]: The co-polarized correlation coefficient is sensitive to permittivity and surface roughness [21][22][23], which is proportional to water content when the roughness is in a certain range.

Freeman-Durden Decomposition
The Freeman-Durden decomposition determines three-component scattering mechanisms, including double-bounce, volume, and surface scattering, based on physical models.The Freeman decomposition is shown as follows [24]: where [C] v , [C] d , and [C] s are the covariance matrices of the three scattering mechanisms, and f v , f d , and f s are the scattering coefficients of double-bounce, volume, and surface, respectively.α and β are the polarized coefficients of double-bounce and surface, respectively.P v , P d , and P s are the scattering powers of the three scattering mechanisms.They can be used as features describing the scattering mechanisms because they are related to echo power from the different mechanisms.

Stacked Sparse Autoencoder
The stacked sparse autoencoder is a deep learning architecture that can extract the deep features of the input data and achieve land classification by the softmax classifier.It has been used widely in image processing and pattern recognition due to its simple process and high classification accuracy.

Sparse Autoencoder
A sparse autoencoder is a feedforward, fully connected, nonrecurrent neural network with an input layer, a hidden layer, and an output layer.The output layer has the same number of nodes as the input layer, and the hidden layer has more nodes than the input layer [24,25].The sparse autoencoder self-learns so that the output get close to the input with a teacher vector identical to the input vector.This means that, the sparse autoencoder is able to extract an optimized feature set from the input data.The action of the sparse autoencoder is described as follows: Step 1: Initialize the weight matrix and the bias vector, and calculate the outputs of the hidden layer and output layer: x = g W (2) y + b (2) , (14) where y and x are the outputs of the hidden layer and output layer, respectively.f (x) and g (x) are the activation functions, and f (x) = 1 1+exp(−x) , g (x) = x.W (1) represents the weight matrix between the input layer and the hidden layer.W (2) is the weight matrix between the hidden layer and the output layer.b (1) and b (2) are the bias vectors of the hidden and output layers, respectively.
Step 2: Calculate the loss function based on outputs.The loss function used in this paper has the following form [18]: where m is the number of training data points.s l (l = 1, 2, 3) represents the number of neurons in the input layer, the hidden layer, and the output layer, respectively.The first term on the right side of Equation ( 15) is the traditional squared error, which stands for the difference between the output and the input.The second one is the regularization penalty which can penalize large values of the parameters to get a simpler network.This can prevent network overfitting.The third part is the sparse penalty term to improve the sparsity of the encoding representation and the Kullback-Leibler divergence KL (ρ ρi ) is defined in Equation ( 16): where ρi is the average activation of the i th hidden unit over the training set.ρ is the sparse parameter which is a small value close to zero.The constraint ρj = ρ is enforced during the training.If the constraint is satisfied, the hidden activities are mostly close to zero, which ensures that the hidden layer can be forced to learn a compressed representation.
Step 3: Minimize the loss function to optimize W (l) and b (l) .It can be seen from step 2 that to get a satisfying network, we need to minimize the loss function by optimizing W (l) and b (l) (l = 1, 2).A stochastic gradient descent algorithm [26] is applied to minimize the loss function, and the process is as follow: where γ is the learning rate.

Stacked Sparse Autoencoder
A hierarchical training strategy is used in SSAE to construct the deep neural network to generate deep features of the input data.First, we establish several SAEs, one by one.Then train the first SAE with input features.After the first SAE is trained, regard the features obtained from the hidden layer as the input of the next SAE, for training.In this way, each SAE can be trained separately, which is called pre-training [27].After pre-training, the decoder layers of all SAEs are removed and a softmax classifier is connected to accomplish the classification.Hence, the trained SAEs and the softmax classifier are connected to establish a whole network, which is called the SSAE.Finally, in order to obtain a more accurate classification, back-propagation is applied to adjust the parameters of the SSAE with the training samples and their labels.This step is called fine-tuning, which treats all layers of the SSAE as a single model.
After pre-training and fine-training, a trained SSAE can be obtained.We can extract optimized features from the input data to obtain favorable results in classification.

Lithology Classification Based on Deep Learning
It is challenging to classify lithology because the surface has suffered from weathering for years.In this paper, a novel classification method is proposed for lithology classification based on dual-frequency Pol-SAR data and the SSAE classifier.The proposed method has the flow chart as shown in Figure 1, which includes pre-processing of Pol-SAR data, polarimetric feature extraction, and classification with the SSAE.SAR images, as a consequence of their coherent nature, are subject to speckle, which decreases image quality and makes image description and interpretation difficult.Thus, a refined Lee filter [28] whose window size is 9 × 9 is applied to the Pol-SAR dataset as a pre-processing process to reduce speckles.

• Superpixel Segmentation
A superpixel is a perceptually consistent unit with the pixels in one group being similar in character and texture.Superpixel segmentation can divide an image into a number of small contiguous regions where the pixels in each region are homogeneous.It is thus a type of oversegmentation of an image.The pixels in the same superpixel are considered as the same category in this paper.
The Pauli RGB image is formed with the intensities of three channels.After oversegmentation, the Pauli RGB image can be segmented as {I 1 , • • • , I N } using the SLIC superpixel generating algorithm [29].After segmentation, all the pixels in a superpixel will be labeled in the same class, which can be expressed by the averaged coherency matrix Tj : where n j is the number of pixels in superpixel I j .Due to the set average processing, superpixel segmentation mitigates the influence of speckle noise compared with the original coherent matrices.To utilize the spatial features of Pol-SAR data, the polarimetric features used in the novel classification are all extracted from the superpixel segmented image.

Lithology Feature Extraction
Despite using advanced remote sensing tools, classification of lithology is challenging due to high similarity caused by weathering and coverage, which can produce confusion in the classification.Thus, in this study, we extract new features to increase the separability among different classes and to improve the classifier performance.
The volume and surface scattering power obtained from Freeman decomposition is used in this study to indicate the differences in structures.The cross polarized ratio measures the roughness of the surface, and the co-polarized correlation coefficient is related to the permittivity.The lithology can be well characterized by combining these features.The lithological differences are sensitive to the wavelength of electromagnetic waves.Therefore, we can improve the classifier by synthesizing different frequency Pol-SAR data.New features are defined in this paper based on the cross polarized ratio, co-polarized correlation coefficient, and the Freeman-Durden decomposition to indicate how features change with different wavelengths: ∆P s = P s_C − P s_L where r xv_i , ρ _i , and P v_i and P s_i (i = L, C) are the cross polarized ratio, co-polarized correlation coefficient, and volume and surface scattering powers under L and C bands.Equations ( 20)-( 23) represent the trend and variation of the target characteristics when the imaging frequency changes.This indicates the differences in the lithology scattering mechanisms.Scattering mechanisms corresponding to target feature variations are shown in Table 1.Compared with the single-frequency Pol-SAR image classification, the proposed method based on dual-frequency Pol-SAR images exploits the trend and variation of the features with frequency change rather than static and isolated scattering mechanisms.Therefore, the new classification feature vector is constructed, including T 11_i , real(T 12_i ), imag(T 12_i ), real(T 13_i ), imag(T 12_i ), T 22_i , real(T 23_i ), imag(T 23_i ), T 33_i , (i = L, C), ∆r xv , ∆ρ, ∆P v , and ∆P s , which is utilized to distinguish different lithology.

Lithology Classification
In this paper, SSAE is used to classify the lithology.Every pixel in the Pol-SAR image can be used for feature vector extraction.Then, we determine an optimal configuration of SSAE according to the classification results.The determined SSAE can be trained by the training dataset.Finally, lithologies can be classified by the trained SSAE.
The optimal architecture of SSAE can be obtained if the number of hidden layers and neurons of each layer are fixed.The steps are as follows.
Step 1: Determine the optimal number of layers.The dataset is used for training on different SSAE topologies with increasing depth but the same neurons.The classification accuracy is analyzed to compare different architectures.In this way, we can obtain the SSAE with most suitable number of layers for this dataset.
Step 2: Determine the suitable number of neurons in each hidden layer.tTe dataset is used for training on different SSAEs.The number of hidden layers in each SSAE has been determined by step 1, and here the number of neurons is varying.Then we can obtain the optimal number of neurons by evaluating the classification accuracy of these SSAEs.

Dataset Description
The experimental data used in this paper are acquired by the NASA/JPL SIR-C system in Qinghe Xinjiang, China, in 1994, including the 4-Look Complex L-Band and C-Band datasets with a resolution of 12.5 m × 12.

Results
In this section, the performance of the proposed method is verified and analyzed.

Experimental Parameters
Our algorithm comes with some parameters that can affect the classification performance, such as the number of hidden layers and nodes in the SSAE, and the number of superpixels and training pixels.Specifically, the number of hidden layers and nodes in each layer can influence the performance of the SSAE.Moreover, superpixels divided from the Pauli RGB image can affect the classification results.A smaller number of superpixels implies that a superpixel contains more pixels and therefore more contextual information of the neighborhood can be considered.However, the classification of boundaries becomes poor when the number of superpixels is too small.Furthermore, the number of training pixels also influences the estimation accuracy.In order to analyze the impact of these parameters, experiments are carried out with the Qinghe dataset.
To determine the optimal number of layers, SSAE topologies with increasing depth were trained with the same number of superpixels and training pixels.In this experiment, the SIR-C image was divided into 2200 superpixels, and 10% of labeled pixels were randomly selected as the training set to train the network.The overall accuracy (OA), which is defined in Equation (24), is calculated to compare different architectures.From the plot shown in Figure 5, it can be seen that a three-layer SSAE is the most suitable for the dataset.Furthermore, the network has 60, 80 ,and 100 nodes in these hidden layers.From the above mentioned results, a three-layer SSAE with 60, 80, and 100 nodes was applied in the following experiments.The OA is: where N is the number of total samples, r denotes the number of classes, and x ii is the diagonal elements in the confusion matrix.Then, the optimal number of superpixels obtained from the SIR-C image can be determined by following experiments.In this step, the SIR-C image is divided into 1000, 1200, 1400, 1600, 1800, 2000, 2200, 2400, 2600, 2800, 3000, and 3200 superpixels, and then classification performance is analyzed to find the optimal number of superpixels.As shown in Figure 6, 1800 superpixels is the favorable parameter.
Finally, we validate the performance of our network architecture with different numbers of training pixels.We respectively choose 0.1%, 0.2%, 0.5%, 1%, 2% , 3%, 5%, 10%, 15%, and 20% of labeled pixels as the training set to train the network.From Figure 7, we find that, when the rate of training samples increases from 0.1% to 5%, the OA of classification increases sharply at first, and then a stable and high OA is obtained.5% of the labeled pixels is enough to train the network, and it can prevent overfitting.Thus, as shown in Table 2, 5% of the labeled pixels are randomly selected as a training dataset, and the remainder are regarded as the testing dataset.

Classification Results
In this experiment, using the proposed method, the Pol-SAR images can be classified.First, the data should be pre-processed.A refined Lee filter with 9 × 9 window size is chosen to reduce the noise.Then, after superpixel segmentation, 1800 superpixels remain, and the result is shown in Figure 8.It can be seen that most pixels belonging to the same superpixel have nearly the same color.This means that these pixels have nearly the same scattering mechanism, and can be classified as one group.
Then, polarimetric features can be extracted from the Pol-SAR data, and the features maps are shown in Figure 9.It can be seen that when just using the original features, some classes can not be distinguished, however significant differences can be found in the features maps between these two classes when the new features are used.
Finally, extracted features are entered into the SSAE which has three hidden layers with 60, 80, and 100 nodes, respectively.The classified map is depicted in Figure 10.To varify the advantages of the new features compared to the original features (P s_i , P v_i , ρ hh_vv_i , and r xv_i (i = L, C)), comparison experiments and analysis are conducted.In this experiment, use d i,j , which is defined as follows, to describe the difference between class i and j: where fi and fj represent the average of class i and j with feature f , respectively.d i,j is a number greater than 1, and it is proportional to the difference between two classes.From Figures 3 and 4, we can directly see that Beitashan Formation (class 2) and Upper Aermantie Formation (class 5) are difficult to distinguish, Diluvium (class 3) is easily confused with Alluvium-Diluvium (class 6).Therefore, d 2,5 and d 3,6 are calculated in the regions marked in Figure 9a with original features and new features.
As can be seen in Table 3, d 2,5 and d 3,6 for the new features are greater than those for the original features, which indicates that the new features are more effective in increasing the class separability than the original features.For example, for the situation that just uses r xv_L , d 3,6 is 1.06.However, d 3,6 is improved to 1.31 when ∆r xv is used.As shown in Figure 10, the proposed method can generate competitive visual effects with good connectivity and smooth borders, compared with the ground truth.Thus, lithology can be classified in a good way with this method.
Although favorable results of lithology classification can be obtained with the proposed method, several different classes of lithology are confused.The confusion matrix is shown in Table 4.The Diluvium is mostly confused with Alluvium-Diluvium, leading to a low classification accuracy of these classes.This is primarily due to the fact that both Diluvium and Alluvium-Diluvium have a low back-scattering power.

Comparison with Other Classifiers
To validate and test the performance of the proposed method, the SIR-C data of Qinghe Xinjiang, China, was used.The experimental verification included two aspects, which are the classification accuracy and computational burden.

Classification Accuracy
To validate the lithology classification accuracy, a comparison between the proposed method and other methods is made in this section.The classifiers compared are three SSAE methods with different inputs, the Support Vector Machine (SVM) classifier, and Hou's method from Reference [18].The three SSAE methods in this paper can be labeled as M1, M2, and M3, whose inputs are the coherency matrix of the L-frequency, the C-frequency, and the dual-frequency Pol-SAR data, respectively.Both the SVM classifier and Hou's method utilized the same features as the proposed method.The superpixel segmentation is introduced for these comparison methods.The SVM is a typical method in classification and Hou's method is a state of the art method in Pol-SAR image classification.In the experiment of SVM, 'rbf' is used as the kernel and c = 100, gamma = 0.01.The classification accuracy is reported in Table 5 and the results are plotted in Figure 11.The OA and kappa coefficient (Kappa) are calculated to evaluate the classification performance.Kappa is defined as: where N is the number of total samples, r denotes the number of classes, x ii is the diagonal elements in confusion matrix, and x i+ and x +i are the number of samples in the ith row and the ith column of confusion matrix, respectively.It can be directly seen from Figure 11 that the classification performance is improved significantly compared with the other five methods.As shown in Table 5, the proposed method has a lithology classification accuracy of 98.90% in this experiment, and it outperforms M1, M2, M3, SVM, and Hou's method, which have accuracies of 72.99%, 80.55%, 94.70%, 86.69%, and 80.05%, respectively.The classification accuracy of lithology, especially for classes with high similarity, has been improved with the proposed method.From the results of M1, M2, and M3, we can conclude that dual-frequency Pol-SAR data can provide more information than single frequency Pol-SAR data; the accuracy of M3 is better than M1 and M2.As can be seen from the classification accuracy of the proposed method and M3, the new features extracted from dual-frequency Pol-SAR data are effective and can improve the classification accuracy of lithology.SSAE is more suitable for lithology classification than SVM, because SSAE can extract deep features of input data and classify the lithology with these features.Therefore, the classification result of the proposed method is better than SVM.By comparing the performance of the proposed method and Hou's method, we can tell that the algorithm in this paper is more effective for lithology classification.Generally, the proposed method outperforms the reference methods because new features extracted from dual-frequency Pol-SAR data can increase the class separability of the input data, and SSAE is suitable for lithology classification.

Computational Burden
The computational burden of the proposed method is also an important issue, and the execution times of all methods are given in Table 6.Please note that the time shown in Table 6 is for the training step, which is the most time consuming of all the steps.The computational burden of the SVM is low when the number and dimension of samples are small.However, in this experiment, the number of samples is too large.Therefore, the computational burden of SVM is heavier than other methods.For the methods using SSAE, the time consumed is almost related to the input dimension because they have the same number of samples and the same architecture of the SSAE.From Table 6, we can see that the proposed method has an acceptable computational burden, and the times consumed for the methods using SSAE are significantly less than the method with SVM.

Discussion
The objective of this study is to classify lithology using dual-frequency Pol-SAR data.In this paper, new features are extracted from dual-frequency Pol-SAR data to increase the class separability of the input data, and the SSAE is selected to obtain deep features of the input data for lithology classification.From Figure 10, we can see that lithology is classified successfully with the proposed method.To verify the performance of the proposed method, several methods including M1, M2, M3, SVM, and Hou's method were applied to the same dataset for classification of lithology.The accuracy of M1 and M2 is 72.99% and 80.55%.The accuracy of these methods is low because they only use single frequency Pol-SAR data.The performance of Hou's method is 80.80%.The classification performance of SVM is 86.69% and the SVM consumes more time than the other methods.Both M3 and the proposed method achieved more than 90% classification accuracy.Accuracy and an acceptable time consumed is essential in the classification with the proposed method.Therefore, the proposed method has the best result of lithology classification.
Classification of lithology is important in geological exploration, thus the method proposed in this paper can be used widely in the future.If we can obtain a few labels spread over a specific region by manual investigation, we can propagate this knowledge over the entire region.Besides, we also can train the network on one set of a region and apply it to another region, if the data of these regions are acquired under the same radar imaging conditions.
The current paper used SSAE to extract deep features of input data and classify lithology.In the future, new classifiers can be applied in the proposed method to improve the classification accuracy.Recently, a combination of results from different classifiers has been proposed [30] and our future study will investigate this approach, which may improve the results obtained in this paper.

Conclusions
The identification of lithology from Pol-SAR data is challenging due to weathering and surface occlusion.In this paper, we have proposed a novel classification method based on dual-frequency Pol-SAR data and a deep learning algorithm.This method includes pre-processing, dual-frequency polarimetric features extraction, and classification with an SSAE.In this method, superpixels are produced to integrate contextual information of the neighborhood.Then, new features are proposed to increase the class separability of lithology by combining dual-frequency Pol-SAR data.These features and the coherency matrices of dual-frequency Pol-SAR data are incorporated into an SSAE for lithology classification.The performance of the classifier is demonstrated with measured SIR-C data.The overall accuracy is about 98.90% for the dataset in the experiment.New features are proved to be necessary by comparing with M1, M2, M3, and Hou's method.Meanwhile, it has been demonstrated that the SSAE outperforms SVM for lithology classification.Experimental results confirm the validation and applicability of the proposed method for lithology classification.
5 m.According to the ground truth acquired by the Chinese Academy of Sciences in 1994, there are eight classes of lithology in the experimental area, including Alluvium (Q a1 4 ), Beitashan Formation (D 2b ), Biotite granite, Leucogranite (γ 2c 4 ), Upper Aermantie Formation (D b 3a ), Nanmingshui Formation (C 1n ), and Mayinebo Formation (D 1m ).The reference map of this area is provided in Figure 2. The L-Band and C-Band Pauli RGB images are shown in Figure 3.

Figure 4
is the ground truth where black color means that the pixels have no class labels.

Figure 5 .Figure 6 .Figure 7 .
Figure 5. Changing of the overall accuracy (OA) with the number of layers.

Figure 9 .
Figure 9. Features maps.(a) Scattering power of the surface in the C-band.(b) Scattering power of the surface in the L-band.(c) The changing trend of the scattering power of the surface.(d) Scattering power of the volume in the C-band.(e) Scattering power of the volume in the L-band.(f) The changing trend of the scattering power of the volume.(g) Co-polarized correlation coefficient of the C-band.(h) Co-polarized correlation coefficient of the L-band.(i) The changing trend of co-polarized correlation coefficient.(j) Cross polarized ratio of the C-band.(k) Cross polarized ratio of the L-band.(l) The changing trend of the cross polarized ratio.

Table 1 .
Scattering mechanisms corresponding to target feature variations.

Table 2 .
Training and testing pixel numbers for the Qinghe dataset.

Table 3 .
The difference between two classes.

s_L P s_C ∆P s P v_L P v_C ∆P v ρ _C ρ _L ∆ρ r xv_C r xv_L ∆r xv
Classified map based on the proposed method.