1. Introduction
Worldwide, tea plants are produced as a significant revenue cro
[1,2]. Because of its distinct flavor and scent, tea is the most consumed nonalcoholic beverage in the world
[3]. The organic matter in the soil has a large influence on how the tea tree grows. The whole soil ecology depends heavily on soil organic matter, which also helps with nitrogen fixation, carbon sequestration, and plant nutrient retention
[4]. As a result, accurate information on the amount of organic matter in the soil is crucial for enhancing tea plants' ability to thrive. However, conventional techniques for determining soil organic matter are labor intensive, costly, time consuming, and leave behind dangerous chemical residues. The quick detection of soil organic materials during the growth of tea trees cannot be met by this
[5].
Machine vision and hyperspectral imaging are two nondestructive testing techniques that are widely utilized in many industries, including forestry
[6], agriculture
[7], animal husbandry
[8], and fisheries
[9]. To explore the impact of moisture variations on soil urease activity, Chenbo Yang et al.
[10] developed a prediction model for soil urease activity based on hyperspectral data, of which the PLSR model was the most successful (R
2V of 0.8564). Rahimi-Ajdadi et al.
[11] employed soil photographs to estimate the soil water content using feature-finding in several color spaces, and their findings had an average absolute error of less than 1.1%. However, data from a single source may not be perfect in defining some qualities, and individuals have begun to employ fused data in the field of NDT to further enhance target identification. Ting An et al. [
12] collected visual (color pictures) and olfactory (sensor array spectra) data on tea leaves during the fermentation of black tea and combined them with hyperspectral data to create a discriminative model for the level of fermentation that had a 95% prediction set accuracy. Yingqian Yin et al. [
13] combined black tea image features with hyperspectral features to characterize different grades of black tea. They then created discriminant models using partial least squares discriminant analysis (PLS-DA), support vector machines (SVM), and probabilistic neural networks (PNN), with the SVM model having the highest validation set accuracy (98.33%). On the other hand, few studies have been performed on the identification of soil organic matter utilizing fused data. The efficacy of the hyperspectral and machine vision data-fused quantitative soil organic matter prediction model is still being assessed.
Based on the aforementioned research state, the soil of typical Chinese tea plantations was chosen as the research object in this work, and the soil's spectrum information and color and texture attributes were retrieved using the obtained hyperspectral pictures. To combine the spectral and color texture features of the soil, middle-level data fusion is used. From there, a quantitative prediction model of soil organic matter is established to find the best algorithmic strategy, which establishes a theoretical foundation for future advancements in the nondestructive detection technology of soil organic matter.
2. Materials and Methods
2.1. Collection of Soil Samples and Estimation of the Organic Matter Content
The soil samples utilized in this study were gathered in the Shandong Province locations of Jufeng Town, Rizhao City; Lanshan District, Qingdao City, and Junan County, Linyi City. In all of the municipal tea gardens, 15 separate soil collection locations with a combined sampling area of 25 m2 and a depth of 0~40 cm were chosen. Soil samples were then taken ten times from each area for a final total of 150 soil samples.
The gathered samples were taken back to the lab at the Shandong Academy of Agricultural Sciences in Jinan, Shandong Province, China, for the analysis of the organic matter content. To ensure that each soil sample weighed the same weight, contaminants such as stones and animal and plant remains were cleaned before drying. Following that, the criteria for assessing soil organic matter content (NY/T 1121.6-2006) were used to evaluate the organic matter content of all soil samples. Weigh each soil sample that passes through a screen with an aperture of 0.25 mm at 0.05~0.5 g, oxidize the soil's organic carbon with an excessive potassium dichromate sulfuric acid solution, and then titrate the excess potassium dichromate with a ferrous sulfate standard solution. The quantity of organic carbon was determined by converting the amount of potassium dichromate consumed into organic carbon using the oxidation correction coefficient and then multiplying the result by the constant 1.724 to obtain the amount of organic matter in the soil. The organic matter content of each sampling area is shown in
Figure 1.
2.2. Hyperspectral Image Acquisition
ISpecHyper-VS1000-Lab hyperspectral imaging equipment from Lyson Optics was used to gather hyperspectral data. The imaging system consists of an integrated hyperspectral dark box, a hyperspectral camera, a high-performance homogeneous light source, a diffuse reflectance standard version (3%/50%) and a linear displacement sample stage. The wavelength acquisition range is 400~1000 nm with a spectral resolution of 2.5 nm. Before sampling, the machine was warmed for 30 minutes to establish a stable atmosphere. After preheating was completed, the camera lens was covered with a lens cap for dark current collection. For consistent tiling, 50 g of each soil sample was weighed and placed in a test dish. The images were placed within the dark box of the hyperspectral image capture system. Before each sample, the system obtains a calibration picture of the reference plate within the dark box. The soil samples were eventually subjected to 150 hyperspectral pictures.
2.3. Extraction of Spectral Features and Picture Features
Using the aforementioned collected hyperspectral pictures, the spectral information of the soil samples was retrieved. On the hyperspectral picture of each soil sample, 10 ROI sections were chosen at random. The spectral data inside the regions are extracted, and the average of the ten sets of data is used as the soil sample's spectral data. The ENVI 5.3 program was used to carry out this procedure. The collected hyperspectral images were adjusted with the following formula to produce precise reflectance profiles of the samples:
where Icor is the hyperspectral image that has been rectified, Iraw is the picture of the sample that was acquired, Idark is the image of the dark current that was acquired, and Iwhite is the image of the diffuse reflective plate.
The image processing program created by the MATLAB GUI module was used to extract the soil's image characteristics (software copyright number: 2014SR149549). Each hyperspectral image contained 10 randomly chosen zones of interest for the soil, from which the color and texture attributes were retrieved. The Dong et al. [
14] algorithm was then used to convert the RGB color system to the HSV color space. To produce 9 color features and 5 texture features for each image, the data from the 10 regions of interest of each image were averaged. The red channel mean (R), green channel mean (G), blue channel mean (B), color point mean (H), saturation mean (S), luminance mean (V), supergreen transform (2G-R-B), ratio of red channel mean to green channel mean (R/G), and color (hab*) are among the color features; the mean gray value (m), standard deviation (
), smoothness (r), consistency (U), and entropy (e) are among the texture features.
2.4. Preprocessing of the Spectral Data and Characteristic Band Screening
The system may be impacted by external conditions during the measurement process. To improve the spectrum features and remove superfluous information such as baseline drift and noise generated during the measurement process, the spectral data are preprocessed. Three preprocessing methods are used, including standard normalized variate (SNV), multisource scattering correction (MSC) and smoothing.
The wavelength range of the hyperspectral imaging equipment employed in this investigation is 400~1000 nm. The wavelength range of the gathered soil spectra is 344~986 nm, with a resolution of less than 2.5 nm. There are 300 wavelength points in all, which is a significant quantity of duplicated data. As a result, distinctive bands were searched for in the gathered spectral data. In this article, three feature band screening techniques, random frog (RF) [
15], variable combination population analysis (VCPA) [
16] and variable combination population analysis and iterative retained information variable algorithm (VCPA-IRIV) [
17], are chosen to filter out irrelevant data and further enhance the model accuracy.
In modeling, data overfitting can occur often, although feature band screening does not always completely prevent it. As a result, further feature extraction from the data is needed using feature extractions. After feature wavelength screening, the principal component analysis (PCA) dimensionality reduction approach was utilized in this work to further reduce the dimensionality of the data. Using mapping, PCA rearranged the data while keeping the primary spectral characteristics and condensing the many feature variables into a handful of key features [
18,
19].
2.5. Data Fusion
To improve the connectivity between several datasets and provide better modeling outcomes, data fusion is used. In this work, the soil's spectral and color texture data are extractedfrom the obtained hyperspectral pictures of the soil, and the two separate data are combined using a data fusion technique to produce superior prediction outcomes. Low-level fusion, middle-level fusion, and high-level fusion are the three different types of data fusion methodologies. Among these, low-level fusion merely splices the two data matrices and does not yield superior outcomes. If the high-level fusion is not handled carefully, some crucial information may be lost [
20]. As a result, the modeling dataset in this work is processed using the intermediate-level fusion approach. The resultant color texture feature matrix is spliced with the spectral data after it has undergone feature wavelength screening and PCA dimensionality reduction, and the newly formed matrix is then normalized to obtain a new data matrix. The newly acquired matrix will be applied to other modeling tasks. The broad idea of data fusion is shown in
Figure 2.
2.6. Model Building and Evaluation Criteria
The dataset was split into training and validation sets using the KS technique at a ratio of 3:1 after acquiring the fusion data of soil spectral characteristics and color texture features. The dataset consists of 150 sets of fusion data, including 112 sets of fusion data in the training set and 38 sets in the validation set. A quantitative prediction model for soil organic matter using linear PLSR and nonlinear SVR models was created on the basis of this information. To see how fused data affected the forecast of soil organic matter content, the model of single source data and the model of fused data were compared. The correlation coefficient (RC), correlation root mean square error (RMSEC), prediction set correlation coefficient (RP), prediction root mean square error (RMSEP), and relative percentage of deviation (RPD) were employed in this work as assessment indices for the correction and prediction sets [
21]. The more closely the RC and RP approach 1, the more accurate the model is; similarly, the more closely the RMESC and RMESP converge, and the more closely they diverge, the more generalizable the model is. Typically, it is believed that if RPD is less than 1.4, the model is subpar and cannot be used for prediction research, that if RPD is between 1.4 and 1.8, the model can achieve its prediction objective, but accuracy still needs to be improved, and that if RPD is larger than 2, the model has good prediction performance [
22].
3. Results
3.1. Spectral Preprocessing Results
The soil spectral data were collected at 300 wavelength locations between 344 and 986 nm. The spectral data for the wavelength range of 344~402 nm were discarded once it was discovered via observation that the spectral data in this region had more significant noise. Three preprocessing methods, MSC, SNV and smooth, were used to preprocess the spectral data to compare the modeling effect under different preprocessing methods. The spectrum data were preprocessed using three preprocessing techniques, MSC, SNV, and smooth, to examine the modeling impact under various preprocessing techniques. The spectral profiles following the three preprocesses are shown in
Figure 3, together with the soil sample's original spectral profile.
3.2. Characteristic Wavelength Screening Results
The RF, VCPA, and VCPA-IRIV algorithms were used for the screening of the unique wavelengths based on various preprocessing of the spectral data. The stochastic selectivity of the variables and the effectiveness of the RF algorithm in selecting spectral feature variables result in a competitive selection mechanism for the variables that in turn ensures their validity [
23]. The VCPA technique makes use of the exponentially declining function (EDF), model population analysis (MPA), and binary matrix sampling (BMS), which may effectively account for potential interactions between random variables [
24]. A novel approach to variable selection that enables the selective omission of unimportant variables is the VCPA-IRIV method [
25]. The number of iterations in the VCPA-IRIV algorithm is set at 50 for the EDF function and 1000 for the BMS function.
The number of wavelengths screened by the three methods after MSC pretreatment was 10, 11, and 43; the number of wavelengths screened by the three methods after SNV pretreatment was 10, 11, and 45; and the number of wavelengths screened by the three methods after smooth pretreatment was 10, 10, and 28. The results of the typical wavelength screening for various preprocessing steps are displayed in
Figure 4a
–c. It is clear that the predominant concentration of the distinctive wavelengths is in the range of 410 nm, 621~832 nm, and 940 nm. The average spectra of the 15 sampling regions are shown in
Figure 4d, and it is clear that there is an obvious absorption peak at 940 nm. This peak is caused by the OH stretching vibration, which creates a slight water vapor absorption band and can be used to characterize the amount of soil organic matter. The primary matching bands of soil organic matter were found to be between 620~810 nm in the investigation by, and there was only a weak link between soil organic matter and the wavelength at 440 nm [
26]. The typical wavelength distribution from this study is consistent.
3.3. Predictive Modeling of Single Spectral Data
Following feature wavelength extraction from the spectral data and PCA downscaling, PLSR and SVR prediction models for soil organic matter were created. Ten sets of main components were utilized as inputs to the models, and
Table 1 displays the results for each model.
The preprocessing techniques MSC and SNV are generally superior to the SMOOTH procedure, according to a comparison of all models. In a side-by-side comparison of MSC and SNV, the SVR model under the MSC approach comes out on top overall, while some of the PLSR models under the SNV method perform better than MSC. Additionally, all SVR models outperformed PLSR models, showing how effective nonlinear spectral data modeling was in predicting soil organic matter. The best nonlinear model combination overall had R
2P, RSMEP, and RPD values of 0.973, 0.693, and 6.119, respectively, and was MSC+VCPA-IRIA+SVR. The best linear model combination had R
2P, RSMEP, and RPD values of 0.953, 0.895, and 4.711, respectively, and was SNV+VCPA+PLSR.
Figure 5 depicts the connection between these two models' training and prediction sets. The findings show that it is possible to anticipate the amount of soil organic matter using hyperspectral technology.
3.4. Predictive Modeling of Fusion Data
SVR and PLSR models were created based on the fused data of spectral and image characteristics, respectively, to investigate the predictive modeling impact of the fused data.
Table 2 displays the accuracy of each model's predictions for soil organic matter. It is evident from the data in the table that the performance of the models with fused data is always better than the performance of the models with single spectral data. With an improvement in R
2C from 0.901 to 0.980, R
2P from 0.914 to 0.962, and RPD from 3.456 to 4.563, SMOOTH+VCPA-IRIV+SVR is the combination of modeling approaches that has the greatest change in impact. In
Figure 6b,d, the training set-prediction set connections using spectral data and fused data are compared.
Figure 6a,c show how spectral data and fused data are optimized for parameters using the SMOOTH+VCPA-IRIV+SVR combination.
Comparing the modeling performance of the combined data reveals that the best linear model combination is SNV+VCPA-PLSR, with R
2C, R
2P, and RPD values of 0.954, 0.965, and 5.448, respectively, while the best nonlinear model combination is MSC+VCPA-IRIV+SVR, with R
2C, R
2P, and RPD values of 0.995, 0.986, and 8.155, respectively. This is in line with the modeling findings from the spectral data, showing that these two combinations can provide the most accurate predictions of soil organic matter and are thus likely to be used. The model performance results for both combinations are displayed in
Figure 7a,b.
4. Conclusions
In this work, the SVR and PLSR models were built using the intermediate layer fusion approach to combine the soil hyperspectral and picture texture characteristic data. The findings demonstrate that the performance of the fused data models is adequate and that both have outperformed the single spectral data models in the same manner. Out of the three, the nonlinear model performs better at estimating the amount of soil organic matter. The most accurate SVR model is the one built using the MSC and VCPA-IRIV algorithms. It had an R2C of 0.995, an R2P of 0.986, and an RPD of 8.155. The results of this study demonstrate the feasibility of soil organic matter content prediction by fusion data, which improves the accuracy of the prediction model and leads to better modeling results compared to single data. It is beneficial to promote accurate control of soil organic matter content in tea plantations. Targeted irrigation and sowing follow from this.
Supplementary Materials
The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Figure S1: title; Table S1: title; Video S1: title.
Author Contributions
Conceptualization, C.D.; methodology, H.Z. and Z.L.; software, C.D.; validation, C.Y. and X.Z.; formal analysis, H.Z.; investigation, H.Z.; resources, C.D. and X.Z.; data curation, H.Z.; writing—original draft preparation, H.Z; writing—review and editing, C.D.; visualization, H.Z.; supervision, M.L.; project administration, C.D.; funding acquisition, X.Z. All authors have read and agreed to the published version of the manuscript.
Data Availability Statement
The data presented in this study are available on request from the corresponding author.
Acknowledgments
This study was supported by The Key R&D Projects in Zhejiang Province (2022C02010,2023C02043), the Research start-up funds-TRI-SAAS (CXGC2023F18, CXGC2023G33) and the Key Projects of Science and Technology Cooperation in Jiangxi Province (20212BDH80025).
Conflicts of Interest
No conflict of interest exits in the submission of this manuscript, and manuscript is approved by all authors for publication. All the authors listed have approved the manuscript that is enclosed.
References
- Chen, J., Wu, S., Dong, F., Li, J., Zeng, L., Tang, J., and Gu, D. (2021). Mechanism Underlying the Shading-Induced Chlorophyll Accumulation in Tea Leaves. Frontiers in Plant Science 12. [CrossRef]
- Dong, C., Ye, Y., Yang, C., An, T., Jiang, Y., Ye, Y., Li, Y., and Yang, Y. (2021). Rapid detection of catechins during black tea fermentation based on electrical properties and chemometrics. Food Bioscience 40. [CrossRef]
- Rahimi-Ajdadi, F., Abbaspour-Gilandeh, Y., Mollazade, K., and Hasanzadeh, R. P. R. (2018). Development of a novel machine vision procedure for rapid and non-contact measurement of soil moisture content. Measurement 121, 179-189. [CrossRef]
- Hoffland, E., Kuyper, T. W., Comans, R. N. J., and Creamer, R. E. (2020). Eco-functionality of organic matter in soils. Plant and Soil 455, 1-22. [CrossRef]
- İnik, O., İnik, Ö., Öztaş, T., Demir, Y., and Yüksel, A. (2023). Prediction of Soil Organic Matter with Deep Learning. Arabian Journal for Science and Engineering 48, 10227-10247. [CrossRef]
- Zwiazek, J. J., Kyaw, T. Y., Siegert, C. M., Dash, P., Poudel, K. P., Pitts, J. J., and Renninger, H. J. (2022). Using hyperspectral leaf reflectance to estimate photosynthetic capacity and nitrogen content across eastern cottonwood and hybrid poplar taxa. Plos One 17. [CrossRef]
- Yu, F., Bai, J., Jin, Z., Zhang, H., Yang, J., and Xu, T. (2023). Estimating the rice nitrogen nutrition index based on hyperspectral transform technology. Frontiers in Plant Science 14. [CrossRef]
- Chen, D., Wu, P., Wang, K., Wang, S., Ji, X., Shen, Q., Yu, Y., Qiu, X., Xu, X., Liu, Y., and Tang, G. (2022). Combining computer vision score and conventional meat quality traits to estimate the intramuscular fat content using machine learning in pigs. Meat Science 185. [CrossRef]
- Chen, F., Xu, J., Wei, Y., and Sun, J. (2019). Establishing an eyeball-weight relationship for Litopenaeus vannamei using machine vision technology. Aquacultural Engineering 87. [CrossRef]
- Yang, C., Feng, M., Song, L., Jing, B., Xie, Y., Wang, C., Qin, M., Yang, W., Xiao, L., Sun, J., Zhang, M., Song, X., and Kubar, M. S. (2022). Hyperspectral monitoring of soil urease activity under different water regulation. Plant and Soil 477, 779-792. [CrossRef]
- Rahimi-Ajdadi, F., Abbaspour-Gilandeh, Y., Mollazade, K., and Hasanzadeh, R. P. R. (2018). Development of a novel machine vision procedure for rapid and non-contact measurement of soil moisture content. Measurement 121, 179-189. [CrossRef]
- An, T., Huang, W., Tian, X., Fan, S., Duan, D., Dong, C., Zhao, C., and Li, G. (2022). Hyperspectral imaging technology coupled with human sensory information to evaluate the fermentation degree of black tea. Sensors and Actuators B: Chemical 366. [CrossRef]
- Yin, Y., Li, J., Ling, C., Zhang, S., Liu, C., Sun, X., and Wu, J. (2023). Fusing spectral and image information for characterization of black tea grade based on hyperspectral technology. Lwt 185. [CrossRef]
- Dong, C., Liang, G., An, T., Wang, J., Zhu, H. (2018). Near-infrared spectroscopy detection model for sensory quality and chemical constituents of black tea, Transactions of the Chinese Society of, Agricultural Engineering 34, 306–313.
- Sun, J., Yang, W., Feng, M., Liu, Q., and Kubar, M. S. (2020). An efficient variable selection method based on random frog for the multivariate calibration of NIR spectra. RSC Advances 10, 16245-16253. [CrossRef]
- Jiang, H., Xu, W., Ding, Y., and Chen, Q. (2020). Quantitative analysis of yeast fermentation process using Raman spectroscopy: Comparison of CARS and VCPA for variable selection. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy 228. [CrossRef]
- Han, M., Wang, X., Xu, Y., Cui, Y., Wang, L., Lv, D., and Cui, L. (2021). Variable selection for the determination of the soluble solid content of potatoes with surface impurities in the visible/near-infrared range. Biosystems Engineering 209, 170-179. [CrossRef]
- Cebi, N., Yilmaz, M. T., and Sagdic, O. (2017). A rapid ATR-FTIR spectroscopic method for detection of sibutramine adulteration in tea and coffee based on hierarchical cluster and principal component analyses. Food Chemistry 229, 517-526. [CrossRef]
- Cruz-Tirado, J. P., Oliveira, M., de Jesus Filho, M., Godoy, H. T., Amigo, J. M., and Barbin, D. F. (2021). Shelf life estimation and kinetic degradation modeling of chia seeds (Salvia hispanica) using principal component analysis based on NIR-hyperspectral imaging. Food Control 123. [CrossRef]
- Liu, Z., Zhang, R., Yang, C., Hu, B., Luo, X., Li, Y., Dong, C. (2022). Research on moisture content detection method during green tea processing based on machine vision and near-infrared spectroscopy technology,Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy 272, 1386-1425. [CrossRef]
- Ren, G., Wang, Y., Ning, J., and Zhang, Z. (2020). Highly identification of keemun black tea rank based on cognitive spectroscopy: Near infrared spectroscopy combined with feature variable selection. Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy 230. [CrossRef]
- Dong, C., Ye, Y., Yang, C., An, T., Jiang, Y., Ye, Y., Li, Y., and Yang, Y. (2021). Rapid detection of catechins during black tea fermentation based on electrical properties and chemometrics. Food Bioscience 40. [CrossRef]
- Ouyang, Q., Wang, L., Zareef, M., Chen, Q., Guo, Z., and Li, H. (2020). A feasibility of nondestructive rapid detection of total volatile basic nitrogen content in frozen pork based on portable near-infrared spectroscopy. Microchemical Journal 157. [CrossRef]
- Jiang, H., He, Y., Xu, W., and Chen, Q. (2021). Quantitative Detection of Acid Value During Edible Oil Storage by Raman Spectroscopy: Comparison of the Optimization Effects of BOSS and VCPA Algorithms on the Characteristic Raman Spectra of Edible Oils. Food Analytical Methods 14, 1826-1835. [CrossRef]
- Yang, C., Zhao, Y., An, T., Liu, Z., Jiang, Y., Li, Y., and Dong, C. (2021). Quantitative prediction and visualization of key physical and chemical components in black tea fermentation using hyperspectral imaging. Lwt 141. [CrossRef]
- Liu, H., Zhang, Y., and Zhang, B. (2008). Novel hyperspectral reflectance models for estimating black-soil organic matter in Northeast China. Environmental Monitoring and Assessment 154, 147-154. [CrossRef]
|
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. |
© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).