Next Article in Journal
A Collaborative Virtual Walkthrough of Matera’s Sassi Using Photogrammetric Reconstruction and Hand Gesture Navigation
Next Article in Special Issue
Classification of Cocoa Beans by Analyzing Spectral Measurements Using Machine Learning and Genetic Algorithm
Previous Article in Journal
Initial Steps towards a Multilevel Functional Principal Components Analysis Model of Dynamical Shape Changes
Previous Article in Special Issue
Spectral Super-Resolution for High Dynamic Range Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hyperspectral Characteristic Band Selection and Estimation Content of Soil Petroleum Hydrocarbon Based on GARF-PLSR

College of Geo-Exploration Science and Technology, Jilin University, Changchun 130026, China
*
Author to whom correspondence should be addressed.
J. Imaging 2023, 9(4), 87; https://doi.org/10.3390/jimaging9040087
Submission received: 20 February 2023 / Revised: 14 April 2023 / Accepted: 18 April 2023 / Published: 20 April 2023
(This article belongs to the Special Issue Multi-Spectral and Color Imaging: Theory and Application)

Abstract

:
With continuous improvements in oil production, the environmental problems caused by oil exploitation are becoming increasingly serious. Rapid and accurate estimation of soil petroleum hydrocarbon content is of great significance to the investigation and restoration of environments in oil-producing areas. In this study, the content of petroleum hydrocarbon and the hyperspectral data of soil samples collected from an oil-producing area were measured. For the hyperspectral data, spectral transforms, including continuum removal (CR), first- and second-order differential (CR-FD, CR-SD), and Napierian logarithm (CR-LN), were applied to eliminate background noise. At present, there are some shortcomings in the method of feature band selection, such as large quantity, time of calculation, and unclear importance of each feature band obtained. Meanwhile, redundant bands easily exist in the feature set, which seriously affects the accuracy of the inversion algorithm. In order to solve the above problems, a new method (GARF) for hyperspectral characteristic band selection was proposed. It combined the advantage that the grouping search algorithm can effectively reduce the calculation time with the advantage that the point-by-point search algorithm can determine the importance of each band, which provided a clearer direction for further spectroscopic research. The 17 selected bands were used as the input data of partial least squares regression (PLSR) and K-nearest neighbor (KNN) algorithms to estimate soil petroleum hydrocarbon content, and the leave-one-out method was used for cross-validation. The root mean squared error (RMSE) and coefficient of determination (R2) of the estimation result were 3.52 and 0.90, which implemented a high accuracy with only 8.37% of the entire bands. The results showed that compared with the traditional characteristic band selection methods, GARF can effectively reduce the redundant bands and screen out the optimal characteristic bands in the hyperspectral data of soil petroleum hydrocarbon with the method of importance assessment, which retained the physical meaning. It provided a new idea for the research of other substances in soil.

1. Introduction

Since the 20th century, with the rapid growth of the global economy, oil consumption has also increased. In oil exploitation, transportation, storage, production, and processing, accidental leakage and the unreasonable discharge of oily wastewater have caused serious pollution to soil [1]. Petroleum hydrocarbon is a type of persistent organic pollutant, and its low reactivity and antidegradation pose a serious threat to the ecological environment and human health [2,3]. To assess the potential risk of soil pollution and monitor oil pollution, researchers try to measure the content of petroleum hydrocarbons in soil. The common methods are gas chromatography, mass spectrometry, liquid chromatography, and solid-phase microextraction [4]. However, the above methods are complex, costly, time-consuming, and unsuitable for large-area detection. Therefore, a model for the rapid and accurate estimation of soil petroleum hydrocarbon content should be established.
With continuous advancements in remote sensing technology, hyperspectral imaging is developed on the basis of imaging spectroscopy. Hyperspectral imaging can record the reflectance information of targets through thousands of closely arranged wavelengths [5,6]. Compared with traditional optical remote sensing (single-band and multi-band) with a low spectral resolution, it has more imaging channels, continuous spectral range, and more flexible spectral selectivity, which greatly improve the detection capability of targets [7]. Hyperspectral imaging can also reduce the phenomena of different objects having the same spectrum and the same objects having different spectra and detect substances that cannot be detected using traditional remote sensing technology.
Valuable achievements have been attained using hyperspectral technology, such as the estimation of soil heavy metal content, and such technology has been applied to other related fields [8,9]. Considerable research has shown that hyperspectral imaging is closely related to the moisture [10,11,12,13], nutrient [14,15,16,17,18,19], and heavy metal contents [20,21,22,23,24,25,26] and other indicators [27,28] of soil. The application of hyperspectral imaging in the oil pollution detection of soil has also achieved accomplishments. For instance, Foudan et al. [29] utilized hyperspectral images to detect oil spill areas and showed that hyperspectral imaging can eliminate the limitations of traditional methods to the greatest extent. The distribution of oil spill areas by using hyperspectral imaging is also more accurate than that by using traditional methods. Horig et al. [30] used Hymap data to identify oil-contaminated soil and found that they are effective for detecting oil pollution in soil via visible-infrared spectroscopy. A spectral band of 1730 nm can be used as the key characteristic band for petroleum hydrocarbon detection. Kuhn et al. [31] proposed the concept of a hydrocarbon index for hyperspectral imaging and considered that the larger the hydrocarbon index is, the higher the degree of hydrocarbon accumulation will be. Fan et al. [32] measured the hyperspectral data and petroleum hydrocarbon content of samples and analyzed the relationship between soil spectral characteristics and petroleum hydrocarbon content by using a univariate prediction model and stepwise regression method. However, there are some shortcomings in the method of feature band selection, such as large quantity, time of calculation, and unclear importance of each feature band obtained. Meanwhile, redundant bands easily exist in the feature set, which seriously affects the accuracy of the inversion algorithm.
Feature engineering is the key to hyperspectral applications. Compared with the traditional hyperspectral feature extraction methods (PCA, LDA, ICA, etc.) [33], this article proposes a new characteristic band selection method for soil petroleum hydrocarbon hyperspectral imaging based on Genetic Algorithm and Random Forest (GARF). GARF regards the importance of hyperspectral bands as the evaluation index and selects a subset of hyperspectral images as characteristic bands, which can contain most of the important information with a small number of bands. It can effectively reduce the dimension of hyperspectral bands and make the characteristic subset have a clear physical meaning, which provides a basis for further study of the physicochemical mechanism in the hyperspectral imaging of soil petroleum hydrocarbon. The characteristic bands selected via GARF were used as input data of partial least squares regression (PLSR) and k-nearest neighbor (KNN) to estimate the content of petroleum hydrocarbons in soil samples. The experimental results showed that this method obtains an excellent estimation result, which provides a novel method for large-scale accurate, rapid, and low-cost detection of soil petroleum hydrocarbon. We conducted all experiments on a computer with Intel® Core™ (Santa Clara, CA, USA) i7-8700K CPU at 3.70 GHz, 16 GB running memory, and all of our algorithms were coded using python and the open-source modules Scikit-learn.

