.Surface cover diversity and the complexity of geological structures can seriously impact the accuracy of mineral mapping. To address this issue, we propose a method for lithological classification and analysis based on random forest (RF) and multiple features. Feature vectors, including spectral, polarization, texture, and terrain features, are constructed to provide multidimensional information. Subsequently, these feature vectors are screened based on their discriminative properties for different lithologies to reduce feature redundancy. Finally, the results of lithological classification can be obtained using the RF algorithm based on the selected features. In the experiments conducted in the Qulong copper deposit area, data from Sentinel-1A, Sentinel-2A, and Terra satellites were used to extract multidimensional features. After calculating the Bhattacharyya distance and analyzing the probability density distribution, 17 features selected were input into the RF classifier, achieving an accuracy of 88.83% in lithological classification. This represents a 7.5% improvement compared to exclusively relying on spectral features, and suggests that the proposed method of combining spectral, polarization, texture, and terrain features provides new possibilities for improving the accuracy of field lithological classification.