Next Article in Journal
Evaluation of Temperature and Humidity Profiles Retrieved from Fengyun-4B and Implications for Typhoon Assimilation and Forecasting
Next Article in Special Issue
Remote Detection of Geothermal Alteration Using Airborne Light Detection and Ranging Return Intensity
Previous Article in Journal
Impact of Satellite-Derived Land Cover Resolution Using Machine Learning and Hydrological Simulations
Previous Article in Special Issue
GF-2 Data for Lithological Classification Using Texture Features and PCA/ICA Methods in Jixi, Heilongjiang, China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improving Rock Classification with 1D Discrete Wavelet Transform Based on Laboratory Reflectance Spectra and Gaofen-5 Hyperspectral Data

College of Geo-Exploration Science and Technology, Jilin University, Changchun 130026, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(22), 5334; https://doi.org/10.3390/rs15225334
Submission received: 14 September 2023 / Revised: 7 November 2023 / Accepted: 9 November 2023 / Published: 13 November 2023

Abstract

:
The high intra-class variability of rock spectra is an important factor affecting classification accuracy. The discrete wavelet transform (DWT) can capture abrupt changes in the signal and obtain subtle differences between the spectra of different rocks. Taking laboratory spectra and hyperspectral data as examples, high-frequency features after DWT were used to improve the discrimination accuracy of rocks. Various decomposition levels, mother wavelet functions, and reconstruction methods were used to compare the accuracy. The intra-class variability was measured using the intra-class Spectral Angle Mapper (SAM). Our results show that the high-frequency features could improve the discrimination accuracy of laboratory spectra by 13.4% (from 46.5% to 59.9%), compared to the original spectral features. The accuracy of image spectra in two study areas increased by 8.6% (from 68.3% to 76.9%) and 7.2% (from 81.3% to 88.5%), respectively. Haar wavelets highlighted the spectral differences between different rocks. After DWT, intra-class SAM reduced and intra-class variability of rocks decreased. The Pearson correlation coefficient indicated a negative correlation between intra-class variability and overall accuracy. It suggested that improving classification accuracy by reducing intra-class variability was feasible. Though the result of lithological mapping still leaves room for improvement, this study provides a new approach to reduce intra-class variability, whether using laboratory spectra or hyperspectral data.

1. Introduction

Remote sensing is an important tool for obtaining lithology maps while saving costs [1,2]. Although some geologists argue that remote sensing may not provide accurate results for creating lithology maps, the combination of field surveys and remote sensing can significantly improve the efficiency of lithological mapping [3].
Airborne imaging spectrometers were used for geological mapping in earlier studies [4]. Greenbaum et al. [5] used 11-band airborne multispectral scanner data to discriminate lithology in central Snowdonia. They found that the Optimal Index Factor (OIF) and principal component images were helpful for sediment identification. With the development of satellite remote sensing, large-scale and high-frequency observations became possible [6]. The Landsat series, Sentinel-2, and other multispectral sensors are widely used for lithology classification [7,8]. Pal et al. [9] proposed an image fusion method and applied the method to Sentinel-2, Landsat 8, and ASTER data to improve rock discrimination accuracy. Rowan et al. [10] performed ASTER ratio imaging, matched-filtering, and spectral-angle mapping for lithological classification. Spectral reflectance and emissivity data were utilized to distinguish between four mafic-ultramafic categories, three categories of alluvial-colluvial deposits, and a quartzite unit. In addition to multispectral data, hyperspectral data like Earth Observing-1 (EO-1) Hyperion and GaoFen-5 (GF-5) Advanced Hyperspectral Imager (AHSI) provide more detailed spectral information for geological bodies [11]. Guo et al. [12] utilized the semi-supervised self-learning (SSL) method to investigate the capability of Hyperion data in lithological classification with limited samples. Zhang et al. [13] mapped areas of acicularite and sericite using GF-5 AHSI data and proposed indicators of uranium mineralization. Shebl et al. [14] tested the ability of the PRecursore IperSpettrale della Missione Applicativa hyperspectral data for lithological mapping, and the Support Vector Machine algorithm achieved the highest accuracy. The abundant spectral information makes mineral mapping possible. The extraction of mineral content helps improve the discrimination accuracy of altered rocks [15]. For example, minerals such as hematite, goethite, and jarosite show significant absorption in the visible-to-near-infrared (VNIR) range of 400–1100 nm. Chlorite and clay minerals exhibit absorption peaks in the shortwave infrared (SWIR) range [16]. A combination of hyperspectral and multi-temporal multispectral data has been proven to improve mineral identification accuracy [17]. Liu et al. [18] used Landsat TM data from the same season in four different years to eliminate the random interference of alteration minerals and also used Hyperion data to identify alteration minerals around granites.
Geological mapping using image spectra was often combined with laboratory spectra [19,20]. Using a handheld spectrometer is the easiest way to obtain spectra for rock samples. The spectra obtained with spectrometers are more suitable for studying the physicochemical properties of rocks and avoiding pixel mixing [21]. Wang et al. [22] collected the reflectance spectra of rocks and minerals in Xinjiang, China, using a spectrometer and established a spectral library, which is of great significance for accelerating regional geological research.
Various spectral transformation algorithms and classification algorithms were used to improve the accuracy of rock or mineral extraction from hyperspectral data. Kumar et al. [23] used different spectral enhancement techniques (Principal Component Analysis and Independent Component Analysis) and machine learning algorithms to improve the accuracy of hyperspectral lithology identification, and the results showed that the optimal band based on Joint Mutual Information Maximization had the best classification results. Tan et al. [24] utilized an efficient clustering algorithm based on affinity propagation to select fewer bands with higher lithological discrimination capability. Unlike vegetation and crops, the same kind of rock may have diversified spectral curves (especially different spectral intensities) due to different mineral composition, mineral content, weathering, etc. [25]. High intra-class variability may lead to spectral confusion and poor classification accuracy [26]. Previous studies used the first derivative of spectra combined with the original spectra for classification [27].
Wavelet transform is a powerful method for extracting features of signals, which can obtain features in both time and frequency domains [28]. The number of features remains the same after decomposition and reconstruction. High-frequency components give details in the signal and can highlight abrupt changes in the spectra [29]. Although rock spectra have high intra-class variability, spectral curves of the same type of rock may have similar absorption characteristics [30,31]. Wavelet transform has advantages in extracting subtle features [32,33]. Chen et al. [34] proposed an adaptive wavelet filter for image denoising, which can effectively reduce the variable band noise of different targets. In the area of geological mapping with spectroscopy and hyperspectral sensors, Lorenz et al. [35] found that wavelet-based sparse reduced-rank regression classified geological bodies more accurately than Principal Component Analysis (PCA) and Minimum Noise Fraction (MNF). Feng et al. [36] used Continuous Wavelet Analysis (CWA) to improve the radiometric quality of images by minimizing random noise and enhancing spectral features, with good results in mapping ultramafic rock units. Mitchley et al. [37] improved the accuracy of mineral content estimation via wavelet denoising. In previous studies, low-frequency features were extracted for signal smoothing and filtering [38], while high-frequency features were often considered as noise. Using high-frequency features to improve rock classification accuracy was rarely mentioned.
A random forest classifier is an efficient machine learning algorithm [39]. Random forests were often used for image classification in remote sensing, whether for multispectral data [40], hyperspectral data [41], or synthetic aperture radar data [42]. A random forest classifier can handle datasets with a large number of features and samples and avoid overfitting. Efficient and stable classifiers are necessary for studying the improvement of classification accuracy with wavelet transform. During the evaluation of different decomposition and reconstruction approaches, a random forest classifier was used to classify rocks or rock units. Spectral Angle Mapper (SAM) was widely used in geological mapping, describing the similarity between the target spectrum and a reference spectrum [43]. SAM has also been used to assess intra-class variability [27].
This paper aims to test the ability of wavelet-transformed high-frequency features to improve hyperspectral lithology discrimination. The specific contents include (1) wavelet decomposition of hand-specimen spectra and image spectra and reconstruction of high-frequency coefficients; (2) detecting the intra-class variability of spectral curves before and after wavelet transform; (3) selecting the best mother wavelet function, decomposition level, and reconstruction method for lithology classification; and (4) geological mapping in areas with rock outcrops and less vegetation.

2. Study Area and Data

2.1. Field Rock Samples for Laboratory Spectral Measurements