2. Materials and Methods

2.1. Soil Sample Collection and Spectral Data Acquisition

The soil sample collection area was in Daqing City, Heilongjiang Province, which is located in the southwest of Daqing Oilfield. The sampling area has a temperate continental climate, with four distinct seasons, annual precipitation about 600 mm, and abundant water resources. The types of soil in the sampling area are chernozem and meadow soil; the pH ranges from 5.4 to 10.7, with an average of 9.0, belonging to weak alkaline soil [34]. When collecting soil samples, we took each sampling point as the center with sampling area 1 m2. Five random locations of soil were stripped of impurities and collected at each sampling grid, and then were put into the sampling bag and mixed fully, with a sampling depth within 15 cm. In consideration of objective factors, such as oil well layout, oil leakage, and traffic, 28 soil samples were collected in the study area (Figure 1). The petroleum hydrocarbon was measured using the gas chromatography method (BS-EN-ISO 16703:2011). Three of the soil samples had abnormal values, in which the concentration of petroleum hydrocarbon was considerably higher than other samples. The analysis showed that the abnormal samples were caused by mixing black oily sludge. Oily sludge is a kind of sludge mixed with heavy oil, such as asphalt and various refined and residue oil, which is not inherent in nature. Therefore, the three abnormal samples were removed, and the spectral data of the remaining 25 soil samples were measured (Table 1). An ASD FieldSpec3 spectro-radiometer (Analytical Spectral Devices Inc., Boulder, CO, USA) was applied to measure the spectra of soil samples, and the spectral range was determined to be 350–2500 nm. A 50 W halogen lamp was used as the light source to simulate sunlight, and the zenith angle was 30°. The field of view angle and the measurement distance of the probe were 15° and 5 cm, respectively. A whiteboard was used for reflection calibration, and measurements of the soil spectrum were carried out in a dark laboratory. Each sample was measured 15 times to ensure the accuracy of the spectral data. The average spectral value was calculated as the reflectance value of the sample. In the ranges of 350–379 nm and 2401–2500 nm, the signal-to-noise ratio was low. Thus, the bands before 380 nm and after 2400 nm were discarded. The sampling interval of the spectrometer was 1 nm; hence, 2021 bands were obtained in a range of 380–2400 nm. Owing to the high spectral resolution, information overlap existed between adjacent bands, which made the results highly vulnerable to noise. Therefore, the spectral data were resampled, and the interval was 10 nm.

2.2. Spectral Data Preprocessing

On the basis of denoising and resampling, the initial spectrum was processed through continuum removal (CR), which can highlight the absorption, reflection, and emission characteristics of the spectrum [35]. First- and second-order differential (CR-FD, CR-SD) and the Napierian logarithm (CR-LN) were used to eliminate the noise in the background of the CR spectrum [36]. The initial spectrum of soil samples is shown in Figure 2a. The initial reflectance of the soil samples was between 0 and 0.25, and the fluctuation shape of the spectral curves of each sample was similar. In the visible bands, with the increase in wavelengths, the reflectance gradually enlarged and tended to be stable at 1200 nm. In the near-infrared bands, the reflectance of the samples fluctuated greatly. Two distinct absorption valleys were distributed at approximately 1400 and 1900 nm, and a slightly sunken one existed at 2200 nm. Figure 2b–e show the results of four different transformations of the initial spectrum. All the four transformations can amplify the initial spectrum. After transformations, the reflectance fluctuated more remarkably at approximately 1400, 1900, and 2200 nm.

2.3. Model Principle

Feature selection is a technology for dimension reduction, which is the process of selecting a feature subset from a feature set. In the hyperspectral field, feature selection is also called band selection [37]. Compared with feature selection, feature extraction may obtain slightly better results in most cases, but the features converted through feature extraction have poor interpretability. The inherent physical meaning of hyperspectral imaging is also lost. However, feature selection directly selects a feature subset from the original feature set, such that it retains the physical meaning as opposed to feature extraction [38]. In view of the above reason, band selection is more widely used in hyperspectral dimensionality reduction than feature extraction. For hyperspectral imaging, dozens of optimal characteristic bands are difficult to select from hundreds or thousands of spectral bands on account of hundreds of millions of band combinations. At present, two strategies, namely group search and point-by-point search, are applied to select a feature subset from a complex feature set. The former method is to generate several candidate feature subsets continuously, reserve relatively good subsets, feed back to the subset generation strategy, evaluate and select new subsets, and repeat the above process until the candidate feature subsets meet the set requirements. The optimal feature subset can be considered as the result of feature selection. The latter method is to start from an empty set and select one (or more) feature from the feature set to join the subset until the subset contains the required number of features or to start from the set including all features and delete the features that meet the conditions one by one (or more) until the number of remaining features meets the set requirements.

2.3.1. Genetic Algorithm

