Keywords

1 Introduction

Ginger (Zingiber officinale Roscoe), commonly known as its high agricultural and medical values, is cultivated widely in southwest, China [1, 2]. Leaf carotenoid concentration (LCC) is an important bio-indicator of plant physiological state for their roles in photosynthesis. It can transfer the light energy to chlorophyll-a for converting to the light energy and protect the plant photosynthetic system [3,4,5]. Being able to estimate the LCC of ginger leaves could help in assessing the physiological status of crops (a) for detecting the leaf stress and senescence [6, 7]. (b) for giving an estimation of the nutrient status indirectly [8], and (c) for providing an effective way to interpret the photosynthetic mechanism of plants and detect the water and fertilizer stress under the crop early growth stage [9]. Traditional method of LCC analysis often requires to take the plant into laboratory and measure it by extraction of ethanol. Although the experiment data is accurate, but it can’t provide the changes of LCC in real time and consume a lot of energy. In recent years, hyperspectral techniques showed a great potential for estimating LCC over a range of spatial scales due to the spectral absorbance properties [10,11,12], but it seldom used to the industrial crop such as gingers. The photosynthetic pigment of plant includes chlorophyll, carotenoid and anthocyanin, however, the carotenoid is much samller compared with other pigments, especially the chlorophyll. So the estimation of LCC faces much challenges. Recently, wavelet transform as a new signal analysis tool has showed great potential in the analysis of hyperspectral data because it can amplify the weak signals and solve the overlap of spectral absorption features between chlorophyll and carotenoid [13]. Therefore, the aim of this paper include: (1) to compare the performance of vegetation indices and wavelet feature vector in estimating LCC (2) to identify the best decomposition scale in estimating LCC in different leave layers.

2 Materials and Methods

2.1 Experimental Design and Data Collection

The leaf scale spectral reflectance measurement was carried out during the 2015 growing season in the Special Seedling Experimental Base, in Chongqing (29º11΄N, 105º50΄E), China (Fig. 1). Two different ginger varieties were selected, and three nitrogen treatments with 0 kg (N0), 23 kg (N1), and 46 kg (N2) per plant were designed for this experiment, which represented three different nitrogen levels (non-fertilized, normal fertilization and excessive fertilization).

Fig. 1.
figure 1

Field experimental base of ginger

2.2 Data Acquisition

The experimental data used in this study were collected from ginger in different leave layers at four crucial growth stages. Leaf spectral measurements were made by using an ASD Fieldspec FR spectroradiometer (Analytical Spectral Devices, Boulder, Colorado, USA) associated with a leaf clip. After that, the collected leaf samples were cut into pieces and extracted for 24 h in the dark at 22° with 80% acetone. The absorbance was measured with a UV-VIS spectrophotometer at 440, 644, and 663 nm. The LCC was determined using the equation provided by Lichtenthaler et al. [14]:

$$ {\text{LCC}} = 4.695 \times {\text{A}}_{440} - 0.268 \times {\text{Chl}}_{{{\text{a}} + {\text{b}}}} $$
(1)

Where Chla+b was chlorophyll concentration, A440 was the absorbance at 440 nm.

2.3 Calculation of Vegetation Indices

To investigate the estimation performance of LCC, five commonly used vegetation indices were selected (Table 1).

Table 1. Vegetation indices used in this study, where Rλ is the reflectance at waveband λ (nm)

2.4 Extraction of the Wavelet Feature Vector

Wavelet transform is an effective mathematic tool that can decompose the original signal into different frequency components and each component is characterized with a resolution appropriate to its decomposition scale [18,19,20]. This method is performed using dilated and shifted versions of the mother wavelet, to produce a set of wavelet basis functions by the following equation:

$$ \psi_{{a{\kern 1pt} ,b}} = \frac{1}{\sqrt a }\psi \left( {\frac{\lambda - b}{a}} \right) $$
(2)

Where \( a \) represents the scaling factor (the wavelet’s width) and \( b \) represents the shifting factor (the wavelet’s position), ensures that wavelets at every scale all have the same energy.

After that, the wavelet feature vector can be calculated using the following methods. The feature vector is calculated from:

$$ F_{j} = \sqrt {\frac{1}{K}\sum\limits_{K = 1}^{K} {{\kern 1pt} {\kern 1pt} \left| {W_{{j,{\kern 1pt} {\kern 1pt} K}} } \right|^{2} } } $$
(3)

Where \( K \) is the number of coefficients at the jth decomposition level, \( W_{{j,{\kern 1pt} k}} \) is the \( kth \) coefficient at level \( j \) and \( p \) is the maximum number of decomposition levels.

3 Results and Discussion

3.1 Estimating LCC by Vegetation Indices

The selected vegetation indices were tested to estimate the LCC in different leave layers of ginger. Their correlation analysis results were shown in Table 2. It turned out that the PSSR had the best result (R2 = 0.46, RMSE = 2.63) for estimating the LCC of lower leaf layer, however, CRI550 and CRI700 showed a poor correlation with the LCC.

Table 2. Correlation analysis between the vegetation indices and LCC

3.2 Estimating LCC by Wavelet Feature Vector

As a widely used wavelet base function, Db5 was used to decompose the hyperspectral reflectance in different leaf layers. To display more clearer, only three hyperspectral reflectance decompositions were shown in Fig. 2.

Fig. 2.
figure 2

Hyperspectral reflectance of three leaf layers decomposed by the DWTA. (A) Hyperspectral signal, (B) First level disassembled spectra, (C) Second level disassembled spectra, (D) Third level disassembled spectra, (E) Fourth level disassembled spectra, (F) Fifth level disassembled spectra, (G) Sixth level disassembled spectra, (H) Seventh level disassembled spectra

It showed that the decomposed spectral reflectance of the first, second and third decomposition scales are very intensive and sharp (Fig. 2B, C and D), which represented the noise separation processes. With the decomposition scales increased, these detailed spectral signals became smooth and varied dramatically in the visible region, that is to say, the most useful information was successfully separated from the carotenoid. There has little change of the decomposed spectral reflectance as the decomposition scale varied from four to seven, that is, the decomposed spectral signals becomes stable, the wavelet feature vector should be extracted (Fig. 2H). Table 3 summarized the performance of each wavelet feature vector for estimating the LCC. It showed that there has almost no correlation between LCC and wavelet feature vectors when the decomposition scales were at the first and second scales. This can be well explained by the signal noise at low decomposition scales. The best estimation accuracy of LCC appeared at sixth decomposion scale and the R2 reached to 0.86. Compared with the selected vegetation indices, the estimation accuracy (R2) improved 46.5%–71.1% in different leave layers.

Table 3. Correlation analysis between the wavelet feature vector and LCC

4 Conclusions

This work aimed to develop a DWTA to estimate the LCC in different leave layers. Meanwhile, five commonly used vegetation indices were selected to compare the estimation accuracy with this method. The BRI, PSND and PSSR were the best vegetation indices for estimating LCC, but their highest estimation accuracy (R2) was only 0.46. However, the wavelet feature vectors with decomposed scales from three to six were sensitive to LCC when the DWTA performed. The most sensitive wavelet feature vector for estimating the LCC was the sixth decomposition scale. The estimation accuracy (R2) improved 46.5%–71.1% compared with vegetation indices. This study proves that DWTA is an effective method for estimating LCC of ginger in different layers, but using it at canopy scale will be a challenge due to the effect of LAI and soil background. Therefore, the successful use of canopy hyperspectral reflectance to estimate LCC will be an interesting extension of this work.