The rock samples were collected from Huludao City, Liaoning Province, China, in 2013 and 2015. This area is located in the North China Craton, which is an important constituent unit of the continental lithosphere and has experienced multiple phases of regional tectonic evolution. The study area is rich in vegetation types dominated by pine, cypress, and pagoda trees. Rocks from the Mesozoic and Mesoproterozoic eras were collected. Weathered surfaces were removed as much as possible during sample collection. The samples were classified into six types based on their mineralogical composition, including 46 samples of dolostone, 55 samples of andesite, 82 samples of tuff, 109 samples of limestone, 66 samples of sandstone, and 121 samples of granite (Table 1).
Dolostone is composed primarily of dolomite, often mixed with quartz, feldspar, calcite, and clay minerals. Andesite is a kind of neutral ejecta with small amounts of plagioclase feldspar, hornblende, pyroxene, and the porphyry of biotite. Tuff is the most common type of volcaniclastic rock and generally contains a variety of clastic materials. Limestone is a class of carbonate rocks with calcite as the main component. Sandstone is mainly composed of various sand particles cemented together, with a particle diameter of 0.05–2 mm. The cement in the sandstone is dominated by quartz, feldspar, calcite, and clay minerals. Granite is an important constituent of the continental crust and is widely distributed. The major minerals in granite are quartz, potassium feldspar, and plagioclase, while minor minerals are smectite and hornblende.

2.2. Study Area for Lithological Mapping

Study areas A and B for hyperspectral image classification were selected in the eastern Tianshan region of Xinjiang, China. There are many rock outcrops and little vegetation in the areas, making them suitable for rock unit classification using remote sensing techniques [44]. The reason for choosing two study areas with different elevations and surface relief was to test whether the ability of the discrete wavelet transform to improve lithological mapping is affected by surface morphology. Compared to study area B, the elevation in study area A is higher, and study area A has a high degree of topographic relief. A large area of moyite is exposed in the southwest, and the rock types in study area A are mainly granite, granodiorite, and volcanic rock. The elevation and topographic relief are lower in study area B. Alluvium and fluvial deposits, granite, volcanic rocks, and continental clastic rocks are mainly composed of the northern part. Samples from field surveys are marked with red circles in Figure 1. To verify the accuracy of the geological map, we conducted field validation in some areas, which helped in the selection of image samples in the next step. For field samples, we set more sample points in complex rock units or rock groups and fewer sample points in rocks with obvious remote sensing interpretation signs, such as granite. Rock units Qp2W and Qp3X in the northern part of study area B have complex compositions (sandstones, loess, alluvial deposits, etc.), and we made intensive field observations. The geological maps of the two study areas are shown in Figure 2.

2.3. Gaofen-5

GF-5 is one of the satellites in China’s High-Resolution Earth Observation System. GF-5 carries six types of payloads, and the payload used in this study is AHSI (Table 2). There are 150 bands in VNIR and 180 bands in SWIR. The spatial resolution of the data is 30 m. The data acquired on 29 October 2019 (study area A) and 31 May 2019 (study area B) were selected. The full width half max (FWHM) was shown in Figure 3.

3. Methods

Figure 4 shows the workflow of this study. The method primarily consists of four parts: (1) spectral measurement of rocks and preprocessing of GF-5 AHSI data; (2) discrete wavelet transform of laboratory spectra and image spectra; (3) lithology classification and intra-class variability measurement; and (4) accuracy evaluation.

3.1. Laboratory Spectral Processing

3.1.1. Spectral Measurement

The spectra of rocks were measured in the laboratory using the Analytical Spectral Device (ASD) FieldSpec 3. The spectrum ranges from 350 to 2500 nm. The sampling interval is 1.4 nm for 350–1000 nm and 2 nm for 1000–2500 nm. The final spectra were interpolated and output at 1 nm intervals. The distance between the rock sample and the probe was approximately 20 cm during the measurement. This distance is comparable to the size of the rocks. Calibration was performed every ten minutes using a whiteboard. During the spectral measurement, tests were taken on fresh surfaces of rocks, and the average of ten tests was calculated. This helped to avoid errors caused by uneven rock surfaces.

3.1.2. One-Dimensional Discrete Wavelet Transform

Discrete wavelet transform (DWT) is a discretized sampling of the scale and translation parameters for the continuous wavelet transform. DWT can extract the detail and approximation features of signals (Figure 5). One-dimensional (1D) discrete wavelet decomposition and reconstruction were performed on all spectra. At each decomposition level, only the low-frequency coefficients of the signal were decomposed. Daubechies (db) series wavelets have good time and frequency localization properties [46]. The mother wavelet functions (MF) used in this study include db1 (haar), db2, db4, db8, and db10. For dbN wavelets, N is the vanishing moment, and the support length is 2N-1. When the vanishing moment is high, the wavelet coefficients of the high-frequency sub-bands are small, and the computational complexity of the wavelet transform is high. The wavelet function and the scaling function are shown in Figure 6. Table 3 shows four different decomposition–reconstruction methods, in which each MF performs all decomposition–reconstruction methods. The naming abbreviation rule is MF_decomposition level_high-frequency coefficient level used during reconstruction. Finally, twenty sets of high-frequency features were obtained. DWT was implemented in the wavelet toolbox of Matlab R2023a, which was developed by MathWorks Inc. (Natick, MA, USA) [47,48].

3.1.3. Random Forest Classification

To retain as much spectral information as possible, we did not resample the spectral curves. The basic unit of the random forest (RF) algorithm is the decision tree. Samples are fed into each tree for classification. The classification results of several weak classifiers are voted on and selected to form a strong classifier. RF has advantages in processing high-dimensional data and can provide feature importance rankings during classification [49]. RF was implemented in R 4.2.3 [50], with mtry and ntree as two important parameters. Ntree is the number of trees to grow. Mtry controls how many features are available to be considered for each new split. Mtry was often set to the square root of the number of features, and in the classification of laboratory spectra, it was set to 46. Ntree was set to 1000, which was commonly considered a stable value.
To avoid overfitting and underfitting, five-fold cross-validation was used in the classification [51]. Each type of rock sample was divided into five equal parts, with four parts used for training the model and one part used for model validation. The overall accuracy (OA), Kappa coefficient, and F1-Score after cross-validation were used to assess the classification accuracy.

3.2. Image Spectral Processing

3.2.1. GF-5 AHSI Data Preprocessing

GF-5 AHSI data need to be preprocessed before extracting the image spectra, and the preprocessing was implemented in Pixel Information Expert-Hyp (PIE-Hyp) of PIESAT Information Technology Co., Ltd., Beijing, China (Figure 7). First, the VNIR and SWIR bands were stacked. We used the Stripe Removal module of the software to reduce the stripes in the image. We selected the column with better imaging quality, which was considered to be free of stripes, calculated its mean and variance, and used this column as a reference to adjust the pixel values of other columns. The oxygen absorption peak at 762 nm was used to detect the smile effect [52]. The need for smile correction was determined by observing whether the difference images of the VNIR Band 87 and VNIR Band 89 had a significant brightness gradient. The GF-5 data used in this study were less affected by the smile effect and did not require smile correction. Next, digital numbers were calibrated to the top of the atmosphere reflectance, and the 6S radiative transfer model was used to generate surface reflectance products. Orthorectification was used to remove the effects of terrain. The resampling method was set to the nearest neighbor method, and the output data had a spatial resolution of 30 m. Finally, the bands affected by atmospheric water absorption and those with significant noise were excluded. Table 4 shows the excluded bands.
The samples were selected with reference to the geological map. We randomly selected 300 pixels for each type of rock unit as samples. The pixel spectra of the rock units after preprocessing are shown in Figure 8. Pξγ and (E3-N1)t have the highest spectral intensity. Compared to study area A, the distribution of image spectra of rock units in study area B is more dispersed. In the range of 1000–2400 nm, the spectral intensities of rock units are significantly different in study area B. For study area A, C1x, C1gd, P1aer, Pγδ, and Pγδπ have overlapping curves that are difficult to distinguish.

3.2.2. One-Dimensional Discrete Wavelet Transform

For a 297-dimensional hyperspectral image, the spectrum of each pixel can be extracted. The extracted spectral curves have 297 bands. For each pixel in the image, a 1D DWT was performed using the same decomposition–reconstruction methods as shown in Table 3.

3.2.3. Random Forest Classification and Accuracy Assessment

The algorithm for geological mapping was the RF algorithm. The number of ntrees was set to 1000, and mtry was set to 17. Similar to laboratory spectral classification, five-fold cross-validation was used for image spectral classification. OA, Kappa coefficient, and F1-Score were also used to evaluate the accuracy. Based on these indicators, the best MF, decomposition level, and reconstruction method were selected. Finally, all pixels in the image were classified. The results using the high-frequency features were compared to those classified using the original features.

3.3. Evaluation of Intra-Class Variability