Genetic algorithm (GA) is a typical group search algorithm, which is based on natural selection and genetic theory. It combines the survival of the fittest rules in the process of biological evolution with the random information exchange mechanism of chromosomes in a population [39,40]. Furthermore, GA is also an efficient global optimization search algorithm, which abandons the traditional search strategy and simulates the biological evolution process in nature to search a feature subset randomly. It regards a possible solution to a problem as an individual or a chromosome of a population and codes each individual into a symbol string to simulate the evolution process of Darwinian genetic selection and natural elimination. In accordance with the evolutionary rules of survival of the fittest, each individual is evaluated using a predetermined objective fitness function, and the better population is retrieved continuously. At the same time, the global parallel search method is used to search the optimal individuals in the optimization population to reserve an optimal solution. The implementation process of GA is as follows:
  • Coding: The transformation of a feasible solution of a practical problem from its solution space to the search space that can be processed using GA is called coding. The most common coding method is binary coding.
  • Population analysis and design: GA randomly generates a certain number of individuals, from which better individuals are selected to form the initial population. In the iterative process, the larger the population size is, the higher the chance to obtain an optimal solution, and the smaller the possibility of the algorithm falling into a local minimum. However, the large population size will lead to an increase in the time consumption of the algorithm.
  • Fitness function: A fitness function is applied to evaluate the optimization process of individuals in the population and estimate the degree close to the optimal solution.
  • Crossover: GA imitates the process of gene recombination into new chromosomes in nature. Some genes in chromosomes are exchanged between two pairs of chromosomes, and a crossover operator is used to form two new individuals.
  • Mutation: Mutation is introduced to induce the formation of new individuals and increase the ability to find the optimal solution.
  • Termination of calculation: The individual with the maximum fitness value reserved in the evolution process is selected as the output of the optimal solution.

2.3.2. Random Forest

Random Forest (RF) is an integrated algorithm composed of numerous decision trees. The idea of RF is to build an excellent tree, which needs to select excellent features. Therefore, the importance of each feature must be judged. RF randomly changes the value of a feature and then compares the error rate of the test set before and after the change. It adopts a point-by-point search strategy. The difference in error rate is regarded as the importance of the feature, and the mean decrease in accuracy (MDA) is considered the index to evaluate the feature importance of RF [41].
M D A = 1 n t = 1 n e r r O O B t e r r O O B t
where n is the number of based learners, and e r r O O B t is the out-of-pocket error after noise is added. The more the M D A index decreases, the greater the effect of the corresponding characteristic on the estimation results, and the higher its importance will be.

2.3.3. Partial Least Squares Regression

Partial least squares regression (PLSR) is a multiple linear regression algorithm, which is one of the most used regression algorithms in hyperspectral imaging [42]. When hyperspectral data are used to estimate soil petroleum hydrocarbon content, the independent variable (X) in the estimation model is soil spectral data, and the dependent variable (Y) is the petroleum hydrocarbon content. The process of PLSR is summarized as follows:
X = T P T + E
Y = U Q T + F
where X is the predictive matrix, Y is the response matrix, and T and U are the projection matrices of X and Y , respectively. P and Q are the orthogonal load matrices, and matrices E and F are the error terms.

2.3.4. K-Nearest Neighbor

K-nearest neighbor (KNN) is a supervised machine learning algorithm. The basic idea of KNN is to traverse the training set, find the k training samples closest to the new sample according to the distance formula, and use the majority voting principle to determine the prediction value of the new sample. It is widely used for dealing with classification and regression problems [43].

2.3.5. Performance Evaluation Scales

To evaluate the performance of models in this study, root mean squared error (RMSE) and coefficient of determination ( R 2 ) were used as indicators to assess the accuracy and stability of the models:
R 2 = 1 i = 1 n y i Y i 2 i = 1 n y i y ¯ 2
R M S E = 1 n i = 1 n y i Y i 2
where n is the number of samples, y i is the measured value, Y i is the predicted value, and y ¯ is the average of the measured values. The lower the RMSE is, the closer the R 2 is to 1, and the higher the accuracy and stability of the estimation model are.

3. Results and Discussion

We used GA to screen the hyperspectral imaging and select the characteristic bands of soil petroleum hydrocarbon. However, the combination of characteristic bands found using GA was still complex. To further eliminate redundant features, RF was used to sort the importance of characteristic bands and conducted secondary screening to select the optimal characteristic bands of soil petroleum hydrocarbon, whose importance was greater than the average. The optimal characteristic bands were regarded as the input data of the PLSR and KNN to estimate the soil petroleum hydrocarbon content. In order to facilitate the determination of parameters, we used the Grid Search method in this experiment. The Grid Search method only needs to input the parameter range, and it can automatically adjust parameters and output the optimal results, avoiding the tedious work of manual adjustment. It is a very suitable automatic parameter adjustment method for small data sets. The data set was divided into two parts by using the leave-one-out method. The advantage of this method is that every iteration used the maximum number of samples as the training set and made the estimation results clearly reflect the accuracy of the model (Figure 3).

3.1. Selection of Optimal Characteristic Bands

GA was used to screen the initial and transformation spectra in accordance with the genetic mechanism and natural selection. The population size, iterations, crossover rate, and mutation rate of GA were 50, 100, 0.6, and 0.02, respectively. Figure 4 shows the distribution of bands marked using GA for five types of spectra. On the basis of the operation of selection, crossover, and mutation with continuous genetic iteration, the bands with better fitness function values were reserved (bands were marked as 1), whereas the bands with worse fitness function values were eliminated (bands were marked as 0). The numbers of characteristic bands for initial, CR, CR-FD, CR-SD, and CR-LN spectra were 108, 91, 98, 108, and 91, respectively (Table 2). Characteristic bands of five spectral forms are listed in Appendix A.
Table 2 shows that the number of bands selected using GA was still large, and data redundancy remained among spectral bands. The band combination was screened for a second time through the function of out-of-bag estimation of RF, and the bands with greater-than-average importance were selected as the optimal characteristic bands. This process not only greatly reduced the dimension of the bands but also retained the most important characteristic information, which provided an effective data basis for the subsequent estimation of soil petroleum hydrocarbon content. Table 3 shows the number and importance of optimal characteristic bands after secondary screening by using RF. Figure 5 depicts the distribution of optimal characteristic bands and the importance value.

