Skip to main content

ORIGINAL RESEARCH article

Front. Plant Sci., 13 July 2022
Sec. Sustainable and Intelligent Phytoprotection
This article is part of the Research Topic Remote Sensing Application for Precision Agriculture View all 22 articles

Estimation Model of Potassium Content in Cotton Leaves Based on Wavelet Decomposition Spectra and Image Combination Features

\nQiushuang YaoQiushuang YaoZe Zhang
Ze Zhang*Xin Lv
Xin Lv*Xiangyu ChenXiangyu ChenLulu MaLulu MaCong SunCong Sun
  • The Key Laboratory of Oasis Eco-Agriculture, College of Agriculture, Shihezi University, Shihezi, China

Potassium (K) is one of the most important elements influencing cotton metabolism, quality, and yield. Due to the characteristics of strong fluidity and fast redistribution of the K in plants, it leads to rapid transformation of the K lack or abundance in plant leaves; therefore, rapid and accurate estimation of potassium content in leaves (LKC, %) is a necessary prerequisite to solve the regulation of plant potassium. In this study, we concentrated on the LKC of cotton in different growth stages, an estimation model based on the combined characteristics of wavelet decomposition spectra and image was proposed, and discussed the potential of different combined features in accurate estimation of the LKC. We collected hyperspectral imaging data of 60 main-stem leaves at the budding, flowering, and boll setting stages of cotton, respectively. The original spectrum (R) is decomposed by continuous wavelet transform (CWT). The competitive adaptive reweighted sampling (CARS) and random frog (RF) algorithms combined with partial least squares regression (PLSR) model were used to determine the optimal decomposition scale and characteristic wavelengths at three growth stages. Based on the best “CWT spectra” model, the grayscale image databases were constructed, and the image features were extracted by using color moment and gray level co-occurrence matrix (GLCM). The results showed that the best decomposition scales of the three growth stages were CWT-1, 3, and 9. The best growth stage for estimating LKC in cotton was the boll setting stage, with the feature combination of “CWT-9 spectra + texture,” and its determination coefficients (R2val) and root mean squared error (RMSEval) values were 0.90 and 0.20. Compared with the single R model (R2val = 0.66, RMSEval = 0.34), the R2val increased by 0.24. Different from our hypothesis, the combined feature based on “CWT spectra + color + texture” cannot significantly improve the estimation accuracy of the model, it means that the performance of the estimation model established with more feature information is not correspondingly better. Moreover, the texture features contributed more to the improvement of model performance than color features did. These results provide a reference for rapid and non-destructive monitoring of the LKC in cotton.

Introduction

Potassium (K) is an essential and favorite nutrient element in the growth of cotton. The level of K directly affects the growth and development of cotton and the quality of fiber (Pettigrew, 2008; Lewis et al., 2021). Affected by the parent material of soil formation, the soil in Xinjiang, China is rich in K. The content of available K in arable soil tends to be high in the north and low in the south, but in recent years, cotton fields in some areas have been deficient in K (Tian et al., 2020; Wang et al., 2021). Furthermore, Xinjiang has high-quality and high-yield cotton and a large demand for soil nutrients. Therefore, an excessive supply of potash fertilizer is usually used to avoid production reduction due to lack of the K, resulting in an increase in cotton production costs. However, the accurate method of element determination is time-consuming and laborious, so it is of great significance to monitor the potassium content in cotton leaves (LKC) quickly and without damage for the healthy growth of cotton, the recommendation of fertilizer application amounts and the reduction in resource waste.