SAM was used to evaluate spectral variability with the following equation [53]:
α = cos 1 x 1 T x 2 x 2 T x 2 x 1 T x 1
where x1 or x2 denotes the laboratory spectrum or image spectrum.
The principle of SAM is to consider each spectrum as a high-dimensional vector and measure the similarity of spectra by calculating the angle between two vectors. For the same type of rock, the intra-class SAM between each rock sample and all other rock samples were calculated. The mean and standard deviation of SAM were used to describe intra-class variability. If the mean value of SAM is small, it indicates that the intra-class variability of this rock type is low. The curves after DWT need to be normalized before calculating the intra-class SAM. The following normalization equation was used for the curves after DWT.
d = d d min d max d min
where d denotes laboratory spectrum or image spectrum after DWT. When d denotes the laboratory spectrum, dmax and dmin represent the maximum and minimum values of all laboratory spectra after DWT. When d denotes the image spectrum, dmax and dmin represent the maximum and minimum values of all image spectra after DWT.
The Pearson correlation coefficient was used to describe the relationship between classification accuracy and intra-class variability [54].

4. Results

4.1. Laboratory Spectral Classification

4.1.1. DWT of Spectral Curves

Figure 9 shows the original spectral curves and reconstructed curves of various rocks. The reflectance ranges from 0.15 to 0.4. Granite has high reflectance in the 500–1500 nm range. The reflectance curve has a high positive slope in the 0.4–0.6 μm range and exhibits an iron absorption peak near 1 μm and water absorption peaks at 1.4 μm and 1.9 μm. The absorption feature near 2.2 μm may be due to clay alteration. Andesite has low overall reflectance and water absorption peaks at 1.4 μm and 1.9 μm. Limestone exhibits an absorption peak at 2.3 μm, which is different from other rocks [55]. In addition to water absorption features, the average reflectance spectrum of limestone has a carbonate vibration absorption peak at 2.3 μm. Sandstone exhibits absorption peaks at 1.4 μm, 1.9 μm, and 2.2 μm. In the reconstructed curves, intervals with flat spectral curves (e.g., 1500–1800 nm) are reconstructed close to zero. The position where the curve undergoes a sudden change forms a peak in the reconstructed curve.
Figure 10a shows the spectral curves and the average curve of all granite samples. It shows high intra-class variability in rocks, which may lead to poor classification accuracy. Figure 10b shows the high-frequency curve after DWT (haar_4_1234). Compared with Figure 10a, the intra-class variability decreases, and the curve distribution is more concentrated.

4.1.2. Intra-Class Variability

We calculated the mean (Figure 11) and standard deviation (Figure 12) of intra-class SAM using different decomposition–reconstruction methods. If the SAM angle is greater than 0.1, then generally the spectra are considered different. The mean is greater than 0.2, which suggests that there is significant variability in the class or that there are multiple classes. The mean value of intra-class SAM for all rocks is greater than 0.1. After wavelet decomposition and reconstruction, regardless of the choice of MF, the mean and standard deviation of SAM decreased. This indicates that after DWT, the curve distribution of the same type of rock becomes more concentrated.
It can be observed that for limestone, original spectra have the highest mean of SAM, and the use of three different wavelet decomposition–reconstruction methods (db4_4_1234, db8_4_1234, and db10_4_1234) all show the lowest mean of intra-class SAM values. This may indicate that the spectral curves of the limestone samples have different intensities but similar absorption peaks. The high-frequency features of granite samples show the lowest standard deviation of intra-class SAM values compared to the original spectra, regardless of which MF was used. This observation might be related to the spectral distribution of granite samples being more concentrated after DWT.

4.1.3. Classification Results

OA using original features was 0.465 (ranking 17), with a Kappa coefficient of 0.335. Haar_4_1234 has the highest OA of 0.599 (Table 5). Means and standard deviations for 5-fold cross-validation of laboratory spectral classification are shown in Figure A1 in Appendix A. Features decomposed and reconstructed using haar and db2 wavelets were ranked high. The OA of four-level decomposition was better than the three-level decomposition. The MF with smaller vanishing moments performed better than those with larger vanishing moments.
We computed the highest F1-Score for each type of rock (Table 6). The highest F1-Score was from granite at 0.743. The lowest F1-Score was from andesite, at 0.418. The F1-Score of limestone was improved by 0.208 after DWT, which was the highest for all rocks. The best accuracy for all rocks came from the high-frequency features produced by the db2 wavelet and the haar wavelet. F1-Scores using the original features ranked low. Limestone spectra transformed by DWT obtained better accuracy than the original spectrum, no matter what kind of decomposition and reconstruction method was used.
We tested the overall accuracy of haar_4_1234 using different values of ntree and mtry, which helped to find suitable parameters for laboratory spectra (Figure 13). When the range of the ntree was 300–900 and the range of the mtry was 35–50, the classification accuracy was relatively good. High ntree values increased the training time and might also lead to low classification accuracy. In our experiments, the model achieved the highest overall accuracy of 0.602 when the ntree was 500 and the mtry was 50 or when the ntree was 300 and the mtry was 40.

4.2. Hyperspectral Image Classification

4.2.1. Classification Accuracy Statistics

The classification results of high-frequency features showed significant improvement compared to using the original features (Table 7). Means and standard deviations for 5-fold cross-validation of image spectral classification are shown in Figure A2 in Appendix A. For study area A, haar_4_1234 has the highest classification accuracy (OA: 0.769, Kappa coefficient: 0.736), with an 8.6% improvement over the original OA. For study area B, haar_4_1234 also exhibited the highest classification accuracy (OA: 0.885, Kappa coefficient: 0.862), with an improvement of 7.2% compared to the original OA. It is noteworthy that, for study area B, all the OAs of high-frequency features were higher than those of the original features. Combining the two study areas, the high-frequency features generated by the haar wavelet had high classification accuracy. The highest OAs were all from haar_4_1234.
For the highest F1-Score statistics, Pξγ and C2wt achieved the highest scores (Table 8). Both of them had a classification accuracy of over 90%. Four types of rock units in study area A and three types of rock units in study area B (Table 9) had the highest F1-Score from the haar_4_1234. The classification accuracy of igneous rocks was better than that of sedimentary rocks.

4.2.2. Intra-Class Variability

The image spectral curves and the wavelet-transformed curves of the rock unit samples were demonstrated in Figure 14 and Figure 15. Due to the lower spectral resolution of hyperspectral data, the distribution of high-frequency features is not as concentrated as in laboratory spectra. The intra-class SAM of image spectra showed a similar trend to laboratory spectra. It decreased with the use of high-frequency features (Figure 16, Figure 17, Figure 18 and Figure 19). For study area A, haar wavelets reduce intra-class SAM best. For study area B, the db2 wavelet reduces intra-class SAM best for most rock units.

4.2.3. Lithological Mapping