3.2. Estimation Accuracies of Soil Petroleum Hydrocarbon Content

PLSR is a nonparametric regression analysis method based on factor analysis, which is highly suitable for modeling under the condition of a small number of high-dimensional samples. Many studies have shown that it has excellent performance in the spectral estimation of soil material content [44,45,46]. Table 3 indicates that the optimal characteristic bands of soil petroleum hydrocarbon hyperspectral imaging selected using GARF can extract most of the important information in only a small amount of bands: CR-GARF had the best effect, in which 82% of the important information was contained in 17 characteristic bands (420, 1220, 1230, 1720, 1760, 1780, 1790, 1830, 2190, 2210, 2260, 2300, 2310, 2340, 2350, 2360, 2390 nm, only 8.37% of entire bands). Partial characteristic bands selected by GARF were similar to the conclusions reached by predecessors: Cloutis et al. studied the reflection characteristics of petroleum hydrocarbon in the visible near-infrared bands and summarized that there were two absorption bands near 1730 nm and 2310 nm [47]; Gao et al. perceived that in the near-infrared bands, crude oil solid had relatively obvious absorption bands near 1700 nm and 2300 nm [48]; nine absorption peaks of oil-contaminated soil were identified in the 1725–14,000 nm range through indoor spectral measurements by Zhu et al., with three of them (1725, 2310, 2348 nm) located in the visible near-infrared bands [49]; Feng et al. screened the above characteristic bands and believed that the double absorption peaks near 1748 nm and 2330 nm were the key wavebands for detecting the petroleum hydrocarbon in soil [50]; based on the sampling analysis of Gudong Oilfield, Fan et al. determined that the 1690–1790 nm range was the main area for estimating the content of petroleum hydrocarbon in soil [32]; Wang et al. conducted spectral analysis of soil samples with different oil contents using the visible near-infrared bands and found that the double absorption near 1748 nm and 2330 nm can be used as characteristic bands for soil oil pollution research [51]. The characteristic bands obtained in the above study obviously coincided with the characteristic bands selected in this article within band ranges of 1720–1790 nm and 2300–2350 nm, which proved the rationality of GARF.
Then, we input the initial, CR, and CR-GARF spectra into PLSR and KNN for comparative analysis to validate the performance in estimation of soil petroleum hydrocarbon content. Table 4 and Figure 6 show that CR-GARF-PLSR can accurately estimate the content of soil petroleum hydrocarbon with fewer bands ( R M S E = 3.52, R 2 = 0.90), which indicated that GARF can validly reduce redundant bands and screen out the optimal characteristic bands of soil petroleum hydrocarbon. Compared with the estimation accuracy of Initial-PLSR ( R M S E = 6.83, R 2 = 0.62), that of CR-PLSR ( R M S E = 5.50, R 2 = 0.75) was improved, which demonstrated that CR can highlight the characteristic information of the initial spectrum, remove background noise, and improve the estimation result. From Figure 6c,d, PLSR was more accurate than KNN in the estimation of soil petroleum hydrocarbon content. For validation of GARF in the band selection of the CR spectrum, the characteristic bands of the CR spectrum selected using GARF calculated the correlation coefficient among bands and between the petroleum hydrocarbon content of soil samples. Figure 7 demonstrated that the optimal characteristic bands selected using GARF had a high correlation with soil petroleum hydrocarbon content (|correlation coefficient| ≥ 0.6) [52]. The correlation coefficient among bands was insignificant except for adjacent intervals (Figure 8).

4. Conclusions

In this article, two tasks were implemented in order to estimate soil petroleum hydrocarbon content: a. laboratory analysis on the collected samples with an ASD spectro-radiometer; b. hyperspectral estimation based on laboratory analysis. For rapid and accurate estimation of soil petroleum hydrocarbon content, CR-GARF-PLSR was proposed. In the experiment, the R M S E and R 2 of the model were 3.52 and 0.90, which illustrated that it was an effective method in estimating soil petroleum hydrocarbon content. The experimental conclusions were as follows: 1. CR can eliminate the background noise in spectral data, highlight the absorption and reflection characteristics of spectral curves, and contribute to a more accurate estimation result. 2. GARF can effectively remove the redundant bands in soil petroleum hydrocarbon hyperspectral imaging and retain the optimal characteristic bands, which is a new method for feature selection based on machine learning. 3. Compared with other models, CR-GARF-PLSR had better performance in estimating soil petroleum hydrocarbon content, which provided a new idea for the research of other substances in soil. However, there are still some difficulties in the detection of oil pollution in soil by using hyperspectral technology. For instance, the differences in the composition of oil and the environment in different oil fields and the water content in crude oil will definitely have an impact on the measured spectra. The corresponding characteristic bands will also change accordingly, so the scalability and applicability of the model need further verification. At the same time, the method proposed in this paper still needs indoor testing, which will increase the cost of the experiment. In the near-infrared band, absorbance is also a common spectral measurement parameter, which may be helpful in improving experimental accuracy.
In this study, we established a new model for hyperspectral characteristic band selection and estimation of soil petroleum hydrocarbon content, which achieved an excellent performance. Furthermore, the research in this paper can be improved in the future. Due to the limitation of policy, we only collected a small number of samples. If the number of samples is enough, we can use more advanced regression methods, such as XGBoost and convolutional neural network, to make the final result more accurate. Similarly, we cannot remove samples containing other contaminations, except for black oily sludge, because of the limited number of samples. If we can fulfill this assumption, our experimental accuracy may be improved. With the resolution of the mixed pixel and the noise generation in field applications, the combination of this research and airborne hyperspectral technology has great potential in large-scale accurate and low-cost rapid detection of soil petroleum hydrocarbon.

