Lithology Discrimination Using Sentinel-1 Dual-Pol Data and SRTM Data

Compared to various optical remote sensing data, studies on the performance of dual-pol Synthetic aperture radar (SAR) on lithology discrimination are scarce. This study aimed at using Sentinel-1 data to distinguish dolomite, andesite, limestone, sandstone, and granite rock types. The backscatter coefficients VV and VH, the ratio VV-VH; the decomposition parameters Entropy, Anisotropy, and Alpha were firstly derived and the Kruskal-Wallis rank sum test was then applied to these polarimetric derived matrices to assess the significance of statistical differences among different rocks. Further, the corresponding gray-level co-occurrence matrices (GLCM) features were calculated. To reduce the redundancy and data dimension, the principal component analysis (PCA) was carried out on the GLCM features. Due to the limited rock samples, before the lithology discrimination, the input variables were selected. Several classifiers were then used for lithology discrimination. The discrimination models were evaluated by overall accuracy, confusion matrices, and the area under the curve-receiver operating characteristics (AUC-ROC). Results show that (1) the statistical differences of the polarimetric derived matrices (backscatter coefficients, ratio, and decomposition parameters) among different rocks was insignificant; (2) texture information derived from Sentinel-1 had great potential for lithology discrimination; (3) partial least square discrimination analysis (PLSDA) had the highest overall accuracy (0.444) among the classification models; (4) though the overall accuracy is unsatisfactory, according to the AUC-ROC and confusion matrices, the predictive ability of PLSDA model for limestone is high with an AUC value of 0.8017, followed by dolomite with an AUC value of 0.7204. From the results, we suggest that the dual-pol Sentinel-1 data are able to correctly distinguish specific rocks and has the potential to capture the variation of different rocks.