Combining the classification accuracies of image samples, a pixel-by-pixel DWT was applied to study areas using haar_4_1234 and implemented lithological mapping. The classifier for geological mapping was the RF algorithm. Ntree was set to 1000, and mtry was set to 17. For study area A (Figure 20) and study area B (Figure 21), the result maps had a high similarity with the geological map, and most of the rock units were correctly identified. The confusion matrices of the rock classification results for all pixels are shown in Table 10 and Table 11.
Study area A has some speckles in the result due to scour marks. Almost all of the scour marks were misclassified as other rock units. High-frequency features have improved the confusion of PQg with C1y in the southeastern part of study area A. Compared to Figure 2A, 69.94% of the pixels were correctly classified, which was higher than the 63.76% using the original features. For PQg, 86.6% of the pixels were correctly classified. With the use of high-frequency features, more pixels were correctly classified for all rock units.
For Figure 21a, C2wt was confused with Cξγ in the classification results (western part of study area B), and some of Cγ were misclassified as (E3-N1)t (central and eastern part of study area B). Figure 21b had less noise and better matched the boundaries of rock units with the geological map. It is clear that, consistent with the results of the original spectra, lithological mapping using wavelet high-frequency features had more pixels correctly classified. Compared to Figure 2B, 84.13% of the pixels are correctly classified, which is higher than the 77.34% using the original features. For (E3-N1)t, 89.2% of the pixels were correctly classified. With the use of high-frequency features, the percentage of correctly classified pixels was over 80% for all rock units.
Table 10. The confusion matrix of the rock classification results for all pixels in study area A. (The number before parentheses represents the predicted or actual percentage using high-frequency features of haar_4_1234. The number in parentheses represents the predicted or actual percentage using original features.) The ground truth came from the geological map (Figure 2A).
Table 10. The confusion matrix of the rock classification results for all pixels in study area A. (The number before parentheses represents the predicted or actual percentage using high-frequency features of haar_4_1234. The number in parentheses represents the predicted or actual percentage using original features.) The ground truth came from the geological map (Figure 2A).
Ground Truth (Percent)
C1xC1gdC1yPξγPQgP1aerPγδPγδπ
Predicted (Percent)C1x80.1 (73.2)0.4 (0.1)7.3 (6.1)0.8 (5.3)0.7 (0.5)2.6 (3.4)0.6 (17.3)16.6 (1.7)
C1gd0.4 (0.3)80.2 (66.4)5.1 (9.4)2.2 (1.4)0.3 (0.8)1.1 (5.6)14.2 (0.5)2.1 (18.0)
C1y6.0 (1.2)0.7 (12.3)48.8 (43.1)6.6 (4.9)1.4 (1.6)7.9 (9.0)2.8 (3.8)5.5 (7.2)
Pξγ1.7 (4.0)2.3 (2.1)9.4 (9.5)76.4 (66.8)8.2 (8.8)7.8 (14.1)8.0 (8.8)1.6 (1.7)
PQg1.6 (0.3)0.3 (1.9)5.3 (5.7)3.4 (8.2)86.6 (85.0)5.7 (6.3)2.9 (3.1)0.9 (1.1)
P1aer1.0 (1.8)0.1 (2.8)6.5 (7.1)6.0 (5.9)2.4 (1.7)70.0 (50.5)1.3 (3.8)2.9 (5.8)
Pγδ0.9 (18.5)15.1 (0.5)11.2 (12.6)4.7 (7.2)0.4 (1.3)4.0 (7.4)70.0 (61.4)3.8 (3.6)
Pγδπ8.3 (0.6)1.0 (13.9)6.4 (6.3)0.1 (0.4)0.0 (0.2)0.9 (3.7)0.3 (1.5)66.7 (61.0)
OA = 69.94% (63.76%)Kappa = 0.62 (0.55)
Table 11. The confusion matrix of the rock classification results for all pixels in study area B. (The number before parentheses represents the predicted or actual percentage using high-frequency features of haar_4_1234. The number in parentheses represents the predicted or actual percentage using original features.) The ground truth came from the geological map (Figure 2B).
Table 11. The confusion matrix of the rock classification results for all pixels in study area B. (The number before parentheses represents the predicted or actual percentage using high-frequency features of haar_4_1234. The number in parentheses represents the predicted or actual percentage using original features.) The ground truth came from the geological map (Figure 2B).
Ground Truth (Percent)
(E3-N1)tQp2WQp3XCξγC2wt
Predicted (Percent)(E3-N1)t89.2 (79.8)13.7 (16.6)5.0 (7.2)1.6 (5.8)0.2 (0.1)0.3 (0.4)
Qp2W8.6 (15.2)81.8 (75.2)4.4 (4.9)0.8 (3.5)0.8 (2.1)0.4 (1.3)
Qp3X1.4 (1.8)3.5 (4.7)81.6 (71.8)3.7 (3.6)4.0 (10.2)1.6 (3.5)
0.9 (3.2)1.0 (3.3)4.7 (9.5)84.9 (79.8)3.7 (4.3)8.1 (7.8)
Cξγ0.0 (0.0)0.0 (0.2)3.8 (5.7)1.3 (1.3)87.0 (77.4)3.0 (5.4)
C2wt0.0 (0.0)0.0 (0.1)0.5 (0.8)7.8 (6.0)4.3 (6.0)86.5 (81.6)
OA = 84.13% (77.34%)Kappa = 0.80 (0.71)

4.3. Relationship between OA and Intra-Class Variability

Figure 22a shows the scatterplot of OA and the mean of SAM, and Figure 22b shows the scatterplot of F1-Score and the mean of SAM. We calculated the Pearson correlation coefficient between OA, or F1-Score, and the mean of intra-class SAM (Table 12 and Table 13). The values of SAM are from Figure 11, Figure 16 and Figure 18. For both laboratory spectra and image spectra, there was a negative correlation between intra-class SAM and OA. As the intra-class SAM decreased, there was a tendency for OA to increase. For study area A, the Pearson correlation coefficient between the mean intra-class SAM and OA shows a weak correlation. For dolostone, the Pearson correlation coefficient between the mean intra-class SAM and F1-Score showed a strong correlation at −0.971. This may indicate that the most important factor affecting the accuracy of dolostone classification is intra-class variability. For sandstone and C1gd, F1-Score and SAM showed a weak positive correlation. This suggests that intra-class variability may not be the cause of classification accuracy for these two rocks (rock units).

5. Discussion

Lithological mapping using remote sensing technology presents significant challenges. This work applied the high-frequency components of DWT to hyperspectral lithological classification. The high-frequency features reduced the intra-class variability of rock spectra and highlighted the bands that were beneficial for rock classification.

5.1. Laboratory Spectral Classification

The high-frequency features highlight the differences between rock samples. Although the inter-class separability decreases in most bands, the dominant bands for rock classification are highlighted. The RF classifier has advantages in processing high-dimensional data, makes feature reduction unnecessary, and makes the DWT-RF method highly generalizable. The tuff samples include breccia-crystalline clastic tuff, andesite tuff, andesite-bearing breccia-crystalline clastic tuff, etc. (Figure 23). The abundance of subclasses in the samples may affect the classification accuracy.
In our study, we only performed wavelet decompositions at three or four levels. As the number of decomposition levels increases, more high-frequency information is extracted, and the time for wavelet decomposition and reconstruction increases. Wavelet transform only decomposes low-frequency coefficients, so higher-level low-frequency coefficients may not contain useful information for classification. Wavelet packet decomposition is a further optimization of wavelet decomposition. The wavelet packet decomposition decomposes not only low-frequency sub-bands but also high-frequency sub-bands at each level of signal decomposition. We consider using the wavelet packet transform to extract more features to improve the classification in the future.

5.2. Hyperspectral Image Classification

Lithological mapping using GF-5 AHSI data obtained clear boundaries for rock units. If visual interpretation is assisted after classification, the boundaries can be directly used for producing geological maps. The image spectra are considered more suitable as samples for lithological mapping, and laboratory spectra are somewhat different from image spectra due to the influence of rock morphology [56]. C1y has the poorest classification accuracy of all the rock units, with only 48.8% of the image elements correctly classified, even after using high-frequency features. In Figure 14e, the image spectra of sample A have high intra-class variability, with spectral intensities varying by more than 0.2. The high intra-class variability in image spectra may be due to rock fragmentation or weathering [57]. In addition, the black powder formed via coal weathering also interferes with lithology discrimination (Figure 24) [44]. Unlike handheld spectrometers, hyperspectral sensors acquire spectra of rock weathering surfaces and may mix multiple types of rocks in one pixel [58].

5.3. Limitations and Strengths of This Study

The DWT-RF method demonstrates good classification ability in both laboratory spectra and image spectra, indicating that subtle differences between rocks are captured by the high-frequency features. OA increased for all fourteen types of rock units, and high-frequency features consistently improved the effectiveness of lithological discrimination. In the laboratory spectral classification of rock samples, granite and dolostone were better identified, but andesite had an F1-Score of only 0.418. After using the high-frequency features, the boundaries of Pξγ in the southern part and C1x in the northwestern part of study area A were well extracted. Pγδπ and P1aer in the northern part of study area A still had pixel confusion after the use of high-frequency features. Study area B with less surface relief had better classification results.
Currently, this method can only be used for geological mapping in arid or semi-arid areas with exposed rock outcrops. In vegetation-covered areas, the image spectrum is a mixture of vegetation and rocks, which is different from the laboratory spectrum [59]. Hewson et al. [60] used a threshold of 0.4 to mask the vegetation before generating the mineral composition map. A variety of rocks are exposed in the area where the hand specimens were collected, and many of them, such as dolostone and andesite, were not found in the lithological mapping area. Unfortunately, the area where the hand-specimen rocks were collected was heavily vegetated, with few samples with no or little vegetation cover. The area chosen for lithological mapping is the Eastern Tianshan region, which is a good area for geological investigations because there is little vegetation cover. The classification of the laboratory spectra and the classification of the image spectra were tested separately, so we chose the two most suitable areas. When DWT has a positive effect on both the hand-specimen spectra and image spectra, methods are considered essentially feasible. If laboratory spectra could assist in geological mapping, this would increase the accuracy of lithological mapping. Establishing the relationship between vegetation spectra and underlying rock spectra enables geological mapping in vegetation-covered areas, which is beneficial for regional geological mapping. More efforts are needed in the future to improve classification accuracy and reduce noise in hyperspectral data. High-frequency features are susceptible to noise, and filtering may lose information useful for classification. Data fusion from multiple sources may be a good choice, which helps to reduce the interference of noise. In addition, SAM is sometimes sensitive to illumination variations. Finding more efficient classifiers and more accurate spectral similarity metric indices is important for geological mapping.

6. Conclusions