Author Contributions

P.S. developed the original idea for the study and completed the manuscript. Q.J. revised the paper. Z.L. contributed to data acquisition and preprocessing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the geological survey project of the China Geological Survey (DD20190705).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to express our gratitude to the reviewers for their insightful comments and discussion, especially regarding the suggestions for spectral measurement.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Characteristic bands of Initial selected by GA.
Table A1. Characteristic bands of Initial selected by GA.
Position of Characteristic Bands (nm)
400420430450460490500510550560580590
600610620630640680690710720740750760
77081082084085089095097098099010301040
105010701100111011201130114011501210122012301260
127012801290130013101320133013501380139014001410
146014801490150015201540155015601580162016301650
166016701730174017601770181018401860188018901910
194020202040206020702080215021602170218021902210
222022402250226022802290230023102350236023902400
Table A2. Characteristic bands of CR selected by GA.
Table A2. Characteristic bands of CR selected by GA.
Position of Characteristic Bands (nm)
390410420430440460490520540560580590
620640650670680690730790800830880890
9009209609709801010104010501080114011701180
122012301240126012901330136013901400146014801500
151015401550156015701610163016601720176017801790
183018501890191019201930194020002010203020402060
207020802090210021102120213021402150219022102250
2260230023102340235023602390
Table A3. Characteristic bands of CR-FD selected by GA.
Table A3. Characteristic bands of CR-FD selected by GA.
Position of Characteristic Bands (nm)
380390410450460500510520530540550560
580620630670710730740760780790810820
84087088090091093094010001030105010601080
109011001120114011501200122012401340135013601380
139014101440146014801490150015201530154016101630
164016501660167017101720176017801790183018401850
188019001960197019902000203020602070208021102120
216021902230224022502260227022802290230023202340
23702400
Table A4. Characteristic bands of CR-SD selected by GA.
Table A4. Characteristic bands of CR-SD selected by GA.
Position of Characteristic Bands (nm)
400420430440450470480510520540570580
590600610620670680690710740760770780
810820830850860870880890910920930940
9609801010104010801100113011801200124012601280
131013301350138014001410143014401450147014901530
155015601580160016101630164016801690174017501760
177017801800181018201840186018701880189019101940
195019701980200020202040205020702080209021702190
220022102220224022602290230023202330235023902400
Table A5. Characteristic bands of CR-LN selected by GA.
Table A5. Characteristic bands of CR-LN selected by GA.
Position of Characteristic Bands (nm)
390410420430440450490500510520540560
580600610620630640680690730750760770
79081083084085090095096097098010101030
104010501060107010801140118012101220123012401250
128013301360139014101440148015001510154015601570
161016301660172017601780183018501890192019301960
197020002030204020502070208021002110212021302140
2150221022502300231023602390

