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).
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].
Figure 3.
The FWHM of GF-5 AHSI data.
Figure 3.
The FWHM of GF-5 AHSI data.
Figure 4.
The method flowchart of this study.
Figure 4.
The method flowchart of this study.
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).
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.
Figure 7.
Preprocessing process of GF-5 AHSI data.
Figure 7.
Preprocessing process of GF-5 AHSI data.
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.
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.)
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.)
Figure 11.
Mean of intra-class SAM in laboratory spectra.
Figure 11.
Mean of intra-class SAM in laboratory spectra.
Figure 12.
The standard deviation of intra-class SAM in laboratory spectra.
Figure 12.
The standard deviation of intra-class SAM in laboratory spectra.
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.
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.)
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.)
Figure 16.
Mean values of intra-class SAM in study area A.
Figure 16.
Mean values of intra-class SAM in study area A.
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.
Figure 18.
Mean values of intra-class SAM in study area B.
Figure 18.
Mean values of intra-class SAM in study area B.
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.
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.
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.
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.)
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.
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].
Table 1.
Names and photos of rock samples collected in the field.
Table 2.
The parameters of GF-5 AHSI.
Table 2.
The parameters of GF-5 AHSI.
Parameters | GF-5 AHSI |
---|
Orbit altitude | 705 km |
Swath width | 60 km |
Spatial resolution | 30 m |
Spectral resolution | VNIR: 5 nm; SWIR: 10 nm |
Number of bands | VNIR: 150; SWIR: 180 |
Spectral range | 0.39–2.51 μm |
SWIR Signal-to-Noise Ratio (SNR) | ~500 |
Dispersive systems | Grating |
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 Method | Abbreviation (Using Haar Wavelets as an Example) |
---|
Reconstruct the high-frequency coefficients of the second and third levels after the three-level decomposition | haar_3_23 |
Reconstruct all three high-frequency coefficients after three-level decomposition | haar_3_123 |
Reconstruct the high-frequency coefficients of the second, third, and fourth levels after four-level decomposition | haar_4_234 |
Reconstruct all four high-frequency coefficients after four-level decomposition | haar_4_1234 |
Table 4.
Bands removed in GF-5 AHSI data preprocessing.
Table 4.
Bands removed in GF-5 AHSI data preprocessing.
Wavelength Range | Band Number | Wavelength |
---|
VNIR | VNIR 1–4 | 390–403 nm |
SWIR | SWIR 1–4 | 1004–1030 nm |
SWIR 43–50 | 1359–1418 nm |
SWIR 96–112 | 1805–1940 nm |
Table 5.
OA and Kappa coefficients of laboratory spectra.
Table 5.
OA and Kappa coefficients of laboratory spectra.
Features | OA | Kappa Coefficients |
---|
haar_4_1234 | 0.599 | 0.502 |
haar_4_234 | 0.595 | 0.496 |
db2_4_234 | 0.593 | 0.493 |
haar_3_23 | 0.580 | 0.476 |
haar_3_123 | 0.574 | 0.467 |
db2_4_1234 | 0.572 | 0.465 |
db4_4_1234 | 0.556 | 0.441 |
db4_4_234 | 0.556 | 0.443 |
db2_3_23 | 0.551 | 0.435 |
db10_4_234 | 0.535 | 0.413 |
db2_3_123 | 0.530 | 0.405 |
db10_4_1234 | 0.526 | 0.403 |
db8_4_234 | 0.514 | 0.384 |
db8_4_1234 | 0.511 | 0.380 |
db4_3_23 | 0.489 | 0.346 |
db4_3_123 | 0.474 | 0.324 |
original spectrum | 0.465 | 0.335 |
db10_3_23 | 0.443 | 0.277 |
db10_3_123 | 0.432 | 0.262 |
db8_3_123 | 0.430 | 0.260 |
db8_3_23 | 0.428 | 0.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 Types | Highest F1-Score | Features (Highest F1-Score) | F1-Score of Haar_4_1234 | Haar_4_1234 Ranking | Original F1-Score | Original F1-Score Ranking |
---|
Dolostone | 0.689 | db2_4_234 | 0.655 | 7 | 0.464 | 18 |
Andesite | 0.418 | db2_4_234 | 0.383 | 2 | 0.271 | 9 |
Tuff | 0.519 | db2_3_123 | 0.462 | 4 | 0.398 | 18 |
Limestone | 0.665 | Haar_4_1234 | 0.665 | 1 | 0.457 | 21 |
Sandstone | 0.465 | Haar_4_234 | 0.456 | 2 | 0.270 | 8 |
Granite | 0.743 | Haar_4_234 | 0.740 | 2 | 0.643 | 16 |
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 A | Study Area B |
---|
Features | OA | Kappa Coefficients | Features | OA | Kappa Coefficients |
---|
haar_4_1234 | 0.769 | 0.736 | haar_4_1234 | 0.885 | 0.862 |
haar_4_234 | 0.768 | 0.735 | haar_3_123 | 0.884 | 0.861 |
haar_3_123 | 0.766 | 0.733 | haar_4_234 | 0.880 | 0.856 |
haar_3_23 | 0.760 | 0.725 | haar_3_23 | 0.878 | 0.853 |
db4_4_234 | 0.721 | 0.681 | db2_4_234 | 0.870 | 0.844 |
db2_3_123 | 0.720 | 0.68 | db2_4_1234 | 0.868 | 0.842 |
db2_4_1234 | 0.720 | 0.68 | db10_4_1234 | 0.858 | 0.829 |
db2_4_234 | 0.718 | 0.677 | db2_3_123 | 0.857 | 0.828 |
db2_3_23 | 0.717 | 0.676 | db2_3_23 | 0.854 | 0.825 |
db10_4_1234 | 0.709 | 0.668 | db8_4_1234 | 0.853 | 0.823 |
db10_4_234 | 0.709 | 0.667 | db4_3_123 | 0.849 | 0.819 |
db4_4_1234 | 0.709 | 0.667 | db10_4_234 | 0.849 | 0.819 |
db8_4_1234 | 0.700 | 0.658 | db8_4_234 | 0.849 | 0.819 |
db4_3_23 | 0.700 | 0.657 | db4_4_1234 | 0.849 | 0.819 |
db8_4_234 | 0.699 | 0.656 | db4_3_23 | 0.848 | 0.818 |
db4_3_123 | 0.698 | 0.655 | db4_4_234 | 0.848 | 0.817 |
db8_3_123 | 0.686 | 0.641 | db8_3_123 | 0.839 | 0.807 |
db8_3_23 | 0.685 | 0.64 | db10_3_123 | 0.836 | 0.803 |
original | 0.683 | 0.637 | db8_3_23 | 0.833 | 0.799 |
db10_3_123 | 0.679 | 0.633 | db10_3_23 | 0.833 | 0.799 |
db10_3_23 | 0.662 | 0.614 | original | 0.813 | 0.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 Units | Highest F1-Score | Features (Highest F1-Score) | F1-Score of Haar_4_1234 | Haar_4_1234 Ranking | Original F1-Score | Original F1-Score Ranking |
---|
C1x | 0.745 | haar_4_234 | 0.737 | 2 | 0.627 | 21 |
C1gd | 0.720 | haar_3_123 | 0.679 | 3 | 0.643 | 7 |
C1y | 0.802 | haar_3_123 | 0.795 | 3 | 0.665 | 21 |
Pξγ | 0.907 | haar_3_123 | 0.905 | 4 | 0.874 | 15 |
PQg | 0.786 | haar_4_1234 | 0.786 | 1 | 0.706 | 19 |
P1aer | 0.651 | haar_4_1234 | 0.651 | 1 | 0.555 | 20 |
Pγδ | 0.813 | haar_4_1234 | 0.813 | 1 | 0.721 | 13 |
Pγδπ | 0.776 | haar_4_1234 | 0.776 | 1 | 0.650 | 16 |
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 Units | Highest F1-Score | Features (Highest F1-Score) | F1-Score of Haar_4_1234 | Haar_4_1234 Ranking | Original F1-Score | Original F1-Score Ranking |
---|
(E3-N1)t | 0.875 | haar_4_1234 | 0.875 | 1 | 0.796 | 21 |
Qp2W | 0.866 | haar_4_1234 | 0.866 | 1 | 0.780 | 11 |
Qp3X | 0.881 | haar_3_123 | 0.863 | 3 | 0.780 | 17 |
Cγ | 0.886 | db2_4_234 | 0.876 | 5 | 0.788 | 21 |
Cξγ | 0.907 | db8_4_1234 | 0.896 | 6 | 0.864 | 21 |
C2wt | 0.933 | haar_4_1234 | 0.933 | 1 | 0.872 | 21 |
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 Target | Pearson 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 Spectra | Study Area A | Study Area B |
---|
Rock | Pearson Correlation Coefficient | Rock Unit | Pearson Correlation Coefficient | Rock Unit | Pearson Correlation Coefficient |
---|
Dolostone | −0.971 | C1x | −0.849 | (E3-N1)t | −0.600 |
Andesite | −0.202 | C1gd | 0.197 | Qp2W | −0.183 |
Tuff | −0.857 | C1y | −0.610 | Qp3X | −0.583 |
Limestone | −0.818 | Pξγ | −0.638 | Cγ | −0.973 |
Sandstone | 0.101 | PQg | −0.600 | Cξγ | −0.827 |
Granite | −0.535 | P1aer | −0.700 | C2wt | −0.889 |
| | Pγδ | −0.163 | | |
| | Pγδπ | −0.611 | | |