The purpose of this study is to discuss the ability of DWT to improve rock classification. Based on laboratory spectra and GF-5 AHSI data, it can be concluded that high-frequency features after DWT reduce the intra-class variability of all rocks and rock units. For both laboratory spectra and image spectra, four-level DWT using the haar wavelet gives the best classification results. The F1-Score of dolostone improved from 0.464 to 0.689 after DWT. For the classification of rock units, after using high-frequency features, the F1-Score of Pξγ is 0.907, and the F1-Score of C2wt is 0.933. The classification accuracy is influenced by the MF, decomposition levels, and reconstruction methods. MF with smaller vanishing moments can better improve classification accuracy. After DWT, the intra-class SAM reduces, and the mean of SAM is negatively correlated with the classification accuracy.
More strategies need to be used to reduce the interference of noise, which is important for hyperspectral remote sensing geology. Once the relationship between vegetation and underlying rocks is determined, regional geological mapping will be possible.

Author Contributions

Conceptualization, S.G.; methodology, S.G. and Q.J.; software, S.G.; validation, S.G. and Q.J.; formal analysis, S.G.; investigation, S.G. and Q.J.; resources, Q.J.; data curation, S.G. and Q.J.; writing—original draft preparation, S.G.; writing—review and editing, S.G. and Q.J.; visualization, S.G.; supervision, Q.J.; project administration, Q.J.; funding acquisition, Q.J. All authors have read and agreed to the published version of the manuscript.

Funding

This study is supported by the China Geological Survey Project (Project Number: DD20191011).

Data Availability Statement

The authors do not have permission to share data. GF-5 AHSI data can be accessed through the website https://data.cresda.cn/#/2dMap (accessed on 20 June 2023) and require contacting the Chinese Resource Satellite Application Center.

Acknowledgments

The authors thank the students at Jilin University for their work in the field and laboratory.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. Means and standard deviations for 5-fold cross-validation of laboratory spectral classification.
Figure A1. Means and standard deviations for 5-fold cross-validation of laboratory spectral classification.
Remotesensing 15 05334 g0a1
Figure A2. Means and standard deviations for 5-fold cross-validation of image spectral classification.
Figure A2. Means and standard deviations for 5-fold cross-validation of image spectral classification.
Remotesensing 15 05334 g0a2