References

  1. Achard, V.; Foucher, P.-Y.; Dubucq, D. Hydrocarbon Pollution Detection and Mapping Based on the Combination of Various Hyperspectral Imaging Processing Tools. Remote Sens. 2021, 13, 1020. [Google Scholar] [CrossRef]
  2. Escandar, G.M.; de la Pena, A.M. Multi-way calibration for the quantification of polycyclic aromatic hydrocarbons in samples of environmental impact. Microchem. J. 2021, 164, 106016. [Google Scholar] [CrossRef]
  3. Truskewycz, A.; Gundry, T.D.; Khudur, L.S.; Kolobaric, A.; Taha, M.; Aburto-Medina, A.; Ball, A.S.; Shahsavari, E. Petroleum Hydrocarbon Contamination in Terrestrial Ecosystems-Fate and Microbial Responses. Molecules 2019, 24, 3400. [Google Scholar] [CrossRef] [PubMed]
  4. Chen, Z.L.; Yin, W.Q.; Liu, H.T.; Liu, Q.; Yang, Y. Review of Monitoring Petroleum-Hydrocarbon Contaminated Soils with Visible and Near-Infrared Spectroscopy. Spectrosc. Spectr. Anal. 2017, 37, 1723–1727. [Google Scholar]
  5. Lin, N.; Liu, H.-Q.; Yang, J.-J.; Wu, M.-H.; Liu, H.-L. Hyperspectral Estimation of Soil Nutrient Content in the Black Soil Region Based on BA-Adaboost. Spectrosc. Spectr. Anal. 2020, 40, 3825–3831. [Google Scholar]
  6. Xie, F.; Lei, C.; Yang, J.; Jin, C. An Effective Classification Scheme for Hyperspectral Image Based on Superpixel and Discontinuity Preserving Relaxation. Remote Sens. 2019, 11, 1149. [Google Scholar] [CrossRef]
  7. Sun, G.; Zhang, A.; Ren, J.; Ma, J.; Wang, P.; Zhang, Y.; Jia, X. Gravitation-Based Edge Detection in Hyperspectral Images. Remote Sens. 2017, 9, 592. [Google Scholar] [CrossRef]
  8. Bai, X.; Xiao, Q.; Zhou, L.; Tang, Y.; He, Y. Detection of Sulfite Dioxide Residue on the Surface of Fresh-Cut Potato Slices Using Near-Infrared Hyperspectral Imaging System and Portable Near-Infrared Spectrometer. Molecules 2020, 25, 1651. [Google Scholar] [CrossRef]
  9. Zhu, S.; Zhang, J.; Chao, M.; Xu, X.; Song, P.; Zhang, J.; Huang, Z. A Rapid and Highly Efficient Method for the Identification of Soybean Seed Varieties: Hyperspectral Images Combined with Transfer Learning. Molecules 2020, 25, 152. [Google Scholar] [CrossRef]
  10. Wu, L.-G.; Wang, S.-L.; He, J.-G. Study on Soil Moisture Mechanism and Establishment of Model Based on Hyperspectral Imaging Technique. Spectrosc. Spectr. Anal. 2018, 38, 2563–2570. [Google Scholar]
  11. Kano, Y.; Mcclure, W.F.; Skaggs, R.W. A Near Infrared Reflectance Soil Moisture Meter. Trans. ASAE—Am. Soc. Agric. Eng. 1985, 28, 1852–1855. [Google Scholar] [CrossRef]
  12. Bowman, G.E.; Hooper, A.W.; Hartshorn, L. A prototype infrared reflectance moisture meter. J. Agric. Eng. Res. 1985, 31, 67–79. [Google Scholar] [CrossRef]
  13. Whalley, W.R.; Leeds-Harrison, P.B.; Bowman, G.E. Estimation of soil moisture status using near infrared reflectance. Hydrol. Process. 1991, 5, 321–327. [Google Scholar] [CrossRef]
  14. He, Y.; Huang, M.; Garcia, A.; Hernandez, A.; Song, H. Prediction of soil macronutrients content using near-infrared spectroscopy. Comput. Electron. Agric. 2007, 58, 144–153. [Google Scholar] [CrossRef]
  15. Daniel, K.W.; Tripathi, X.K.; Honda, K. Artificial neural network analysis of laboratory and in situ spectra for the estimation of macronutrients in soils of Lop Buri (Thailand). Aust. J. Soil Res. 2003, 41, 47–59. [Google Scholar] [CrossRef]
  16. Mouazen, A.M.; Kuang, B.; Baerdemaeker, J.D.; Ramon, H. Comparison among principal component, partial least squares and back propagation neural network analyses for accuracy of measurement of selected soil properties with visible and near infrared spectroscopy. Geoderma 2010, 158, 23–31. [Google Scholar] [CrossRef]
  17. Yang, H.; Kuang, B.; Mouazen, A.M. Quantitative analysis of soil nitrogen and carbon at a farm scale using visible and near infrared spectroscopy coupled with wavelength reduction. Eur. J. Soil Sci. 2012, 63, 410–420. [Google Scholar] [CrossRef]
  18. Luan, F.-M.; Zhang, X.-L.; Xiong, H.-G.; Zhang, F.; Wang, F. Comparative Analysis of Soil Organic Matter Content Based on Different Hyperspectral Inversion Models. Spectrosc. Spectr. Anal. 2013, 33, 196–200. [Google Scholar]
  19. Ma, Y.; Jiang, Q.-G.; Meng, Z.-G.; Liu, H.-X. Black Soil Organic Matter Content Estimation Using Hybrid Selection Method Based on RF and GABPSO. Spectrosc. Spectr. Anal. 2018, 38, 181–187. [Google Scholar]
  20. Wang, Y.; Ma, H.; Wang, J.; Liu, L.; Pietikainen, M.; Zhang, Z.; Chen, X. Hyperspectral monitor of soil chromium contaminant based on deep learning network model in the Eastern Junggar coalfield. Spectrochim. Acta Part A—Mol. Biomol. Spectrosc. 2021, 257, 119739. [Google Scholar] [CrossRef]
  21. Zeng, R.; Zhao, Y.G.; Li, D.C.; Wu, D.W.; Wei, C.L.; Zhang, G.L. Selection of “Local” Models for Prediction of Soil Organic Matter Using a Regional Soil Vis-NIR Spectral Library. Soil Sci. 2016, 181, 13–19. [Google Scholar] [CrossRef]
  22. Riedel, F.; Denk, M.; Mueller, I.; Barth, N.; Glaesser, C. Prediction of soil parameters using the spectral range between 350 and 15,000 nm: A case study based on the Permanent Soil Monitoring Program in Saxony, Germany. Geoderma 2018, 315, 188–198. [Google Scholar] [CrossRef]
  23. Kooistra, L.; Wehrens, R.; Leuven, R.; Buydens, L.M.C. Possibilities of visible-near-infrared spectroscopy for the assessment of soil contamination in river floodplains. Anal. Chim. Acta 2001, 446, 97–105. [Google Scholar] [CrossRef]
  24. Lian, S.; Ji, J.; De-Jun, T.; Hong-Bing, X.; Zhen-Fu, L.; Bo, G. Estimate of heavy metals in soil and streams using combined geochemistry and field spectroscopy in Wan-sheng mining area, Chongqing, China. Int. J. Appl. Earth Obs. Geoinf. 2015, 34, 9. [Google Scholar] [CrossRef]
  25. Choe, E.; van der Meer, F.; van Ruitenbeek, F.; van der Werff, H.; de Smeth, B.; Kim, Y.-W. Mapping of heavy metal pollution in stream sediments using combined geochemistry, field spectroscopy, and hyperspectral remote sensing: A case study of the Rodalquilar mining area, SE Spain. Remote Sens. Environ. 2008, 112, 3222–3233. [Google Scholar] [CrossRef]
  26. Wang, J.-F.; Wang, S.-J.; Bai, X.-Y.; Liu, F.; Lu, Q.; Tian, S.-Q.; Wang, M.-M. Prediction Soil Heavy Metal Zinc Based on Spectral Reflectance in Karst Area. Spectrosc. Spectr. Anal. 2019, 39, 3873–3879. [Google Scholar]
  27. Waiser, T.H.; Morgan, C.L.S.; Brown, D.J.; Hallmark, C.T. In situ characterization of soil clay content with visible near-infrared diffuse reflectance spectroscopy. Soil Sci. Soc. Am. J. 2007, 71, 389–396. [Google Scholar] [CrossRef]
  28. Bilgili, A.V.; Cullu, M.A.; Van Es, H.; Aydemir, A.; Aydemir, S. The Use of Hyperspectral Visible and Near Infrared Reflectance Spectroscopy for the Characterization of Salt-Affected Soils in the Harran Plain, Turkey. Arid. Land Res. Manag. 2011, 25, 19–37. [Google Scholar] [CrossRef]
  29. Salem, F.; Kafatos, M.; El-Ghazawi, T.; Gomez, R.; Yang, R.X. Hyperspectral image assessment of oil-contaminated wetland. Int. J. Remote Sens. 2005, 26, 811–821. [Google Scholar] [CrossRef]
  30. Horig, B.; Kuhn, F.; Oschutz, F.; Lehmann, F. HyMap hyperspectral remote sensing to detect hydrocarbons. Int. J. Remote Sens. 2001, 22, 1413–1422. [Google Scholar] [CrossRef]
  31. Kuhn, F.; Oppermann, K.; Horig, B. Hydrocarbon Index—An algorithm for hyperspectral detection of hydrocarbons. Int. J. Remote Sens. 2004, 25, 2467–2473. [Google Scholar] [CrossRef]
  32. Fan, Y.; Zhang, L. Soil oil content hyperspectral model in Gudong Oilfield. J. Remote Sens. 2012, 16, 378–389. [Google Scholar]
  33. Kumar, B.; Dikshit, O.; Gupta, A.; Singh, M.K. Feature extraction for hyperspectral image classification: A review. International Journal of Remote Sensing 2020, 41, 6248–6287. [Google Scholar] [CrossRef]
  34. Zhang, H.; Zheng, Z.Z.; Yang, H. Discrimination of Heavy Metal Sources in Topsoil in Zhaoyuan County Based on Multivariate Statistics and Geostatistical. Soil 2017, 49, 819–827. [Google Scholar] [CrossRef]
  35. Li, Y.; Yin, Y.; Yu, H.; Yuan, Y. Fast detection of water loss and hardness for cucumber using hyperspectral imaging technology. J. Food Meas. Charact. 2021, 16, 76–84. [Google Scholar] [CrossRef]
  36. Goodin, D.G.; Han, L.; Fraser, R.N.; Rundquist, D.C.; Stebbins, W.A. Analysis of suspended solids in water using remotely sensed high resolution derivative spectra. Photogramm. Eng. Remote Sens. 1993, 59, 505–510. [Google Scholar] [CrossRef]
  37. Feng, J.; Jiao, L.C.; Liu, F.; Sun, T.; Zhang, X.R. Unsupervised feature selection based on maximum information and minimum redundancy for hyperspectral images. Pattern Recognit. 2016, 51, 295–309. [Google Scholar] [CrossRef]
  38. Zhang, W.; Li, X.; Dou, Y.; Zhao, L. A Geometry-Based Band Selection Approach for Hyperspectral Image Analysis. IEEE Trans. Geosci. Remote Sens. 2018, 56, 4318–4333. [Google Scholar] [CrossRef]
  39. Yang, Z.; Xiao, H.; Zhang, L.; Feng, D.; Zhang, F.; Jiang, M.; Sui, Q.; Jia, L. Fast determination of oxides content in cement raw meal using NIR-spectroscopy and backward interval PLS with genetic algorithm. Spectrochim. Acta Part A—Mol. Biomol. Spectrosc. 2019, 223, 117327. [Google Scholar] [CrossRef]
  40. Aghelpour, P.; Mohammadi, B.; Biazar, S.M.; Kisi, O.; Sourmirinezhad, Z. A Theoretical Approach for Forecasting Different Types of Drought Simultaneously, Using Entropy Theory and Machine-Learning Methods. ISPRS Int. J. Geo-Inf. 2020, 9, 701. [Google Scholar] [CrossRef]
  41. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  42. Zhang, W.; Cao, A.; Shi, P.; Cai, L. Rapid evaluation of freshness of largemouth bass under different thawing methods using hyperspectral imaging. Food Control 2021, 125, 108023. [Google Scholar] [CrossRef]
  43. Pang, L.; Wang, J.; Men, S.; Yan, L.; Xiao, J. Hyperspectral imaging coupled with multivariate methods for seed vitality estimation and forecast for Quercus variabilis. Spectrochim. Acta Part A—Mol. Biomol. Spectrosc. 2021, 245, 118888. [Google Scholar] [CrossRef]
  44. Zhang, Y.; Hartemink, A.E.; Huang, J.; Townsend, P.A. Synergistic use of hyperspectral imagery, Sentinel-1 and LiDAR improves mapping of soil physical and geochemical properties at the farm-scale. Eur. J. Soil Sci. 2021, 72, 1690–1717. [Google Scholar] [CrossRef]
  45. Crusiol, L.G.T.; Nanni, M.R.; Furlanetto, R.H.; Sibaldelli, R.N.R.; Cezar, E.; Sun, L.; Foloni, J.S.S.; Mertz-Henning, L.M.; Nepomuceno, A.L.; Neumaier, N.; et al. Yield Prediction in Soybean Crop Grown under Different Levels of Water Availability Using Reflectance Spectroscopy and Partial Least Squares Regression. Remote Sens. 2021, 13, 977. [Google Scholar] [CrossRef]
  46. Xu, X.; Chen, S.; Xu, Z.; Yu, Y.; Zhang, S.; Dai, R. Exploring Appropriate Preprocessing Techniques for Hyperspectral Soil Organic Matter Content Estimation in Black Soil Area. Remote Sens. 2020, 12, 3765. [Google Scholar] [CrossRef]
  47. Cloutis, E.A. Spectral reflectance properties of hydrocarbons: Remote-sensing implications. Science 1989, 245, 165–168. [Google Scholar] [CrossRef]
  48. Gao, L.; Yang, B. A Study on Near Infrared Spectral Characteristics of Petroleum Matter Applied to Remote Sensing of Oil Gas Resources. Remote Sens. Land Resour. 1991, 4, 9–12+29. [Google Scholar]
  49. Zhu, Z. Hydrocarbon Microseepage Theory and Oil-Gas Reservoir Detecting by Remote Sensing. Remote Sens. Technol. Appl. 1994, 31, 10. [Google Scholar]
  50. Feng, X.; Shi, Y. Near Infrared Spectroscopy and Its Application in the Analysis of Petroleum Products; China Petrochemical Press: Beijing, China, 2022; Volume 1, p. 238. ISBN 7-80164-239-2. [Google Scholar]
  51. Wang, X.; Tian, Q.; Guan, Z. The Extraction of Oiland Gas Information by Using Hyperion Imagery in the Sebei Gas Field. Remote Sens. Nat. Resour. 2006, 71, 36–40+107. [Google Scholar]
  52. Chen, C.; Jiang, Q.; Zhang, Z.; Shi, P.; Xu, Y.; Liu, B.; Xi, J.; Chang, S. Hyperspectral Inversion of Petroleum Hydrocarbon Contents in Soil Based on Continuum Removal and Wavelet Packet Decomposition. Sustainability 2020, 12, 4218. [Google Scholar] [CrossRef]