Proximal hyperspectral remote sensing technology has become an effective means to evaluate precision agriculture (Pandey et al., 2017; Li et al., 2019), which can be divided into imaging spectra and non-imaging spectra. They can collect hyperspectral reflectance data from the visible, near-infrared (NIR) and short-wave infrared (SWIR) regions of the electromagnetic spectrum (Mertens et al., 2021), so that a wide variety of physiological traits of crops can be studied, such as crop nutrient deficiency (Furlanetto et al., 2021; Jiang et al., 2021; Mahajan et al., 2021, photosynthetic efficiency (El-Hendawy et al., 2017), water stress (Sun et al., 2021; Zhou et al., 2021), chlorophyll fluorescence (Zhao et al., 2021), heavy metal pollution (Lin et al., 2021) and early plant disease detection (El-Hendawy et al., 2017; Barros et al., 2020). On the other hand, hyperspectral imaging can simultaneously obtain the target spectrum and image information, and is regarded as a technique with high-throughput plant phenotype potential (Pandey et al., 2017). Although there are many studies on nutrition monitoring using near-end hyperspectral imaging, most of them focus on quantitative monitoring and diagnosis of crop nitrogen (N), such as wheat (Mahajan et al., 2014; Jiang et al., 2021), rice (Men et al., 2021), maize (Furlanetto et al., 2021), cotton (Oliveira et al., 2020), rape (Liu et al., 2020a), soybean (Chen et al., 2019), orange (Osco et al., 2019, 2020a), tea (Wang et al., 2020) and mango (Mahajan et al., 2021). At present, the quantitative monitoring research on crop K is also gradually carried out, but more studies often analyze the K together with other elements (Liu et al., 2020b; Osco et al., 2020a,b; Mahajan et al., 2021), and there are few studies only on the characteristics of single the K nutrient element. Indeed, a large group of K+ transporters and channels has been identified in plants (Gierth and Maser, 2007), and cytoplasmic concentration of K+ is maintained around 80–150 MM (Ahmad and Maathuis, 2014). Preserving this concentration range is important for many physiological processes as the enzyme activations, and stabilization of protein synthesis (Villette et al., 2020). These processes are present in all tissues and subcellular compartments of cells, which enables the precise quantification of foliar K attributes of the foliage. It has been shown that the 550–700 and 1,390–1,880-nanometer (nm) wavelengths were the best wavelengths to explain the difference in nutrient levels of N, P, and K in cotton (Oliveira et al., 2020; Wang et al., 2020). Thus, the research utilizing sensitive characteristic wavelengths or vegetation indexes to identify and estimate the K deficiency are common method in rice (Das et al., 2020), wheat (Hussain et al., 2017), and maize (Furlanetto et al., 2021). However, to which extent the K can be estimated using hyperspectral requires further investigations.

Continuous wavelet transform (CWT) has attracted increasing attention in image and spectral signal decomposition due to its rich wavelet basis function, multi-resolution, and time-frequency locality (Chen et al., 2010; Yue et al., 2020). Because the CWT can perform multi-scale decomposition of spectrum and has good performance in characteristic wavelength selection and fine spectral signal extraction (Chen et al., 2019), it has been widely used in crop biochemical parameter inversion of hyperspectral data (Zhang et al., 2014), including estimating the above-ground biomass of wheat (Yao et al., 2018; Yue et al., 2020), analyzing the relationship between leaf copper content and spectrum (Lin et al., 2021) and rapidly detecting the chlorophyll fluorescence parameters of potato leaves (Zhao et al., 2021). Therefore, it is of great significance to improve the accuracy of spectral monitoring to construct a quantitative regression relationship between the wavelet coefficients and nutrient parameters (Mahajan et al., 2014).

As an imaging spectrometer can provide very high spatial and spectral resolution data (Pandey et al., 2017), it is necessary to consider the spatial information (e.g., color and texture) in hyperspectral images in addition to the spectral information to estimate crop nutrients. Image color can express the color distribution and range of image, while image texture reflects the information of uniformity, sharpness and spatial arrangement of image gray distribution. Although there are few studies on the role of image features in hyperspectral nutrient monitoring, it has important application potential in the field of hyperspectral imaging (Jiang et al., 2021). Zheng et al. (2017) extracted 14 vegetation indices related to color features to segment corn, and the accuracy rate over 90.19%. Zou et al. (2019) segmented broccoli seedlings from weeds and soil by extracting GLCM features and color features, and achieved higher accuracy. In the existing studies have demonstrated that the K deficiency causes discoloration of crop leaf tips and edges [such as wheat (Mahajan et al., 2014), rice (Sun et al., 2018), soybean (Ghosal et al., 2018), and cotton (Oliveira et al., 2020)], then gradually spread to the center of the leaf, develop into brown spots, and finally wither and necrosis, resulting in changes in leaf color and texture (Laddi et al., 2013). Also, vegetation coverage and NDVI value are significantly reduced (Severtson et al., 2016). Besides, through the calculation of crop RGB image, it was found that the extension rate of the K deficient leaves slowed down and the wilting rate accelerated (Sun et al., 2018). However, the potential for the image features of leaf hyperspectral imaging data for estimating crop nutrients stress (e.g., K) is not well documented.

Hence, using the high-resolution proximal hyperspectral imaging data of cotton leaves in different growth stages, this study proposed an estimation model of the LKC in cotton based on the combined characteristics of “CWT spectra + image.” The main objectives of this study were to (1) clarify the characteristics of hyperspectral response of cotton LKC at different growth stages, and the effective characteristic wavelengths of the best decomposition scale was determined combined with CWT and PLSR, (2) construct a gray image database of characteristic wavelengths in different growth periods to extract and screen sensitive image features, and (3) evaluate the potential of different “CWT spectra + image” combination features to estimate the LKC of cotton at different growth stages.

Materials and Methods

Experimental Design

The research area was located in Erlian (85°59′41″E, 44°19′54″N), the teaching experiment field of Shihezi University. Sunshine duration is 2,721–2,818 h, ≥ 0°C active accumulated temperature is 4,023–4,118°C, ≥10°C active accumulated temperature is 3,570–3,729°C and frost-free period is 168–171 days. The soil texture was loam, and the 0–20 cm soil layer contained 19.06 g·kg−1 organic matter, 12.8 mg·kg−1 total nitrogen, 20.8 mg·kg−1 available phosphorus, and 165.1 mg·kg−1 available potassium. The soil pH is 8.17 and electrical conductivity (EC) is 0.42 ms·cm−1. During the whole growth period of cotton, nitrogen, phosphorus, and potassium fertilizer were applied with water drops. The urea (N, 46%) of 276 kg·hm−2, monoammonium phosphate (P2O5, 61%) of 174 kg·hm−2 and potassium sulfate (K2O, 50%) was used as a potassium fertilizer. A total of 9-times drips were given during the whole growth period, and the fertilization ratios of the three fertilizers (N, P, and K) were 2.5, 7.7, and 0% (June 7), 7.5, 11.7, and 6.7% (June 15), 7.5, 11.7, and 6.7% (June 24), 12.5, 19.2, and 20% (July 2), 20, 19.2, and 20% (July 18), 25, 15.4, and 13.3% (July 26), 15, 15.4, and 13.3% (Aug 5), 10, 0, and 13.3% (Aug 15), 0, 0, and 6.7% (Aug 25). At each fertilization, the three fertilizers weighed in proportion are poured into the corresponding differential pressure fertilization tank to dissolve, and then drip irrigation was applied to the plot.

The experiment was carried out in the study area from April to September 2020. The variety Xinluzao 53 was selected for the experiment. The planting pattern was “one film, three tubes, and six rows” and the plant spacing was 10 + 66 + 10 cm. Four K application levels were set, namely, blank (0 kg·hm−2), low K (75 kg·hm−2), conventional K (225 kg·hm−2), and high K (375 kg·hm−2). We used a random block design with 3 replicates on a total of 12 plots with a single plot area of 25 m2 (Figure 1). The sowing date was 18 April 18 2020, the topping date was 9 July 202 and the sampling periods were the budding stage (30 June 2020), flowering stage (12 July 2020) and boll setting stage (30 July 2020). Five pieces of cotton main-stem leaves with similar growth in the middle and upper parts were randomly collected from each plot, a total of 60 main-stem leaves were collected in one growth period, and a total of 180 main-stem leaves were collected in three growth periods.

FIGURE 1
www.frontiersin.org

Figure 1. The study site and the location of the experiments.

Hyperspectral Image Data Acquisition

The SOC710-VP portable visible-near-infrared hyperspectral imaging spectrometer (Surface Optics Corporation, USA) was used for data acquisition. The spectral resolution is 5 nm, the image resolution is 692 × 520 and the spectral range is 376–1,044 nm, with a total of 128 bands. After removing the front and rear spectral noises, each hyperspectral image cube selected a wavelength in the range of 400–950 nm, with a total of 106 image bands. To reduce the influence of natural light, all images were captured in a dark box (Figure 2).

FIGURE 2
www.frontiersin.org

Figure 2. Hyperspectral image acquisition system. System consisting of a hyperspectral imager, dark box, lighting system, lifting platform, and a computer.

Cotton main-stem leaves are the main source organs providing assimilates to cotton bolls (Pace et al., 1999). Fresh main-stem leaves at the three key growth stages were selected to be tiled in a dark box with a low-reflectivity black background plate according to the order of leaf position. At the same time, a standard gray plate was placed 5 cm away from the leaf edge to assist black-and-white correction. To reduce the influence of light source intensity, exposure time, and dark current in the sensor during spectral scanning, the instrument should be preheated for 30 min. The scanning parameters of the hyperspectral imager were as follows: object distance, 88 cm, scanning rate, 150–200 frames·s−1, aperture, 5.6. The collected spectral data were digital (DN) and were converted into spectral reflectance through spectral calibration and radiometric calibration in SRAnal 710 software according to the grayscale reference panel in each original image. The average spectral reflectance of the whole leaf was extracted as the original spectral data of this sample (Lin et al., 2021).

Determination of Total Potassium Content in Plants

The total potassium content in leaves was determined using a H2SO4-H2O2 flame photometer (Bao, 2000). Fresh leaves were dried at 85°C for 30 min and then at 105°C until reaching a constant weight. The dried leaf samples were ground, weighed and then digested with H2O2-H2SO4, and the K was determined using a laboratory flame photometer (FP640, Yidian Co., Ltd, Shanghai, China). The total LKC was calculated according to the following formula:

K(%)=ρ×V×ts×10-4m    (1)

where ρ is the mass concentration of K obtained from the standard curve (ug·mL−1), V is the measuring liquid volume (ml), ts is the separation multiple, and m is the dry sample mass (g).

Data Processing

Continuous Wavelet Transform

Continuous wavelet transform (CWT) is an effective signal processing tool to decompose an original signal into multidimensional signals, mainly including discrete wavelet transform and continuous wavelet transform (Liu et al., 2020b). In the CWT, the algorithm uses the selected mother wavelet to decompose the hyperspectral data into a series of wavelet coefficients of different scales, which is a linear transformation. Its transformation formula is as follows:

Wf(a,b)=-f(λ)ψa,bdλ    (2)
ψa,b(λ)=1aψ(λ-ba)    (3)

where f (λ) is the leaf hyperspectral reflectance, λ is the wavelength within 400–950 nm, ψa,b is the wavelet basis function, a is the scale factor, b is the translocation factor, and Wf(a, b) is a two-dimensional matrix, including i and j, where i represents the decomposition scale (i = 1, 2, 3, …, m) and j represents the band range of the spectrum (j = 1, 2, 3, …, n), forming an m by n matrix. In this study, the spectral data of cotton leaves at three growth stages were obtained, and each leaf sample included 106 bands. As the setting of the decomposition scale has a certain influence on spectral feature recognition (Liu et al., 2020; Lin et al., 2021), the decomposition scale in this study was set as 21, 22, 23, …, 210, scales 1–10. Among them, the decomposition scales 1–3 and 4–7 belong to low frequency and middle frequency, respectively, and the rest belong to high frequency. Sym2 was selected as the wavelet basis function, and then, a PLSR model was used to quantitatively analyze the relationship between the wavelet coefficients of each decomposition scale and the LKC so as to determine the optimal decomposition scale and effective wavelength.

Selection Method of Characteristic Variables

Different feature selection methods lead to different features being selected. To select spectral and image features sensitive and stable to the LKC in cotton, the competitive adaptive reweighted sampling (CARS), and random frog (RF) algorithms were used in this study to screen features.

Competitive adaptive reweighted sampling selects wavelength points with a large coefficient absolute value in the model through the Monte Carlo strategy and removes wavelength points with a low weight (Sun et al., 2021). The subset with the lowest root mean squared error of cross validation (RMSECV) value is retained as the feature selection result by cross-validation. In this study, the Monte Carlo strategy was set to run 50 times, using 5-fold cross-validation.

Random frog is an algorithm to measure the importance of variables (El-Hendawy et al., 2017). The main steps are as follows: (1) a subset of initial variables containing m variables is randomly initialized, (2) variables in the initial variable subset are continuously selected into the candidate subset, and the number of variables in the candidate subset increases and decreases with the number of iterations, (3) the selection probability of each variable is calculated as a measure of the importance of the variable, and (4) the characteristic wavelength is selected according to the probability of the occurrence of recorded variables in each iteration. In this study, the selection probability of each wavelength was used to screen the feature information, and the running results are presented in descending order. The number of iterations was set as 10,000, and the selection probability thresholds of the three growth periods were 0.40, 0.21, and 0.23, respectively.

Image Feature Extraction

The most common gray level co-occurrence matrix (GLCM) algorithm was adopted to extract texture features (Yang et al., 2021). In this study, the energy (ENE), entropy (ENT), contrast (CON), correlation (COR), and their mean (MEA) and variance (VAR) in four directions were calculated by using the gray comatrix function. The calculation equation is shown in Table 1, where the P(i, j) is the value of the GLCM in the ith row and jth column, k is the number of gray levels in the GLCM. The gray level is 256, the step size is 1, the angle is 0°, 45°, 90°, and 135°. Finally, each characteristic wavelength grayscale image will eventually produce 24 (4 × 4 + 8) texture features.

TABLE 1
www.frontiersin.org

Table 1. The calculation equations for the characteristics of the GLCM.

Color moments are used to represent the color distribution in the image (Ge et al., 2021). Since the color information is mainly distributed in low-order moments, first-order moments (mean, MEA), second-order moments (variance, VAR), and third-order moments (skewness, SKE) are sufficient to express the color distribution of the image. Its formula is as follows:

MEA=j=1n1nPij    (4)
VAR=1nj=1n(PijMEAi)2    (5)
SKE=(1nj=1nP(PijMEAi)33    (6)

where Pij is the color value of the jth pixel on the ith color channel, i is the number of color channels of the image. The image in this study is grayscale image, so i = 1; MEAi is the color mean of the ith color channel of all pixels. Finally, each characteristic wavelength grayscale image will eventually produce three color features.

Notably, the named representation of combined features is as follows: (1) The characteristic wavelength-texture feature-direction, such as 400 nm-ENE-0°, which means the texture feature is ENE in the 0° direction of the 400-nm grayscale image and (2) characteristic wavelength-color feature, such as 400-nm MEA, which means the color feature is MEA of the 400-nm grayscale image.

Modeling and Analysis Methods

Partial least squares regression (PLSR) is one of the most widely used modeling methods in spectral analysis, which can be used for dimensionality reduction and comprehensive screening of spectral data, with high modeling stability and reliability. The PLSR is widely favored in hyperspectral analysis (Lin et al., 2021; Zhao et al., 2021) because it can solve the collinearity and overfitting characteristics of hyperspectral data compared with other multivariate models.

The determination coefficients (R2) and RMSE values were used to evaluate the performance of the model. In general, better performing models have higher R2 and lower RMSE values. Original hyperspectral data were extracted by ENVI5 3. The CWT, PLSR, and GLCM analyses of leaf spectral data were carried out by Matlab R2018a (The MathWorks, Inc., Natick MA, USA). Origin 2020 was used for creating graphs (OriginLab Corporation, Northampton, MA, USA).

R2=1-i=1n(yi-ŷi)2i=1n(yi-y¯)2    (7)
RMSE=1ni=1n(yi-ŷi)2    (8)

where R2cal is expressed as the determination coefficient of calibration sets, R2val is expressed as the determination coefficient of validation sets, n is the number of samples, yi and ŷi, are, respectively, the measured and estimated values of sample i in the corresponding sample set, ȳ is the average value of yi.

Results

Analysis of Spectral Characteristics

Statistical Data of Cotton Leaf Sample Set

The total LKC in three key growth stages of cotton was measured (Table 2). The concentration gradient method (Liu et al., 2015) was used to divide the total samples into 40 calibration sets and 20 validation sets in a ratio of 2:1. The range of the K content in the calibration set including the validation set was 3.56–0.51, 2.80–0.57, and 2.30–0.40% in the three growth stages, respectively, indicating that the calibration set could well represent the entire datasets. The coefficients of variation in the calibration sets and validation sets were both between 37 and 50%, show that the LKC in cotton studied had a wide range and had good representativeness and coverage.

TABLE 2
www.frontiersin.org

Table 2. Statistical results of calibration and validation sets.

Spectral Reflectance Analysis

The single-band threshold segmentation method (Figure 3) was used to extract the average spectrum of the whole cotton leaves from hyperspectral images as the original spectrum (R) and a region of interest (ROI). The threshold was set to 0.25–0.3098, and the sample area lower than 0.25 was the background plate.

FIGURE 3
www.frontiersin.org

Figure 3. ROI extraction process by single band threshold segmentation. (A–D) Budding stage. (E–H) Flowering stage. (I–L) Boll setting stage. (A,E,I) The true color picture synthesized by ENVI cannot completely represent the true color of leaf. (B,F,J) The grayscale image of 800 nm. This grayscale image can display the outline of cotton leaves more completely in terms of brightness and clarity, so 800 nm is selected as the segmentation band. (C,G,K) Mask image. (D,H,L) Gray histogram.

The spectral curves of cotton leaves at the three growth stages were consistent with the spectral characteristics of green plants (Figure 4). There were strong absorption peaks at 450 and 680 nm and a strong reflection peak at 550 nm. Due to the many cavities in the mesophyll sponge structure, the reflectance increases sharply near the red edge region (690–760 nm), and a highly reflective platform appears in the NIR region (760–950 nm). Among them, the spectral reflectance of cotton leaves at different growth stages differed significantly in the NIR region, which showed as boll setting stage > flowering stage > budding stage. This may be because after mid-July, the redistribution and utilization of the K nutrients during fruit development and the fluidity of potassium make the K in leaves gradually transfer to cotton bolls at the flowering and boll setting stages (Singh et al., 2019), leading to the decrease in LKC, while the lack of the K in leaves increases the thickness of leaves, and the palisade tissue and parenchymal cells shrink and partially break (Zhao et al., 2001; Ramírez-Soler et al., 2021). Finally, the spectral reflectance of the boll setting stage was higher than that of the budding stage in the NIR region.

FIGURE 4
www.frontiersin.org

Figure 4. R spectral curves of the LKC in cotton at different growth stages.

Correlation Analysis Between CWT Spectra and the LKC

The R spectrum of cotton leaves was decomposed by the CWT at 10 scales. Correlation analysis was performed between the wavelet coefficients generated under each decomposition scale and the LKC, and the results were expressed as the absolute value of the correlation coefficient (|R|). The correlation between the wavelet coefficient and the LKC was relatively high, especially at the flowering stage of cotton growth (Figure 5). Under the calculation of different scales and movement factors, the regions with high correlation at the budding stage are mainly focused in the range of 500–550 and 640–660 nm on the mesoscales 3–5 and 7, with the highest correlation |R| = 0.74 (Figure 5A). The flowering stage showed obvious regional distribution, with a high correlation in the low dimension on scale 3, with the highest correlation |R| = 0.86 (Figure 5B). At the boll setting stage, the high-correlation area was mainly distributed in mesoscales 3–6, and the highest correlation |R| = 0.70 (Figure 5C). In addition, the wavelengths of 500–530, 640–660, and 740–760 nm showed higher correlation in the three growth stages.

FIGURE 5
www.frontiersin.org

Figure 5. Correlation analysis between the LKC and CWT coefficient. (A) Budding stage. (B) Flowering stage. (C) Boll setting stage.

Characteristic Wavelength Screening of R Spectrum and CWT Spectra

To further screen out characteristic wavelengths for the rapid estimation of the LKC in cotton, reduce the analytical dimension of spectral data and highlight the timeliness and convenience of spectral monitoring, the CARS and RF algorithms were selected to screen the R spectra and the CWT spectra (scales 1–10) of the three growth periods, and the selected characteristic wavelengths did not exceed 10. In the whole spectrum, the selected characteristic wavelengths of the three growth stages were similar, but there are significant differences in the screening methods of different characteristic wavelengths. The characteristic wavelengths selected based on the CARS algorithm were evenly distributed in the range of 400–950 nm, mainly located in the visible (500 nm), red edge (700 nm), and NIR regions (900 nm) (Figures 6A–C). The characteristic wavelengths screened by the RF algorithm were mainly concentrated in the visible and NIR regions in the whole spectrum, but more characteristic wavelengths appear in the visible region (Figures 6D–F). The characteristic wavelengths located in the visible region reflect the information of leaf pigment, especially the characteristic bands distributed near the strong absorption and reflection of chlorophyll. Red edge is closely related to the physical and chemical parameters of plants, is generally used to describe the health status of plants and is affected by leaf pigment and leaf area index. The characteristic wavelengths were in the range of 800–950 nm, reflecting the structure of cotton leaves and some water absorption.

FIGURE 6
www.frontiersin.org

Figure 6. Characteristic wavelengths of the R spectrum and CWT spectra screened by CARS and RF algorithms. The red and blue dots represent the characteristic wavelengths of CWT spectra and R spectrum, respectively. (A) Budding stage, CARS. (B) Flowering stage, CARS. (C) Boll setting stage, CARS. (D) Budding stage, RF. (E) Flowering stage, RF. (F) Boll setting stage, RF.

The PLSR Model Based on R Spectrum and CWT Spectra

To explore the quantitative regression relationship between cotton leaf spectral data and the LKC, a quantitative estimation model was established to realize the quantification of spectral monitoring. Using the R spectrum and CWT spectra composed of characteristic wavelengths screened by CARS and RF as independent variables and the LKC as a dependent variable, PLSR estimation models of the LKC at different growth stages were established (Table 3). The calibration and validation sets perform differently in estimating the effect of the LKC model at different growth stages and decomposition scales. At the budding stage, the model R2cal and R2val both were >0.6136 (RF-CWT-7) and 0.2675 (CARS-CWT-4). At the flowering stage, the model R2cal and R2val were >0.7311 (RF-CWT-10) and 0.6158 (CARS-CWT-2). In the boll setting period, the model R2cal and R2val were both >0.5717 (RF-CWT-3) and 0.5430 (RF-CWT-6).

TABLE 3
www.frontiersin.org

Table 3. CARS and RF algorithms are used to screen the characteristic wavelengths of the wavelet coefficient spectra, the PLSR estimation model of cotton LKC in different growth stages is established, and emphasize the relatively better performances of these the wavelet coefficient spectra in LKC estimation.

The all model R2val and RMSEval results are shown in Figure 7. The results showed that compared with the R spectrum, the CWT spectra could significantly improve the prediction performance of the LKC (Figures 7D–F). The optimal estimation model of the R spectrum at the budding, flowering and boll growth stages was constructed using the characteristic wavelengths selected by the CARS algorithm, indicating that CARS had a better estimation performance than the RF algorithm, this is similar to previous studies (Sun et al., 2021). The R2val values were 0.6613, 0.753, and 0.6643 and the RMSEval values were 0.5292, 0.2202, and 0.3405, respectively. Using the multi-decomposition scale CWT method, the best decomposition scales of the three growth stages were found to be CWT-1, CWT-3 and CWT-9 spectrum (Table 3). The best estimation models for the LKC were CARS-CWT-1 at the budding stage, RF-CWT-3 at the flowering stage and CARS-CWT-9 at the boll setting stage. The R2val values were 0.7918, 0.79, and 0.808 and the RMSEval values were 0.368, 0.1987 and 0.2508, respectively (Figures 7A–C). Compared with the single R spectrum model, the improved R2 values at the three growth stages were 0.13, 0.04, and 0.15, respectively. Higher R2val values and lower RMSEval values indicate that these models have a good fitting degree and accuracy, and the decomposed CWT spectra can effectively extract weak information, but there are significant differences in the prediction accuracy of the CWT spectra at different decomposition scales, and the screening results of the CARS and RF algorithms also show different model effects.

FIGURE 7
www.frontiersin.org

Figure 7. PLSR estimation model of the LKC in cotton based on R and CWT spectra (validation sets). (A) CARS-Budding stage. (B) CARS-Flowering stage. (C) CARS-Boll setting stage. (D) RF-Budding stage. (E) RF-Flowering stage. (F) RF-Boll setting stage.

Image Feature Analysis

Leaf Grayscale Image Database

To find image features that could optimize the estimation model of the LKC in cotton, gray image databases of leaves at the budding, flowering, and boll setting stages were constructed according to the selected characteristic wavelengths (Figure 8). Among them, the number of effective characteristic wavelengths in the three growth stages are 10, 10, and 9, respectively. Therefore, the total number of grayscale images in budding, flowering, and boll setting stages are 60 × 10, 60 × 10 and 60 × 9, respectively.

FIGURE 8
www.frontiersin.org

Figure 8. Characteristic wavelength grayscale image databases based on CWT-1, CWT-3, and CWT-9 wavelet coefficients.

Image Feature Extraction and Determination

After the mask processing for all images, three color features and 24 texture features of each characteristic wavelength grayscale image were successively calculated. Finally, the color and texture data extracted from the three growth periods were stored in 3 × 30 and 3 × 240 matrices, and a correlation analysis was performed with the LKC, respectively. The results are shown in Figure 10, where a square represents a feature.

For color features, the correlation between MEA and VAR was higher than that with SKE at the budding, flowering, and boll setting stages. The highest correlations of color features in the three growth stages were obtained for 532 nm-VAR (R = 0.39, p < 0.01), which reached extremely significant correlation, 461 nm-MEA (R = 0.18, p < 0.05) and 522 nm-MEA (R = −0.33, p < 0.05). The results showed that the overall brightness and color distribution of images were closely related to the LKC. The texture features with a high correlation between the budding, flowering, and boll setting stage were different. The high correlation at the budding stage was mainly the CON of 471, 512, and 641 nm grayscale images in the 90° and 135° directions, and the highest correlation was 641 nm-CON-135° (R = −0.48). The results showed that the furrow depth of the cotton leaf surface at the budding stage was negatively correlated with K content. The texture features with high correlation at the flowering stage was mainly the COR of the characteristic wavelength of visible light (e.g., 476 nm) in the direction of 135°, especially 476 nm-COR-135° (R = −0.38). During the boll setting stage, highly correlated texture features were concentrated in the CON at the characteristic wavelength of NIR (e.g., 799 nm), and the highest correlation was obtained at 799 nm-VAR (CON) (R = −0.44), indicating that the uniformity of the leaf surface texture was significantly negatively correlated with the K at the flowering and boll setting stages.

As can be seen from Figure 9, the color and texture feature dimensions of each growth period have higher dimensions and contain a large amount of invalid information. Therefore, we choose the CARS algorithm with better performance and combined correlation analysis to determine the effective image features. The screening results are shown in Table 4.

FIGURE 9
www.frontiersin.org

Figure 9. Correlation analysis of the LKC with grayscale image features. (A) Color features. There were 10 × 3, 10 × 3 and 9 × 3 color features in three growth stages, respectively. (B) Texture features. There were 10 × 24, 10 × 24, and 9 × 24 texture features in three growth stages, respectively.

TABLE 4
www.frontiersin.org

Table 4. Grayscale image feature extraction results.

Model Establishment and Validation

Estimation Model of the LKC Based on Combination Feature

Table 5 shows the results of the PLSR estimation model of the LKC in cotton at different growth stages was created based on the combined characteristics of “CWT spectra + image.” Compared with the R spectrum and CWT spectra models, the fusion of image features can improve the accuracy of the K estimation model at the three growth stages, but the performance of models constructed based on CWT spectra fusion with either “color” or “texture” or with “color + texture” features is different. Among them, the best estimation model was “CWT-1 + texture” for the budding stage, “CWT-3 + color” for the flowering stage and “CWT-9 + texture” for the boll setting stage. Based on the “spectra + image” combination, the best estimation models of the LKC in cotton at the budding, flowering, and boll setting stages had 16, 14, and 12 features, respectively. Moreover, texture features contribute more to the model performance improvement than color features do. These results provide a reference for rapid and non-destructive monitoring of the K.

TABLE 5
www.frontiersin.org

Table 5. PLSR estimation models of the LKC with different combinations of features.

Comparison of Model Between Single Spectrum and “Spectrum + Image” Feature

Figure 10 shows the best inversion model of the LKC in cotton at three growth stages. The optimal LKC estimation model constructed based on the characteristic wavelengths of R spectra at budding, flowering, and boll setting stages had R2val s of 0.6613, 0.7530, and 0.6643; and RMSEvals of 0.5292, 0.2202, and 0.3405, respectively (Figures 10A–C). Based on the “CWT spectra + image” feature, the accuracy R2val values of the best LKC estimation model for the three growth stages had R2val s of 0.8652, 0.8261, and 0.8952 and the RMSEval of 0.3009, 0.1821, and 0.2019, respectively (Figures 10D,F). Compared with the single R spectrum, the model accuracy R2 increased by 0.2, 0.08, and 0.23, respectively, indicating the feasibility of estimating the LKC in cotton based on CWT analysis and image feature fusion and achieving high-precision and rapid monitoring. The best estimation of the LKC in cotton was at the boll setting stage. Furthermore, because the fusion image features at the flowering stage contributed little to the improvement of the model accuracy, it indicated that the quantitative inversion of the K in the growth stage could be satisfied only by R spectral data, and the model had high stability.

FIGURE 10
www.frontiersin.org

Figure 10. Optimal PLSR estimation model of the LKC in cotton at different growth stages. (A) Budding stage, R spectrum. (B) Flowering stage, R spectrum. (C) Boll setting stage, R spectrum. (D) Budding stage, CWT-1 + texture. (E) Flowering stage, CWT-3 + color. (F) Boll setting stage, CWT-9 + texture.

Discussion

Selection of Characteristic Wavelength

In our study, there was a high correlation between CWT spectra with the LKC at the three key growth stages of cotton. However, it should be noted that LKC in cotton is not highly correlated with wavelet coefficients in the whole band of 400–950 nm but only in some important spectral regions (Figure 5). Therefore, to reduce collinearity between spectral data dimensions and adjacent wavelengths, screening several effective wavelengths containing maximum spectral information plays an important role in reducing model complexity and improving estimation ability (Lu et al., 2019b; Ruffing et al., 2021). In this study, characteristic wavelengths of the R spectrum and CWT spectra at three growth periods were screened based on the CARS and RF algorithms (Table 3). In general, the characteristic wavelengths were mainly concentrated in the visible and NIR regions in the whole spectrum. The research found that the characteristic wavelengths of the LKC of six degraded vegetation types in the green, red and NIR regions (Peng et al., 2020). Studies have shown that the K deficiency in leaves has a significant impact on the content of photosynthetic pigments (e.g., chlorophyll, carotenoid, and lutein), and the cell structure of leaves (e.g., leaf area, leaf thickness, and cell space) (Curran, 1989; Hu et al., 2020a,b), which is a key factor affecting the light absorption and utilization of plant leaves, leading to the change in reflectance (Peuelas and Filella, 1998). When crops are subjected to the K stress, the spectral reflectance of the visible and near-infrared regions increases, while chlorophyll concentration decreases (Zhao et al., 2001). At the same time, the chlorophyll ultrastructure is significantly damaged, leaf thickness increases, palisade tissue and parenchyma cells contract, and local rupture occurs (Lu et al., 2019a). In conclusion, the K deficiency symptoms can significantly affect the absorption and reflection of light by cotton leaves, change the path of light reflection and refraction and produce different spectral reflectance curves. In addition, studies have shown that the spectral reflectance of the SWIR (1,300–2,000 nm) band of rice was sensitive to K level and significantly correlated with the LKC in rice (Lu et al., 2019a). Pimstein et al. (2011) pointed out that the SWIR (1,450 nm) reflectance was significantly correlated with LKC in wheat. Sibanda et al. (2015) also showed that the SWIR spectroscopy can be used to determine K value defense on steppe.

Estimation of the LKC Based on CWT Algorithm

The CWT can decompose hyperspectral data in the time domain and frequency domain simultaneously and estimate the physiological and biochemical components of plants by looking for the best signals at different decomposition scales. The estimation model of the LKC in cotton constructed in this study using CTW spectral data has good prediction accuracy. Based on spectral data, the accuracy R2 value of the best CWT spectral models constructed for the three growth stages was 0.13, 0.04, and 0.15 higher, than that of the R spectral data model (Figure 7). Since CWT can further continuously decompose spectral data, the decomposed wavelet coefficients can correspond to the R spectrum so as to extract subtle signals in spectral data more effectively and improve spectral monitoring accuracy (Li et al., 2019). However, it should be noted that when using the CWT method, the mother wavelet function should be selected first, rather than the commonly used mother wavelet function (Sun et al., 2021; Zhao et al., 2021). A large number of studies have shown that spectral data transformed by CWT have achieved a high accuracy in the inversion of crop nutrients, chlorophyll, and agronomic traits that is superior to models obtained by traditional conversion methods (Yue et al., 2020; Lin et al., 2021; Zhao et al., 2021). Therefore, when crops are under nutrient stress, CWT can effectively mine more complete spectral information, which has great potential in feature selection, noise elimination and weak information extraction.

Estimation Model of the LKC in Cotton

The single feature extracted from the hyperspectral image has limited abilities to estimate nutrient content. Our comparison of the model performance evaluation resulting from different feature combinations (Table 5) showed that the estimation of the LKC in cotton based on “CWT spectra + image” features had high accuracy and stability at the three key growth stages, but there were differences in the modeling results of different feature combinations. Contrary to our hypothesis, it was not the case that more feature information led to better model performance. In this study, the estimation model constructed based on the features of “CWT spectra + color + texture” did not show significantly improved accuracy. Instead, the performance of models based on “CWT spectra + color” or “CWT spectra + texture” features is improved. It indicates that when constructing the model based on the feature information of “CWT spectra + color + texture,” some invalid information is added, which is interference for improving the accuracy and stability of the prediction model (Li et al., 2019). It should be noted that although the combination feature of the best estimation model of the LKC at the flowering stage is “cwt-3 + color,” and the model accuracy R2val and RMSEval are 0.8261 and 0.1821, the absolute values of the difference between R2val R2val and RMSEval based on the “CWT-3 + color” and “CWT-3 + texture” models are 0.0249 and 0.0131 (Table 5). This slight difference may be influenced by the feature parameter selection algorithm, as can be seen from Table 4, the number of color features selected by CARS during the flowering stage is relatively a little more than in other growth periods. On the other hand, statistically speaking, the difference between the two models is negligible. Therefore, it cannot be fully stated that the combined characteristics of the best estimation model for the LKC at the flowering stage is “CWT-3 + color.” Moreover, this study aimed at extracting image features of grayscale images with feature wavelengths, and different feature selection methods may obtain different feature wavelengths (Sun et al., 2021), which will lead to extracting different image feature values. Therefore, future studies will further explore the relationship between the selection of feature wavelengths and image features.

Compared with texture, color features did not contribute significantly to the estimation of the LKC. This is similar to the results of the Jamil et al. (2015) study, which showed that in the taxonomic identification of 455 Chinese herbal medicines, single texture features were superior to color or shape features, with a recognition rate of 92%. In our study, the reason for this discrepancy may be that the image we studied was a single-band grayscale image that contained different information about color and texture characteristics. For color features, we use low-dimensional color moments composed of mean, variance, and skewness to represent the color characteristics of single-band grayscale images, and because the number of channels is 1, the number of color features obtained by calculation is relatively small. In addition, for the color features we extracted, it can also be considered as another expression of spectral information, because each grayscale image has a corresponding wavelength, and perhaps there is interference information between them, reducing the additional effect of color features. However, further investigation needs to be done to positively confirmed the claim. For texture features, we used the GLCM algorithm to extract 240, 240, and 216 high-dimensional texture features from three reproductive periods, including 4 texture features in 4 directions, and its advantage may be that a larger number of feature parameters are conducive to the selection of sensitive parameters. Further, images of different wavelengths of grayscale can clearly show the veins and mesophyll parts of cotton leaves and their degree of brightness and shade (Figure 8), while texture is another feature that can be used in plant identification to describe the vein structure or leaf's surface, and it is considered as an additional feature to better describe properties of the leaves (Jamil et al., 2015). These reasons explain, as far as possible, why when we build cotton LKC estimation models, the texture features are better than the color features, and the mechanism of the relationship between LKC and texture features needs further study. Therefore, near-range hyperspectral images with high spatial and hyperspectral resolution can provide more details (Pandey et al., 2017), the model based on “CWT spectra + image” features provide a potential method for estimating the LKC in cotton.

Previous studies mainly focused on using hyperspectral data (Das et al., 2020; Furlanetto et al., 2021) or RGB images (Ghosal et al., 2018; Sun et al., 2018) for analytical modeling of the estimation of the LKC in crop. Ge et al. (2019) constructed a low-cost, non-destructive, and high-throughput maize multiphysiological parameter (including K) estimation model based on the full spectrum band (VIS-NIR-SWIR) through PLSR and SVM methods. The results for K nutrients show that the modeling results of PLSR are similar to those of SVM, and the and the model accuracy R2val is 0.586 and 0.543, respectively. The performance of the model largely depends on the sensitivity of input parameters. Compared with their research, we established the model based on the characteristic wavelength sensitive to the LKC in cotton, rather than the full wavelength, but we do not have the SWIR region, which is a deficiency. Although our study used a PLSR model to estimate the LKC with good robustness, the analysis of hyperspectral image data for large samples needs to be further explored, especially as the continuous optimization of deep learning algorithms may be beneficial to LKC estimation (Mahajan et al., 2021; Mertens et al., 2021). Das et al. (2020) studied the content changes of eight nutrient elements (K, Na, Ca, Mg, Fe, Mn, Zn, and Cu) in rice leaves under salt stress and constructed different coupled machine learning models. The results showed that the most accurate estimation of the LKC based on the PLSR-ELNET model (r = 0.928). Liu et al. (2020a) proposed a novel ensemble-modeling framework to transform the rape canopy reflectance data of the selected bands into more distinguishable probability features and identify the N, P, and K deficiency levels using the probabilities. The overall accuracy of nutritional deficiency analysis of this framework is 80.76%, it shows a competitive advantage in severe and moderate potassium deficiency. In this regard, we will study more predictive model algorithms in future work to provide a reference for rapid and non-destructive monitoring of the LKC in cotton.

Conclusions

The CWT method and PLSR model were used to estimate the LKC in cotton, which had high spectral prediction accuracy and feasibility. The CARS and RF algorithms combined with the PLSR model were used to determine the optimal decomposition scales of CWT at three growth stages, which were CWT-1, CWT-3, and CWT-9 spectra. Also, the CARS is better than RF in Characteristics selection. Compared with the single R spectrum model, the R2 values were improved by 0.13, 0.04, and 0.15, respectively.

Compared with the best estimation model of the R spectrum and CWT spectra, the PLSR model accuracy was improved after the fusion of image features at the three growth stages. The best feature combination of estimation models was “CWT-1 + texture” at the budding stage, “CWT-3 + color” at the flowering stage and “CWT-9 + texture” at the boll setting stage. Compared with the single R spectrum model, the accuracy R2 values increased by 0.2, 0.08, and 0.24, respectively.

Using characteristic wavelength to fuse image information can optimize the performance of the LKC estimation model and improve the stability and accuracy of the model. Based on the combination of “CWT spectra + image,” the best growth stage for assessing LKC in cotton was the boll setting stage, with the feature combination of “CWT-9 spectra + texture,” and its R2val and RMSEval values were 0.90 and 0.20. However, the optimal growth stage for estimating LKC only by R spectral estimation is the flowering stage. The model did not show significantly improved prediction accuracy and had high stability after integrating image features, indicating that the quantitative estimation of the LKC based on spectral data can be satisfied in this growth stage.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

QY conceived the study, analyzed the data, and wrote the manuscript. CS finished the laboratory experiment. XC and LM trained the algorithms. ZZ and XL proofread and made comments on the manuscript. All authors discussed in the review process and contributed to the article and approved the submitted version.

Funding

This study was supported by the National Natural Science Foundation of China (42061058) and the Corps' key areas of research projects (2020AB005).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

QY thank CS for supporting and helping in the experiment. ZZ and XL for revising the article.

References

Ahmad, I., and Maathuis, F. (2014). Cellular and tissue distribution of potassium: physiological relevance, mechanisms and regulation. J. Plant Physiol. 171, 708–714. doi: 10.1016/j.jplph.2013.10.016

PubMed Abstract | CrossRef Full Text | Google Scholar

Bao, S. D. (2000). Soil and Agricultural Chemistry Analysis, 3rd Edn. Beijing: China Agriculture Press.

Barros, P., Schutze, I. X., Filho, F., Yamamoto, P. T., Fiorio, P. R., and Demattê, J. A. M. (2020). Monitoring bemisia tabaci gennadius (hemiptera: aleyrodidae) infestation in soybean using by proximal sensing. Insects 12, 47. doi: 10.3390/insects12010047

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, J. X., Li, F., Wang, R., Fan, Y. F., Muhammad, A. R., Liu, Q. L., et al. (2019). Estimation of nitrogen and carbon content from soybean leaf reflectance spectra using wavelet analysis under shade stress. Comput. Electron. Agric. 156, 482–489. doi: 10.1016/j.compag.2018.12.003

CrossRef Full Text | Google Scholar

Chen, P., Haboudane, D., Tremblay, N., Wang, J., Vigneault, P., and Li, B. (2010). New spectral indicator assessing the efficiency of crop nitrogen treatment in corn and wheat. Remote Sens. Environ. 114, 1987–1997. doi: 10.1016/j.rse.2010.04.006

CrossRef Full Text | Google Scholar

Curran, P. J. (1989). Remote sensing of foliar chemistry. Remote Sens. Environ. 30, 271–278. doi: 10.1016/0034-4257(89)90069-2

CrossRef Full Text | Google Scholar

Das, B., Manohara, K. K., Mahajan, G. R., and Sahoo, R. N. (2020). Spectroscopy based novel spectral indices, pca- and plsr-coupled machine learning models for salinity stress phenotyping of rice. Spectrochim. Acta A Mol. Biomol. Spectrosc. 229, 117983. doi: 10.1016/j.saa.2019.117983

PubMed Abstract | CrossRef Full Text | Google Scholar

El-Hendawy, S., Al-Suhaibani, N., Hassan, W., Tahir, M., and Schmidhalter, U. (2017). Hyperspectral reflectance sensing to assess the growth and photosynthetic properties of wheat cultivars exposed to different irrigation rates in an irrigated arid region. PLoS ONE 12:e0183262. doi: 10.1371/journal.pone.0183262

PubMed Abstract | CrossRef Full Text | Google Scholar

Furlanetto, R. H., Nanni, M. R., Crusiol, L. G. T., Silva, G. F. C., and Sibaldelli, R. N. R. (2021). Identification and quantification of potassium (K+) deficiency in maize plants using an unmanned aerial vehicle and visible/near-infrared semi-professional digital camera. Int. J. Remote Sens. 42, 8783–8804. doi: 10.1080/01431161.2020.1871091

CrossRef Full Text | Google Scholar

Ge, H., Xiang, H., Ma, F., Li, Z., and Du, C. (2021). Estimating plant nitrogen concentration of rice through fusing vegetation indices and color moments derived from uav-rgb images. Remote Sens. 13, 1620. doi: 10.3390/rs13091620

CrossRef Full Text | Google Scholar

Ge, Y., Atefi, A., Zhang, H., Miao, C., Ramamurthy, R. K., Sigmon, B., et al. (2019). High-throughput analysis of leaf physiological and chemical traits with VIS–NIR–SWIR spectroscopy: a case study with a maize diversity panel. Plant Methods 15, 66. doi: 10.1186/s13007-019-0450-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosal, S., Blystone, D., and Singh, A. K. (2018). An explainable deep machine vision framework for plant stress phenotyping. Proc. Nat. Acad. Sci. U. S. A. 115, 4613–4618. doi: 10.1073/pnas.1716999115

PubMed Abstract | CrossRef Full Text | Google Scholar

Gierth, M., and Maser, P. (2007). Potassium transporters in plants – involvement in k + acquisition, redistribution and homeostasis. FEBS Lett. 581, 2348–2356. doi: 10.1016/j.febslet.2007.03.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, N., Li, W., Du, C., Zhang, Z., and Wang, Z. (2020a). Predicting micronutrients of wheat using hyperspectral imaging. Food Chem. 343, 128473. doi: 10.1016/j.foodchem.2020.128473

PubMed Abstract | CrossRef Full Text | Google Scholar

Hu, W. S., Lu, Z. F., Meng, F. J., Li, X. K., Cong, R. H., Ren, T., et al. (2020b). The reduction in leaf area precedes that in photosynthesis under potassium deficiency: the importance of leaf anatomy. New Phytol. 227, 1749–1763. doi: 10.1111/nph.16644

PubMed Abstract | CrossRef Full Text | Google Scholar

Hussain, A., Sahoo, R. N., Kumar, D., and Pradhan, S. (2017). Relationship of hyperspectral reflectance indices with leaf n and p concentration, dry matter accumulation and grain yield of wheat. J. Indian Soc. Remote Sens. 45, 773–784. doi: 10.1007/s12524-016-0633-y

CrossRef Full Text | Google Scholar

Jamil, N., Hussin, N., Nordin, S., and Awang, K. (2015). Automatic plant identification: is shape the key feature? Proc. Comput. Sci. 76, 436–442. doi: 10.1016/j.procs.2015.12.287

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, J., Zhu, J., Wang, X., Cheng, T., Tian, Y., Zhu, Y., et al. (2021). Estimating the leaf nitrogen content with a new feature extracted from the ultra-high spectral and spatial resolution images in wheat. Remote Sens. 13, 739. doi: 10.3390/rs13040739

CrossRef Full Text | Google Scholar

Laddi, A., Sharma, S., Kumar, A., and Kapur, P. (2013). Classification of tea grains based upon image texture feature analysis under different illumination conditions. J. Food Eng. 115, 226–231. doi: 10.1016/j.jfoodeng.2012.10.018

CrossRef Full Text | Google Scholar

Lewis, K., Morgan, G., Frame, W. H., Fromme, D., and Nichols, R. L. (2021). Cotton yield response to soil applied potassium across the United States cotton belt. Agron. J. 113, 3600–3614. doi: 10.1002/agj2.20719

CrossRef Full Text | Google Scholar

Li, F. L., Wang, L., Liu, J., Wang, Y., and Chang, Q. (2019). Evaluation of leaf n concentration in winter wheat based on discrete wavelet transform analysis. Remote Sens. 11, 1331. doi: 10.3390/rs11111331

CrossRef Full Text | Google Scholar

Lin, D., Li, G., Zhu, Y., Liu, H., and Jiao, Q. (2021). Predicting copper content in chicory leaves using hyperspectral data with continuous wavelet transforms and partial least squares. Comput. Electron. Agric. 187, 106293. doi: 10.1016/j.compag.2021.106293

CrossRef Full Text | Google Scholar

Liu, N., Xing, Z., Zhao, R., Qiao, L., Li, M., Liu, G., et al. (2020b). Analysis of chlorophyll concentration in potato crop by coupling continuous wavelet transform and spectral variable optimization. Remote Sens. 12, 2826. doi: 10.3390/rs12172826

CrossRef Full Text | Google Scholar

Liu, S., Zhai, R., and Hui, P. (2015). “Non-destructive detection of the pH value of cold fresh pork using hyperspectral imaging technique,” in 9th IFIP WG 5.14 International Conference on Computer and Computing Technologies in Agriculture (CCTA), Vol. 478, 266–274. doi: 10.1007/978-3-319-48357-3_26

CrossRef Full Text | Google Scholar

Liu, S. S., Yang, X., Guan, Q. F., Lu, Z., and Lu, J. (2020a). An ensemble modeling framework for distinguishing nitrogen, phosphorous and potassium deficiencies in winter oilseed rape (Brassica napus L.) using hyperspectral data. Remote Sens. 12, 4060. doi: 10.3390/rs12244060

CrossRef Full Text | Google Scholar

Lu, J., Yang, T., Su, X., Qi, H., Yao, X., Cheng, T., et al. (2019a). Monitoring leaf potassium content using hyperspectral vegetation indices in rice leaves. Precis. Agric. 21, 324–348. doi: 10.1007/s11119-019-09670-w

CrossRef Full Text | Google Scholar

Lu, N., Wang, W., Zhang, Q., Li, D., and Cheng, T. (2019b). Estimation of nitrogen nutrition status in winter wheat from unmanned aerial vehicle based multi-angular multispectral imagery. Front. Plant Sci. 10:1601. doi: 10.3389/fpls.2019.01601

PubMed Abstract | CrossRef Full Text | Google Scholar

Mahajan, G. R., Das, B., Murgaokar, D., Herrmann, I., and Kulkarni, R. M. (2021). Monitoring the foliar nutrients status of mango using spectroscopy-based spectral indices and PLSR-combined machine learning models. Remote Sensing. 13, 641. doi: 10.3390/rs13040641

CrossRef Full Text | Google Scholar

Mahajan, G. R., Sahoo, R. N., Pandey, R. N., Gupta, V. K., and Kumar, D. (2014). Using hyperspectral remote sensing techniques to monitor nitrogen, phosphorus, sulphur and potassium in wheat (Triticum aestivum L.). Precis. Agric. 15, 499–522. doi: 10.1007/s11119-014-9348-7

CrossRef Full Text | Google Scholar

Men, H., Yuan, H., Shi, Y., Liu, M., and Liu, J. (2021). A residual network with attention module for hyperspectral information of recognition to trace the origin of rice. Spectrochim. Acta A Mol. Biomol. Spectrosc. 263, 120155. doi: 10.1016/j.saa.2021.120155

PubMed Abstract | CrossRef Full Text | Google Scholar

Mertens, S., Verbraeken, L., Sprenger, H., Demuynck, K., and Wuyts, N. (2021). Proximal hyperspectral imaging detects diurnal and drought-induced changes in maize physiology. Front. Plant Sci. 12:640914. doi: 10.3389/fpls.2021.640914

PubMed Abstract | CrossRef Full Text | Google Scholar

Oliveira, M. R. R., Queiroz, T. R. G., Teixeira, A. D., Moreira, L. C. J., and Leao, R. A. D. (2020). Reflectance spectrometry applied to the analysis of nitrogen and potassium deficiency in cotton. Rev. Cienc. Agron. 51, e20196705. doi: 10.5935/1806-6690.20200074

CrossRef Full Text | Google Scholar

Osco, L. P., Ramos, A. P. M., Faita Pinheiro, M. M., Moriya, É. A. S., Imai, N. N., Estrabis, N., et al. (2020a). A machine learning framework to predict nutrient content in valencia-orange leaf hyperspectral measurements. Remote Sens. 12, 906. doi: 10.3390/rs12060906

CrossRef Full Text | Google Scholar

Osco, L. P., Ramos, A. P. M., Moriya, E. A. S., de Souza, M., Marcato Junior, J., Matsubara, E. T., et al. (2020b). Improvement of leaf nitrogen content inference in valencia-orange trees applying spectral analysis algorithms in UAV mounted-sensor images-sciencedirect. Int. J. Appl. Earth Observ. Geoinform. 83, 101907. doi: 10.1016/j.jag.2019.101907

CrossRef Full Text | Google Scholar

Osco, L. P., Ramos, A. P. M., Pereira, D. R., Moriya, E. A. S., Imai, N. N., Matsubara, E. T., et al. (2019). Predicting canopy nitrogen content in citrus-trees using random forest algorithm associated to spectral vegetation indices from uav-imagery. Remote Sens. 11, 2925. doi: 10.3390/rs11242925

CrossRef Full Text | Google Scholar

Pace, P. F., Cralle, H. T., Cothren, J. T., and Senseman, S. A. (1999). Photosynthate and dry matter partitioning in short and long-season cotton cultivars. Crop Sci. 39, 1065–1069. doi: 10.2135/cropsci1999.0011183X003900040018x

CrossRef Full Text | Google Scholar

Pandey, P., Ge, Y. F., Stoerger, V., and Schnable, J. C. (2017). high throughput in vivo analysis of plant leaf chemical properties using hyperspectral imaging. Front. Plant Sci. 8:1348. doi: 10.3389/fpls.2017.01348

PubMed Abstract | CrossRef Full Text | Google Scholar

Peng, Y., Zhang, M., Xu, Z., Yang, T., Su, Y., Zhou, T., et al. (2020). Estimation of leaf nutrition status in degraded vegetation based on field survey and hyperspectral data. Sci. Rep.10, 4361. doi: 10.1038/s41598-020-61294-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Pettigrew, W. T. (2008). Potassium influences on yield and quality production for maize, wheat, soybean and cotton. Physiol. Plant. 133, 670–681. doi: 10.1111/j.1399-3054.2008.01073.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Peuelas, J., and Filella, I. (1998). Visible and near-infrared reflectance techniques for diagnosing plant physiological status. Trends Plant Sci. 3, 151–156. doi: 10.1016/S1360-1385(98)01213-8

CrossRef Full Text | Google Scholar

Pimstein, A., Karnieli, A., Bansal, S. K., and Bonfil, D. J. (2011). Exploring remotely sensed technologies for monitoring wheat potassium and phosphorus using field spectroscopy. Field Crop Res. 121, 125–135. doi: 10.1016/j.fcr.2010.12.001

CrossRef Full Text | Google Scholar

Ramírez-Soler, C. H., Magnitskiy, M., Martínez, S. S. E., Alvarez-Flórez, F., and Melgarejo, L. M. (2021). Photosynthesis, biochemical activity, and leaf anatomy of tree tomato (Solanum betaceum cav.) plants under potassium deficiency. J. Appl. Bot. Food Qual. 94, 75–81. doi: 10.5073/JABFQ.2021.094.009

CrossRef Full Text | Google Scholar

Ruffing, A. M., Anthony, S. M., Strickland, L. M., Lubkin, I., and Dietz, C. R. (2021). Identification of metal stresses in arabidopsis thaliana using hyperspectral reflectance imaging. Front. Plant Sci. 12:624656. doi: 10.3389/fpls.2021.624656

PubMed Abstract | CrossRef Full Text | Google Scholar

Severtson, D., Callow, N., Flower, K., Neuhaus, A., Olejnik, M., and Nansen, C. (2016). Unmanned aerial vehicle canopy reflectance data detects potassium deficiency and green peach aphid susceptibility in canola. Precis. Agric. 17, 659–677. doi: 10.1007/s11119-016-9442-0

CrossRef Full Text | Google Scholar

Sibanda, M., Mutanga, O., Rouget, M., and Odindi, J. (2015). Exploring the potential of in situ hyperspectral data and multivariate techniques in discriminating different fertilizer treatments in grasslands. J. Appl. Remote Sens. 9, 1931–3195. doi: 10.1107/S010876738708303X

CrossRef Full Text | Google Scholar

Singh, J., Singh, N. P., Sangwan, A. K., Gill, P. S., Gill, B. S., Kaur, S., et al. (2019). Standardization of litchi leaf nutrient composition for tissue analysis in sub tropics of India. Commun. Soil Sci. Plant Anal. 50, 2889–2898. doi: 10.1080/00103624.2019.1689253

CrossRef Full Text | Google Scholar

Sun, J., Yang, W., Zhang, M., Feng, M., and Ding, G. (2021). Estimation of water content in corn leaves using hyperspectral data based on fractional order savitzky-golay derivation coupled with wavelength selection. Comput. Electron. Agric. 182, 105989. doi: 10.1016/j.compag.2021.105989

CrossRef Full Text | Google Scholar

Sun, Y., Tong, C., He, S., Wang, K., and Chen, L. (2018). Identification of nitrogen, phosphorus, and potassium deficiencies based on temporal dynamics of leaf morphology and color. Sustainability 10, 762. doi: 10.3390/su10030762

CrossRef Full Text | Google Scholar

Tian, L. W., Qi, Y. C., Dai, L., Cui, J. P., Guo, R. S., Xu, H. J., et al. (2020). Evaluation of soil nutrient content and its distribution of cultivated land in south of Xinjiang: taking Aksu prefecture as an example. J. Nucl. Agric. Sci. 34, 214–223.

Villette, J., Cuéllar, T., Verdeil, J. L., Delrot, S., and Gaillard, I. (2020). Grapevine potassium nutrition and fruit quality in the context of climate change. Front. Plant Sci. 11:123. doi: 10.3389/fpls.2020.00123

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, A. Y., Liu, Y. H., and Sheng, J. D. (2021). Spatial-temporal variation of soil nutrients in typical drip irrigation cotton fields in Xinjiang. Tianjin Agric. Sci. 27, 86–90. doi: 10.3969/j.issn.1006-6500.2021.11.019

CrossRef Full Text | Google Scholar

Wang, Y. J., Jin, G., Li, L. Q., Liu, Y., Kalkhajeh, Y. K., Ning, J. M., et al. (2020). NIR hyperspectral imaging coupled with chemometrics for non-destructive assessment of phosphorus and potassium contents in tea leaves. Infrared Phys. Technol. 108, 103365. doi: 10.1016/j.infrared.2020.103365and

CrossRef Full Text | Google Scholar

Yang, Y., Nie, J., Kan, Z., Yang, S., Zhao, H., and Li, J. (2021). Cotton stubble detection based on wavelet decomposition and texture features. Plant Methods 17, 113. doi: 10.1186/s13007-021-00809-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, X., Si, H. Y., Cheng, T., Jia, M., Chen, Q., and Tian, Y. C. (2018). Hyperspectral estimation of canopy leaf biomass phenotype per ground area using a continuous wavelet analysis in wheat. Front. Plant Sci. 9:1360. doi: 10.3389/fpls.2018.01360

PubMed Abstract | CrossRef Full Text | Google Scholar

Yue, J., Zhou, C., Guo, W., Feng, H., and Xu, K. (2020). Estimation of winter-wheat above-ground biomass using the wavelet analysis of unmanned aerial vehicle-based digital images and hyperspectral crop canopy images. Int. J. Remote Sens. 42, 1602–1622. doi: 10.1080/01431161.2020.1826057

CrossRef Full Text | Google Scholar

Zhang, J. C., Yuan, L., Pu, R. L., Loraamm, R. W., Yang, G. J., and Wang, J. H. (2014). Comparison between wavelet spectral features and conventional spectral features in detecting yellow rust for winter wheat. Comput. Electron. Agric. 100, 79–87. doi: 10.1016/j.compag.2013.11.001

CrossRef Full Text | Google Scholar

Zhao, D., Oosterhuis, D. M., and Bednarz, C. W. (2001). Influence of potassium deficiency on photosynthesis, chlorophyll content, and chloroplast ultrastructure of cotton plants. Photosynthetica 39, 103–109. doi: 10.1023/A:1012404204910

CrossRef Full Text | Google Scholar

Zhao, R., An, L., Song, D., Li, M., and Sun, H. (2021). Detection of chlorophyll fluorescence parameters of potato leaves based on continuous wavelet transform and spectral analysis. Spectrochim. Acta A Mol. Biomol. Spectrosc. 259, 119768. doi: 10.1016/j.saa.2021.119768

PubMed Abstract | CrossRef Full Text | Google Scholar

Zheng, Y., Zhu, Q. B., Huang, M., Guo, Y., and Qin, J. W. (2017). Maize and weed classifcation using color indices with support vector data description in outdoor fields. Comput. Electron. Agric. 141, 215–222. doi: 10.1016/j.compag.2017.07.028

CrossRef Full Text | Google Scholar

Zhou, J. J., Zhang, Y. H., Han, Z. M., Liu, X. Y., Jian, Y. F., and Hu, C. G. (2021). Evaluating the performance of hyperspectral leaf reflectance to detect water stress and estimation of photosynthetic capacities. Remote Sens. 13, 2160. doi: 10.3390/rs13112160

CrossRef Full Text | Google Scholar

Zou, K. L., Ge, L. Z., Zhang, C. L., Yuan, T., and Li, W. (2019). Broccoli seedling segmentation based on support vector machine combined with color texture features. IEEE Access 7, 168565–168574. doi: 10.1109/ACCESS.2019.2954587

CrossRef Full Text | Google Scholar

Keywords: hyperspectral imaging, potassium content in leaves, continuous wavelet transform, gray level co-occurrence matrix, cotton, growth stage

Citation: Yao Q, Zhang Z, Lv X, Chen X, Ma L and Sun C (2022) Estimation Model of Potassium Content in Cotton Leaves Based on Wavelet Decomposition Spectra and Image Combination Features. Front. Plant Sci. 13:920532. doi: 10.3389/fpls.2022.920532

Received: 14 April 2022; Accepted: 07 June 2022;
Published: 13 July 2022.

Edited by:

Zhenhai Li, Shandong University of Science and Technology, China

Reviewed by:

Saeed Saadat, Soil & Water Research Institute, Iran
Abbas Atefi, California Polytechnic State University, United States

Copyright © 2022 Yao, Zhang, Lv, Chen, Ma and Sun. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ze Zhang, emhhbmd6ZTEyMjcmI3gwMDA0MDtzaHp1LmVkdS5jbg==; Xin Lv, bHV4aW4mI3gwMDA0MDtzaHp1LmVkdS5jbg==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.