References

  1. Radford, D.D.G.; Cracknell, M.J.; Roach, M.J.; Cumming, G.V. Geological Mapping in Western Tasmania Using Radar and Random Forests. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2018, 11, 3075–3087. [Google Scholar] [CrossRef]
  2. Rezaei, A.; Hassani, H.; Moarefvand, P.; Golmohammadi, A. Lithological mapping in Sangan region in Northeast Iran using ASTER satellite data and image processing methods. Geol. Ecol. Landsc. 2020, 4, 59–70. [Google Scholar] [CrossRef]
  3. El-Omairi, M.A.; El Garouani, A. A review on advancements in lithological mapping utilizing machine learning algorithms and remote sensing data. Heliyon 2023, 9, e20168. [Google Scholar] [CrossRef] [PubMed]
  4. Windeler, D.S. Garnet-pyroxene alteration mapping in the Ludwig skarn (Yerington, Nevada) with geoscan airborne multispectral data. Photogramm. Eng. Remote Sens. 1993, 59, 1277–1286. [Google Scholar]
  5. Greenbaum, D. Lithological discrimination in central Snowdonia using airborne multispectral scanner imagery. Int. J. Remote Sens. 1987, 8, 799–816. [Google Scholar] [CrossRef]
  6. Nguemhe Fils, S.C.; Bekele Mongo, C.H.; Nkouathio, D.G.; Mimba, M.E.; Etouna, J.; Njandjock Nouck, P.; Nyeck, B. Radarsat-1 image processing for regional-scale geological mapping with mining vocation under dense vegetation and equatorial climate environment, Southwestern Cameroon. Egypt. J. Remote Sens. Space Sci. 2018, 21, S43–S54. [Google Scholar] [CrossRef]
  7. Bentahar, I.; Raji, M. Comparison of Landsat OLI, ASTER, and Sentinel 2A data in lithological mapping: A Case study of Rich area (Central High Atlas, Morocco). Adv. Space Res. 2021, 67, 945–963. [Google Scholar] [CrossRef]
  8. Alibegovic, G.; Schut, A.G.T.; Wardell-Johnson, G.W.; Robinson, T.P. Seasonal differences assist in mapping granite outcrops using Landsat TM imagery across the Southwest Australian Floristic Region. J. Spat. Sci. 2015, 60, 37–49. [Google Scholar] [CrossRef]
  9. Pal, M.K.; Rasmussen, T.M.; Abdolmaleki, M. Multiple multi-spectral remote sensing data fusion and integration for geological mapping. In Proceedings of the 2019 10th Workshop on Hyperspectral Imaging and Signal Processing: Evolution in Remote Sensing (WHISPERS), Amsterdam, The Netherlands, 24–26 September 2019; pp. 1–5. [Google Scholar]
  10. Rowan, L.C.; Mars, J.C.; Simpson, C.J. Lithologic mapping of the Mordor, NT, Australia ultramafic complex by using the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER). Remote Sens. Environ. 2005, 99, 105–126. [Google Scholar] [CrossRef]
  11. Thannoun, R.G. Mapping lithological and mineralogical units using hyperspectral imagery. Malays. J. Sci. 2021, 40, 93–106. [Google Scholar] [CrossRef]
  12. Guo, X.; Li, P.; Li, J. Lithological mapping using EO-1 Hyperion hyperspectral data and semisupervised self-learning method. J. Appl. Remote Sens. 2021, 15, 032209. [Google Scholar] [CrossRef]
  13. Zhang, Y.; Pan, W.; Yu, Z. Application of Gaofen-5 hyperspectral data in uranium exploration: A case study of Weijing in Inner Mongolia, China. In Proceedings of the SPIE, Ninth Symposium on Novel Photoelectronic Detection Technology and Applications, Hefei, China, 2–4 November 2022; Volume 12617. [Google Scholar]
  14. Shebl, A.; Abriha, D.; Fahil, A.S.; El-Dokouny, H.A.; Elrasheed, A.A.; Csámer, Á. PRISMA hyperspectral data for lithological mapping in the Egyptian Eastern Desert: Evaluating the support vector machine, random forest, and XG boost machine learning algorithms. Ore Geol. Rev. 2023, 161, 105652. [Google Scholar] [CrossRef]
  15. Peyghambari, S.; Zhang, Y. Hyperspectral remote sensing in lithological mapping, mineral exploration, and environmental geology: An updated review. J. Appl. Remote Sens. 2021, 15, 031501. [Google Scholar] [CrossRef]
  16. Pour, A.B.; Hashim, M. ASTER, ALI and Hyperion sensors data for lithological mapping and ore minerals exploration. SpringerPlus 2014, 3, 130. [Google Scholar] [CrossRef]
  17. Hu, B.; Xu, Y.; Wan, B.; Wu, X.; Yi, G. Hydrothermally altered mineral mapping using synthetic application of Sentinel-2A MSI, ASTER and Hyperion data in the Duolong area, Tibetan Plateau, China. Ore Geol. Rev. 2018, 101, 384–397. [Google Scholar] [CrossRef]
  18. Liu, L.; Zhou, J.; Han, L.; Xu, X. Mineral mapping and ore prospecting using Landsat TM and Hyperion data, Wushitala, Xinjiang, northwestern China. Ore Geol. Rev. 2017, 81, 280–295. [Google Scholar] [CrossRef]
  19. Chen, X.; Warner, T.A.; Campagna, D.J. Integrating visible, near-infrared and short-wave infrared hyperspectral and multispectral thermal imagery for geological mapping at Cuprite, Nevada: A rule-based system. Int. J. Remote Sens. 2010, 31, 1733–1752. [Google Scholar] [CrossRef]
  20. El Sobky, M.A.; Madani, A.A.; Surour, A.A. Spectral characterization of the Batuga granite pluton, South Eastern Desert, Egypt: Influence of lithological and mineralogical variation on ASD Terraspec data. Arab. J. Geosci. 2020, 13, 1246. [Google Scholar] [CrossRef]
  21. Chiari, R.; Longhi, I.; Sgavetti, M.A.; Gelli, A.; Orsi, A.; Pecoraro, F. Spectral classification of rocks: Analysis of laboratory 0.4- to 2.5-um reflectance and 2.5- to 25-um transmittance spectra of sedimentary and metamorphic rocks. In Proceedings of the Earth Surface Remote Sensing, London, UK, 22–25 September 1997; pp. 284–294. [Google Scholar]
  22. Wang, S.; Zhou, K.; Zhang, N.; Wang, J. Spectral data analysis of rock and mineral in Hatu Western Junggar Region, Xinjiang. In Proceedings of the SPIE Asia Pacific Remote Sensing, Beijing, China, 13–16 October 2014. [Google Scholar]
  23. Kumar, C.; Chatterjee, S.; Oommen, T.; Guha, A. Automated lithological mapping by integrating spectral enhancement techniques and machine learning algorithms using AVIRIS-NG hyperspectral data in Gold-bearing granite-greenstone rocks in Hutti, India. Int. J. Appl. Earth Obs. Geoinf. 2020, 86, 102006. [Google Scholar] [CrossRef]
  24. Tan, Y.; Lu, L.; Bruzzone, L.; Guan, R.; Chang, Z.; Yang, C. Hyperspectral Band Selection for Lithologic Discrimination and Geological Mapping. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 471–486. [Google Scholar] [CrossRef]
  25. Murphy, R.J.; Monteiro, S.T.; Schneider, S. Evaluating classification techniques for mapping vertical geology using field-based hyperspectral sensors. IEEE Trans. Geosci. Remote Sens. 2012, 50, 3066–3080. [Google Scholar] [CrossRef]
  26. Sgavetti, M.; Pompilio, L.; Meli, S. Reflectance spectroscopy (0.3–2.5 µm) at various scales for bulk-rock identification. Geosphere 2006, 2, 142–160. [Google Scholar] [CrossRef]
  27. Zhang, X.; Li, P. Lithological mapping from hyperspectral data by improved use of spectral angle mapper. Int. J. Appl. Earth Obs. Geoinf. 2014, 31, 95–109. [Google Scholar] [CrossRef]
  28. Mallat, S.G. A Theory for Multiresolution Signal Decomposition: The Wavelet Representation. IEEE Trans. Pattern Anal. Mach. Intell. 1989, 11, 674–693. [Google Scholar] [CrossRef]
  29. Yang, C.B.; Gao, W.B.; Hou, G.Y.; Li, X.Z.; Gao, M.T. Response Relationship between Feldspar Content and Characteristic Spectra in Igneous Rocks. Guang Pu Xue Yu Guang Pu Fen Xi/Spectrosc. Spectr. Anal. 2019, 39, 2077–2082. [Google Scholar] [CrossRef]
  30. Amer, R.; Kusky, T.; Ghulam, A. Lithological mapping in the Central Eastern Desert of Egypt using ASTER data. J. Afr. Earth Sci. 2010, 56, 75–82. [Google Scholar] [CrossRef]
  31. Hunt, G.R. Near-infrared (1.3–2.4 micrometre) spectra of alteration minerals—Potential for use in remote sensing. Geophysics 1979, 44, 1974–1986. [Google Scholar] [CrossRef]
  32. Guo, S.; Yang, C.; He, R.; Li, Y. Improvement of Lithological Mapping Using Discrete Wavelet Transformation from Sentinel-1 SAR Data. Remote Sens. 2022, 14, 5824. [Google Scholar] [CrossRef]
  33. Yang, C.B.; Liu, N.; Kuai, K.F. Research on Relationship between Spectral Characteristics, Physical Parameters and Metal Elements of Rocks in Xingcheng Area. Guang Pu Xue Yu Guang Pu Fen Xi/Spectrosc. Spectr. Anal. 2019, 39, 2953–2965. [Google Scholar] [CrossRef]
  34. Chen, B.; Feng, X.; Wu, R.; Guo, Q.; Wang, X.; Ge, S. Adaptive Wavelet Filter With Edge Compensation for Remote Sensing Image Denoising. IEEE Access 2019, 7, 91966–91979. [Google Scholar] [CrossRef]
  35. Lorenz, S.; Ghamisi, P.; Kirsch, M.; Jackisch, R.; Rasti, B.; Gloaguen, R. Feature extraction for hyperspectral mineral domain mapping: A test of conventional and innovative methods. Remote Sens. Environ. 2021, 252, 112129. [Google Scholar] [CrossRef]
  36. Feng, J.; Rogge, D.; Rivard, B. Comparison of lithological mapping results from airborne hyperspectral VNIR-SWIR, LWIR and combined data. Int. J. Appl. Earth Obs. Geoinf. 2018, 64, 340–353. [Google Scholar] [CrossRef]
  37. Mitchley, M.; Sears, M.; Damelin, S. Target detection in hyperspectral mineral data using wavelet analysis. In Proceedings of the 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 12–17 July 2009; pp. IV-881–IV-884. [Google Scholar]
  38. Guo, B.; Guo, X.; Zhang, B.; Suo, L.; Bai, H.; Luo, P. Using a Two-Stage Scheme to Map Toxic Metal Distributions Based on GF-5 Satellite Hyperspectral Images at a Northern Chinese Opencast Coal Mine. Remote Sens. 2022, 14, 5804. [Google Scholar] [CrossRef]
  39. Pal, M. Random forest classifier for remote sensing classification. Int. J. Remote Sens. 2005, 26, 217–222. [Google Scholar] [CrossRef]
  40. Xi, Y.; Mohamed Taha, A.M.; Hu, A.; Liu, X. Accuracy comparison of various remote sensing data in lithological classification based on random forest algorithm. Geocarto Int. 2022, 37, 14451–14479. [Google Scholar] [CrossRef]
  41. Puletti, N.; Camarretta, N.; Corona, P. Evaluating EO1-Hyperion capability for mapping conifer and broadleaved forests. Eur. J. Remote Sens. 2016, 49, 157–169. [Google Scholar] [CrossRef]
  42. Zhang, X.; Xu, J.; Chen, Y.; Xu, K.; Wang, D. Coastal Wetland Classification with GF-3 Polarimetric SAR Imagery by Using Object-Oriented Random Forest Algorithm. Sensors 2021, 21, 3395. [Google Scholar] [CrossRef] [PubMed]
  43. Rao, D.A.; Guha, A. Potential Utility of Spectral Angle Mapper and Spectral Information Divergence Methods for mapping lower Vindhyan Rocks and Their Accuracy Assessment with Respect to Conventional Lithological Map in Jharkhand, India. J. Indian Soc. Remote Sens. 2018, 46, 737–747. [Google Scholar] [CrossRef]
  44. Lu, Y.; Yang, C.; He, R. Towards lithology mapping in semi-arid areas using time-series Landsat-8 data. Ore Geol. Rev. 2022, 150, 105163. [Google Scholar] [CrossRef]
  45. Zhang, Y.; Wang, J.; Yu, J.; Tian, J.; Zhou, J. Study of Late Paleozoic Mineralization and Target Area Selection in the Jorotag Metallogenic Belt of the East Tianshan Mountains; Geological Survey Institute of Jilin University: Changchun, China, 2018. [Google Scholar]
  46. Daubechies, I. Orthonormal bases of compactly supported wavelets. Commun. Pure Appl. Math. 1988, 41, 909–996. [Google Scholar] [CrossRef]
  47. MATLAB, 9.14.0.2239454 (R2023a); The MathWorks Inc.: Natick, MA, USA, 2023.
  48. Misiti, Y.; Misiti, M.; Oppenheim, G.; Poggi, J.-M. Micronde: A Matlab Wavelet Toolbox for Signals and Images. In Wavelets and Statistics; Antoniadis, A., Oppenheim, G., Eds.; Springer: New York, NY, USA, 1995; pp. 239–259. [Google Scholar]
  49. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  50. Liaw, A.; Wiener, M. Classification and regression by randomForest. R News 2002, 2, 18–22. [Google Scholar]
  51. Wong, T.T. Performance evaluation of classification algorithms by k-fold and leave-one-out cross validation. Pattern Recognit. 2015, 48, 2839–2846. [Google Scholar] [CrossRef]
  52. Yokoya, N.; Miyamura, N.; Iwasaki, A. Detection and correction of spectral and spatial misregistrations for hyperspectral data. In Proceedings of the 2010 IEEE International Geoscience and Remote Sensing Symposium, Honolulu, HI, USA, 25–30 July 2010; pp. 1003–1006. [Google Scholar]
  53. Kruse, F.A.; Lefkoff, A.B.; Boardman, J.W.; Heidebrecht, K.B.; Shapiro, A.T.; Barloon, P.J.; Goetz, A.F.H. The spectral image processing system (SIPS)-interactive visualization and analysis of imaging spectrometer data. Remote Sens. Environ. 1993, 44, 145–163. [Google Scholar] [CrossRef]
  54. Pearson, K. Notes on the history of correlation. Biometrika 1920, 13, 25–45. [Google Scholar] [CrossRef]
  55. Rasouli Beirami, M.; Tangestani, M.H. A New Band Ratio Approach for Discriminating Calcite and Dolomite by ASTER Imagery in Arid and Semiarid Regions. Nat. Resour. Res. 2020, 29, 2949–2965. [Google Scholar] [CrossRef]
  56. Hecker, C.; Meijde, M.v.d.; Werff, H.v.d.; Meer, F.D. Assessing the Influence of Reference Spectra on Synthetic SAM Classification Results. IEEE Trans. Geosci. Remote Sens. 2008, 46, 4162–4172. [Google Scholar] [CrossRef]
  57. Lu, Y.; Yang, C.; Jiang, Q. Evaluation of the Performance of Time-Series Sentinel-1 Data for Discriminating Rock Units. Remote Sens. 2021, 13, 4824. [Google Scholar] [CrossRef]
  58. van der Meer, F.D.; van der Werff, H.M.; van Ruitenbeek, F.J.; Hecker, C.A.; Bakker, W.H.; Noomen, M.F.; van der Meijde, M.; Carranza, E.J.M.; de Smeth, J.B.; Woldai, T. Multi- and hyperspectral geologic remote sensing: A review. Int. J. Appl. Earth Obs. Geoinf. 2012, 14, 112–128. [Google Scholar] [CrossRef]
  59. Grebby, S.; Cunningham, D.; Tansey, K.; Naden, J. The impact of vegetation on lithological mapping using airborne multispectral data: A case study for the north Troodos region, Cyprus. Remote Sens. 2014, 6, 10860–10887. [Google Scholar] [CrossRef]
  60. Hewson, R.; Mshiu, E.; Hecker, C.; van der Werff, H.; van Ruitenbeek, F.; Alkema, D.; van der Meer, F. The application of day and night time ASTER satellite imagery for geothermal and mineral mapping in East Africa. Int. J. Appl. Earth Obs. Geoinf. 2020, 85, 101991. [Google Scholar] [CrossRef]