Figure 1. The position of sampling sites in the study area.
Figure 1. The position of sampling sites in the study area.
Jimaging 09 00087 g001
Figure 2. Five spectral forms of samples: (a) initial reflectance; (b) continuum removal of reflectance; (c) first derivative of reflectance; (d) second derivative of reflectance; (e) napierian logarithm of reflectance.
Figure 2. Five spectral forms of samples: (a) initial reflectance; (b) continuum removal of reflectance; (c) first derivative of reflectance; (d) second derivative of reflectance; (e) napierian logarithm of reflectance.
Jimaging 09 00087 g002
Figure 3. Schematic diagram of leave-one-out method.
Figure 3. Schematic diagram of leave-one-out method.
Jimaging 09 00087 g003
Figure 4. Distribution of bands marked using GA.
Figure 4. Distribution of bands marked using GA.
Jimaging 09 00087 g004
Figure 5. Distribution of optimal characteristic bands after RF selection: (a) Initial-GARF; (b) CR-GARF; (c) CR-FD-GARF; (d) CR-SD-GARF; (e) CR-LN-GARF.
Figure 5. Distribution of optimal characteristic bands after RF selection: (a) Initial-GARF; (b) CR-GARF; (c) CR-FD-GARF; (d) CR-SD-GARF; (e) CR-LN-GARF.
Jimaging 09 00087 g005
Figure 6. Scatter plots of the predicted versus measured soil petroleum hydrocarbon content: (a) Initial-PLSR; (b) CR-PLSR; (c) CR-GARF-PLSR; (d) CR-GARF-KNN.
Figure 6. Scatter plots of the predicted versus measured soil petroleum hydrocarbon content: (a) Initial-PLSR; (b) CR-PLSR; (c) CR-GARF-PLSR; (d) CR-GARF-KNN.
Jimaging 09 00087 g006
Figure 7. Correlation coefficient between soil petroleum hydrocarbon content and CR-GARF optimal characteristic bands.
Figure 7. Correlation coefficient between soil petroleum hydrocarbon content and CR-GARF optimal characteristic bands.
Jimaging 09 00087 g007
Figure 8. Heatmap of the correlation coefficient of CR-GARF optimal characteristic bands (wavelength/nm).
Figure 8. Heatmap of the correlation coefficient of CR-GARF optimal characteristic bands (wavelength/nm).
Jimaging 09 00087 g008
Table 1. Statistical description of the measured petroleum hydrocarbon content of soil samples (Min: minimum; max: maximum; ave: average; SD: standard deviation; CV: coefficient of variation).
Table 1. Statistical description of the measured petroleum hydrocarbon content of soil samples (Min: minimum; max: maximum; ave: average; SD: standard deviation; CV: coefficient of variation).
No.Min
(g/kg−1)
Max
(g/kg−1)
Ave
(g/kg−1)
SD
(g/kg−1)
CV
250.008134.18.4611.02130.23%
Table 2. Statistical table of bands selected using GA.
Table 2. Statistical table of bands selected using GA.
Spectral FormNumber of Resampled BandsNumber of Selected Bands
Initial203108
CR20391
CR-FD20398
CR-SD203108
CR-LN20391
Table 3. Importance of optimal characteristic bands after secondary screening by using RF.
Table 3. Importance of optimal characteristic bands after secondary screening by using RF.
Spectral FormNo.AveSum
Initial-GA470.00930.721
CR-GA170.01020.822
CR-FD-GA330.01100.810
CR-SD-GA340.00930.818
CR-LN-GA300.01100.790
Table 4. Estimation accuracies of soil petroleum hydrocarbon content.
Table 4. Estimation accuracies of soil petroleum hydrocarbon content.
ModelRMSER2
Initial-PLSR6.830.62
CR-PLSR5.500.75
CR-GARF-PLSR3.520.90
CR-GARF-KNN4.700.82
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

Shi, P.; Jiang, Q.; Li, Z. Hyperspectral Characteristic Band Selection and Estimation Content of Soil Petroleum Hydrocarbon Based on GARF-PLSR. J. Imaging 2023, 9, 87. https://doi.org/10.3390/jimaging9040087

AMA Style

Shi P, Jiang Q, Li Z. Hyperspectral Characteristic Band Selection and Estimation Content of Soil Petroleum Hydrocarbon Based on GARF-PLSR. Journal of Imaging. 2023; 9(4):87. https://doi.org/10.3390/jimaging9040087

Chicago/Turabian Style

Shi, Pengfei, Qigang Jiang, and Zhilian Li. 2023. "Hyperspectral Characteristic Band Selection and Estimation Content of Soil Petroleum Hydrocarbon Based on GARF-PLSR" Journal of Imaging 9, no. 4: 87. https://doi.org/10.3390/jimaging9040087

APA Style

Shi, P., Jiang, Q., & Li, Z. (2023). Hyperspectral Characteristic Band Selection and Estimation Content of Soil Petroleum Hydrocarbon Based on GARF-PLSR. Journal of Imaging, 9(4), 87. https://doi.org/10.3390/jimaging9040087

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