Comparison of Landsat OLI, ASTER, and Sentinel 2A data in lithological mapping : A Case study of Rich area (Central High Atlas, Morocco)

Abstract

The eastern part of the Rich area consists of the massive Paleozoic and Meso-Cenozoic cover formations that present the geodynamic development of the study area, where is characterized by various carbonate facies of Jurassic age. The geographical characteristic of the study area leaves the zone difficult to map by conventional methods. The objective of this work focuses on the mapping of the constituent lithological units of the study area using multispectral data of Landsat OLI, ASTER, and Sentinel 2A MSI. The processing of these data is based on a precise methodology that distinguishs and highlights the limits of the different lithological units that have an approximate similarity of spectral signature. Three techniques were used to enhance the image including Principal Component Analysis (PCA), Minimum Noise Fraction (MNF), and Independent Component Analysis (ICA). Lithological mapping was performed using two types of supervised classification : Maximum likelihood classifier (MLC) and Support Vector Machine (SVM).The results of processing data show the effectiveness of Sentinel 2A data in mapping of lithological units than the ASTER and Landsat OLI data. The classification evaluation of two methods of the Sentinel 2A MSI image showed that the SVM method give a better classification with an overall accuracy of 93,93% and a Kappa coefficient of 0.93, while the MLC method present an overall accuracy of 82,86% and a Kappa coefficient of 0.80. The results of mapping obtained show a good correlation with the geological map of the study area as well as the efficiency of remote sensing in identification of different lithological units in the Central High Atlas.