Figure 1. The location of study areas (A,B) and the distribution of field samples. The background images are from Gaofen-5 AHSI data© China Centre for Resources Satellite Data and Application (R: VNIR Band 59, G: VNIR Band 38, and B: VNIR Band 20). Elevation data are extracted from the shuttle radar topography mission (SRTM).
Figure 1. The location of study areas (A,B) and the distribution of field samples. The background images are from Gaofen-5 AHSI data© China Centre for Resources Satellite Data and Application (R: VNIR Band 59, G: VNIR Band 38, and B: VNIR Band 20). Elevation data are extracted from the shuttle radar topography mission (SRTM).
Remotesensing 15 05334 g001
Figure 2. Geological maps of the study areas (A,B). The scale of the mapping format is 1:500,000 [45].
Figure 2. Geological maps of the study areas (A,B). The scale of the mapping format is 1:500,000 [45].
Remotesensing 15 05334 g002
Figure 3. The FWHM of GF-5 AHSI data.
Figure 3. The FWHM of GF-5 AHSI data.
Remotesensing 15 05334 g003
Figure 4. The method flowchart of this study.
Figure 4. The method flowchart of this study.
Remotesensing 15 05334 g004
Figure 5. 1D-DWT of spectral signals (An represents the low-frequency coefficients after DWT, Dn represents the high-frequency coefficients after DWT, and n represents the number of decomposition levels).
Figure 5. 1D-DWT of spectral signals (An represents the low-frequency coefficients after DWT, Dn represents the high-frequency coefficients after DWT, and n represents the number of decomposition levels).
Remotesensing 15 05334 g005
Figure 6. Wavelet function and scaling function of (a) haar, (b) db2, (c) db4, (d) db8, and (e) db10.
Figure 6. Wavelet function and scaling function of (a) haar, (b) db2, (c) db4, (d) db8, and (e) db10.
Remotesensing 15 05334 g006
Figure 7. Preprocessing process of GF-5 AHSI data.
Figure 7. Preprocessing process of GF-5 AHSI data.
Remotesensing 15 05334 g007
Figure 8. The mean image spectra for each rock unit in (a) study area A and (b) study area B.
Figure 8. The mean image spectra for each rock unit in (a) study area A and (b) study area B.
Remotesensing 15 05334 g008
Figure 9. Original spectral curves and reconstructed curves of (a) dolostone, (b) andesite, (c) tuff, (d) limestone, (e) sandstone, and (f) granite. Curve db2_3_23 was shifted upward by 0.01, and curve db2_4_234 was shifted upward by 0.02. (Rock spectral curves using the left y-axis and high-frequency features after DWT using the right y-axis.)
Figure 9. Original spectral curves and reconstructed curves of (a) dolostone, (b) andesite, (c) tuff, (d) limestone, (e) sandstone, and (f) granite. Curve db2_3_23 was shifted upward by 0.01, and curve db2_4_234 was shifted upward by 0.02. (Rock spectral curves using the left y-axis and high-frequency features after DWT using the right y-axis.)
Remotesensing 15 05334 g009aRemotesensing 15 05334 g009b
Figure 10. (a) The original laboratory spectral curve and (b) the high-frequency curve after DWT (db2_4_1234) of granite samples. (Spectra were from ASD FieldSpec 3.)
Figure 10. (a) The original laboratory spectral curve and (b) the high-frequency curve after DWT (db2_4_1234) of granite samples. (Spectra were from ASD FieldSpec 3.)
Remotesensing 15 05334 g010
Figure 11. Mean of intra-class SAM in laboratory spectra.
Figure 11. Mean of intra-class SAM in laboratory spectra.
Remotesensing 15 05334 g011
Figure 12. The standard deviation of intra-class SAM in laboratory spectra.
Figure 12. The standard deviation of intra-class SAM in laboratory spectra.
Remotesensing 15 05334 g012
Figure 13. OA corresponding to different ntree and mtry in the RF classifier.
Figure 13. OA corresponding to different ntree and mtry in the RF classifier.
Remotesensing 15 05334 g013
Figure 14. (a,c,e,g,i,k,m,o) are the original spectral curves, and (b,d,f,h,j,l,n,p) are the high-frequency curves after DWT (db2_4_1234) of study area A. (Spectra were from GF-5 AHSI data.)
Figure 14. (a,c,e,g,i,k,m,o) are the original spectral curves, and (b,d,f,h,j,l,n,p) are the high-frequency curves after DWT (db2_4_1234) of study area A. (Spectra were from GF-5 AHSI data.)
Remotesensing 15 05334 g014aRemotesensing 15 05334 g014bRemotesensing 15 05334 g014c
Figure 15. (a,c,e,g,i,k) are the original spectral curves, and (b,d,f,h,j,l) are the high-frequency curves after DWT (db2_4_1234) of study area B. (Spectra were from GF-5 AHSI data.)
Figure 15. (a,c,e,g,i,k) are the original spectral curves, and (b,d,f,h,j,l) are the high-frequency curves after DWT (db2_4_1234) of study area B. (Spectra were from GF-5 AHSI data.)
Remotesensing 15 05334 g015aRemotesensing 15 05334 g015b
Figure 16. Mean values of intra-class SAM in study area A.
Figure 16. Mean values of intra-class SAM in study area A.
Remotesensing 15 05334 g016
Figure 17. The standard deviation of intra-class SAM in study area A.
Figure 17. The standard deviation of intra-class SAM in study area A.
Remotesensing 15 05334 g017
Figure 18. Mean values of intra-class SAM in study area B.
Figure 18. Mean values of intra-class SAM in study area B.
Remotesensing 15 05334 g018
Figure 19. The standard deviation of intra-class SAM in study area B.
Figure 19. The standard deviation of intra-class SAM in study area B.
Remotesensing 15 05334 g019
Figure 20. Lithological mapping results of study area A using (a) original spectra and (b) high-frequency features (haar_4_1234), with blue boxes indicating locations with significant differences in results.
Figure 20. Lithological mapping results of study area A using (a) original spectra and (b) high-frequency features (haar_4_1234), with blue boxes indicating locations with significant differences in results.
Remotesensing 15 05334 g020
Figure 21. Lithological mapping results of study area B using (a) original spectra and (b) high-frequency features (haar_4_1234), with blue boxes indicating locations with significant differences in results.
Figure 21. Lithological mapping results of study area B using (a) original spectra and (b) high-frequency features (haar_4_1234), with blue boxes indicating locations with significant differences in results.
Remotesensing 15 05334 g021
Figure 22. (a) The scatterplot of OA and the mean of SAM. (b) The scatterplot of F1-Score and the mean of SAM. (The lines in the figure are linear fitting lines.)
Figure 22. (a) The scatterplot of OA and the mean of SAM. (b) The scatterplot of F1-Score and the mean of SAM. (The lines in the figure are linear fitting lines.)
Remotesensing 15 05334 g022
Figure 23. Part of the collected tuff samples. (a) Breccia-crystalline clastic tuff, (b) andesite tuff, (c) breccia-crystalline clastic tuff, (d) and andesite-bearing breccia-crystalline clastic tuff.
Figure 23. Part of the collected tuff samples. (a) Breccia-crystalline clastic tuff, (b) andesite tuff, (c) breccia-crystalline clastic tuff, (d) and andesite-bearing breccia-crystalline clastic tuff.
Remotesensing 15 05334 g023
Figure 24. Black powder formed via coal weathering was found in some areas [44].
Figure 24. Black powder formed via coal weathering was found in some areas [44].
Remotesensing 15 05334 g024
Table 1. Names and photos of rock samples collected in the field.
Table 1. Names and photos of rock samples collected in the field.
NamePhotoNamePhoto
dolostoneRemotesensing 15 05334 i001andesiteRemotesensing 15 05334 i002
tuffRemotesensing 15 05334 i003limestoneRemotesensing 15 05334 i004
sandstoneRemotesensing 15 05334 i005graniteRemotesensing 15 05334 i006
Table 2. The parameters of GF-5 AHSI.
Table 2. The parameters of GF-5 AHSI.
ParametersGF-5 AHSI
Orbit altitude705 km
Swath width60 km
Spatial resolution30 m
Spectral resolutionVNIR: 5 nm; SWIR: 10 nm
Number of bandsVNIR: 150; SWIR: 180
Spectral range0.39–2.51 μm
SWIR Signal-to-Noise Ratio (SNR)~500
Dispersive systemsGrating
Table 3. Four decomposition–reconstruction methods and abbreviations were used in the study.
Table 3. Four decomposition–reconstruction methods and abbreviations were used in the study.
Decomposition–Reconstruction MethodAbbreviation (Using Haar Wavelets as an Example)
Reconstruct the high-frequency coefficients of the second and third levels after the three-level decompositionhaar_3_23
Reconstruct all three high-frequency coefficients after three-level decompositionhaar_3_123
Reconstruct the high-frequency coefficients of the second, third, and fourth levels after four-level decompositionhaar_4_234
Reconstruct all four high-frequency coefficients after four-level decompositionhaar_4_1234
Table 4. Bands removed in GF-5 AHSI data preprocessing.
Table 4. Bands removed in GF-5 AHSI data preprocessing.
Wavelength RangeBand NumberWavelength
VNIRVNIR 1–4390–403 nm
SWIRSWIR 1–41004–1030 nm
SWIR 43–501359–1418 nm
SWIR 96–1121805–1940 nm
Table 5. OA and Kappa coefficients of laboratory spectra.
Table 5. OA and Kappa coefficients of laboratory spectra.
FeaturesOAKappa Coefficients
haar_4_12340.5990.502
haar_4_2340.5950.496
db2_4_2340.5930.493
haar_3_230.5800.476
haar_3_1230.5740.467
db2_4_12340.5720.465
db4_4_12340.5560.441
db4_4_2340.5560.443
db2_3_230.5510.435
db10_4_2340.5350.413
db2_3_1230.5300.405
db10_4_12340.5260.403
db8_4_2340.5140.384
db8_4_12340.5110.380
db4_3_230.4890.346
db4_3_1230.4740.324
original spectrum0.4650.335
db10_3_230.4430.277
db10_3_1230.4320.262
db8_3_1230.4300.260
db8_3_230.4280.257
Table 6. The highest value of F1-Score among the rock samples on laboratory spectra using high-frequency features.
Table 6. The highest value of F1-Score among the rock samples on laboratory spectra using high-frequency features.
Rock TypesHighest F1-ScoreFeatures
(Highest F1-Score)
F1-Score of Haar_4_1234Haar_4_1234 RankingOriginal F1-ScoreOriginal F1-Score Ranking
Dolostone0.689db2_4_2340.65570.46418
Andesite0.418db2_4_2340.38320.2719
Tuff0.519db2_3_1230.46240.39818
Limestone0.665Haar_4_12340.66510.45721
Sandstone0.465Haar_4_2340.45620.2708
Granite0.743Haar_4_2340.74020.64316
Table 7. OA and Kappa coefficient of study area A and study area B. The order of features was sorted from high to low in OA.
Table 7. OA and Kappa coefficient of study area A and study area B. The order of features was sorted from high to low in OA.
Study Area AStudy Area B
FeaturesOAKappa CoefficientsFeaturesOAKappa Coefficients
haar_4_12340.7690.736haar_4_12340.8850.862
haar_4_2340.7680.735haar_3_1230.8840.861
haar_3_1230.7660.733haar_4_2340.8800.856
haar_3_230.7600.725haar_3_230.8780.853
db4_4_2340.7210.681db2_4_2340.8700.844
db2_3_1230.7200.68db2_4_12340.8680.842
db2_4_12340.7200.68db10_4_12340.8580.829
db2_4_2340.7180.677db2_3_1230.8570.828
db2_3_230.7170.676db2_3_230.8540.825
db10_4_12340.7090.668db8_4_12340.8530.823
db10_4_2340.7090.667db4_3_1230.8490.819
db4_4_12340.7090.667db10_4_2340.8490.819
db8_4_12340.7000.658db8_4_2340.8490.819
db4_3_230.7000.657db4_4_12340.8490.819
db8_4_2340.6990.656db4_3_230.8480.818
db4_3_1230.6980.655db4_4_2340.8480.817
db8_3_1230.6860.641db8_3_1230.8390.807
db8_3_230.6850.64db10_3_1230.8360.803
original0.6830.637db8_3_230.8330.799
db10_3_1230.6790.633db10_3_230.8330.799
db10_3_230.6620.614original0.8130.776
Table 8. The highest value of F1-Score among the rock units on the image spectra in study area A using high-frequency features.
Table 8. The highest value of F1-Score among the rock units on the image spectra in study area A using high-frequency features.
Rock UnitsHighest F1-ScoreFeatures
(Highest F1-Score)
F1-Score of Haar_4_1234Haar_4_1234 RankingOriginal F1-ScoreOriginal F1-Score Ranking
C1x0.745haar_4_2340.73720.62721
C1gd0.720haar_3_1230.67930.6437
C1y0.802haar_3_1230.79530.66521
Pξγ0.907haar_3_1230.90540.87415
PQg0.786haar_4_12340.78610.70619
P1aer0.651haar_4_12340.65110.55520
Pγδ0.813haar_4_12340.81310.72113
Pγδπ0.776haar_4_12340.77610.65016
Table 9. The highest value of F1-Score among the rock units on image spectra in study area B using high-frequency features.
Table 9. The highest value of F1-Score among the rock units on image spectra in study area B using high-frequency features.
Rock UnitsHighest F1-ScoreFeatures
(Highest F1-Score)
F1-Score of Haar_4_1234Haar_4_1234 RankingOriginal F1-ScoreOriginal F1-Score Ranking
(E3-N1)t0.875haar_4_12340.87510.79621
Qp2W0.866haar_4_12340.86610.78011
Qp3X0.881haar_3_1230.86330.78017
0.886db2_4_2340.87650.78821
Cξγ0.907db8_4_12340.89660.86421
C2wt0.933haar_4_12340.93310.87221
Table 12. Pearson correlation coefficient between the mean intra-class SAM and OA.
Table 12. Pearson correlation coefficient between the mean intra-class SAM and OA.
Research TargetPearson Correlation Coefficient
Laboratory spectra−0.734
Image spectra (study area A)−0.588
Image spectra (study area B)−0.751
Table 13. Pearson correlation coefficient between the mean intra-class SAM and F1-Score.
Table 13. Pearson correlation coefficient between the mean intra-class SAM and F1-Score.
Laboratory SpectraStudy Area AStudy Area B
RockPearson Correlation CoefficientRock UnitPearson Correlation CoefficientRock UnitPearson Correlation Coefficient
Dolostone−0.971C1x−0.849(E3-N1)t−0.600
Andesite−0.202C1gd0.197Qp2W−0.183
Tuff−0.857C1y−0.610Qp3X−0.583
Limestone−0.818Pξγ−0.638−0.973
Sandstone0.101PQg−0.600Cξγ−0.827
Granite−0.535P1aer−0.700C2wt−0.889
Pγδ−0.163
Pγδπ−0.611
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, S.; Jiang, Q. Improving Rock Classification with 1D Discrete Wavelet Transform Based on Laboratory Reflectance Spectra and Gaofen-5 Hyperspectral Data. Remote Sens. 2023, 15, 5334. https://doi.org/10.3390/rs15225334

AMA Style

Guo S, Jiang Q. Improving Rock Classification with 1D Discrete Wavelet Transform Based on Laboratory Reflectance Spectra and Gaofen-5 Hyperspectral Data. Remote Sensing. 2023; 15(22):5334. https://doi.org/10.3390/rs15225334

Chicago/Turabian Style

Guo, Senmiao, and Qigang Jiang. 2023. "Improving Rock Classification with 1D Discrete Wavelet Transform Based on Laboratory Reflectance Spectra and Gaofen-5 Hyperspectral Data" Remote Sensing 15, no. 22: 5334. https://doi.org/10.3390/rs15225334

APA Style

Guo, S., & Jiang, Q. (2023). Improving Rock Classification with 1D Discrete Wavelet Transform Based on Laboratory Reflectance Spectra and Gaofen-5 Hyperspectral Data. Remote Sensing, 15(22), 5334. https://doi.org/10.3390/rs15225334

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop