1. Introduction
In sugar beet (
Beta vulgaris L.), nitrogen plays a major role in the spreading of leaves to capture sunlight. It is considered a decisive factor in the growth rate of both leaves and the storage root [
1]. Nitrogen concentration in leaves increases in the first 70 days of the sugar beet’s growth cycle, and then decreases as the growth cycle progresses [
2]. It is known that leaf chlorophyll content is related to N status [
3], and a decrease in chlorophyll content and an acceleration of canopy senescence towards the end of the crop cycle is often reported [
4]. Low levels of N result in a pale green foliage due to low chlorophyll concentration [
5], although late N application increases chlorophyll concentration in leaves [
6].
The determination of leaf N levels in the last stage of sugar beet’s growth cycle becomes relevant since it has been demonstrated that late N incorporations or releases from the soil decrease sucrose content [
7]. Draycott et al. [
5] and Malnou et al. [
6] showed that, beyond an optimum level, N has a negative effect on sugar yield. Soils that release a lot of N late in the summer feature lower yields, given that polarization (i.e., the apparent sucrose content) in that growth stage develops inversely to N availability. Gordo-Ingelmo [
2] reported that sugar beet reacts to N fertilization increases with a larger development of leaves and roots, which in turn cause an excessive use of sucrose and an increase of nonsugars. This mainly happens in cases of excessive organic fertilization, because part of the N is released belatedly, causing a stop of root ripeness.
This crucial importance of N to multiple aspects of sugar beet growth (and other crops) has led to the development of different methods to determine N levels, including destructive chemical analyses [
8], subjective leaf color charts, fluorescence techniques (including the fluorescence excitation ratio method [
9], fluorescence spectral records [
10] and fluorescence kinetic records [
11]) and chlorophyll monitoring methods [
12,
13,
14]. Most of the fluorescence-based techniques either require costly instrumentation or are labor-intensive, in a similar way as chlorophyll meter use. In spite of this drawback, measurements conducted with chlorophyll meters have proven to be highly correlated with chemical analyses in the case of sugar beet [
6,
15].
Another accepted method for N estimation is based on the use of remote sensing image analysis [
16,
17]. In this technique, images captured at different scales with different types of sensors are used to calculate vegetation indices. These indices allow leaf chlorophyll content evaluation [
18], yield prediction [
19,
20,
21], nutrient status estimation [
22,
23], disease [
24,
25] and weed detection [
26,
27], crop management [
28] or crop growth monitoring [
29].
Vegetation indices based on narrow-band imaging spectrometers (also called hyperspectral sensors) have been proposed, which—in the case of sugar beet—have been shown to have great potential to remotely estimate LAI and canopy chlorophyll content (RMSE ≤ 10%) and to retrieve canopy nitrogen content (RMSE = 10%) [
30]. Nonetheless, they remain expensive and create very large data volumes [
31], while farmers require short-term, low-expense solutions to manage their fields. To comply with these requirements, Kawashima et al. [
32] put forward a facile and low-cost diagnostic method to assess the nutrient status of plants based on the estimation of chlorophyll content of wheat and rye leaves, using a portable color video camera and a personal computer. The authors showed that leaf chlorophyll content could be estimated with sufficient accuracy using basic equipment. This research line has been continued in recent years, with the proposal of indices calculated from RGB bands of the visible light spectrum, using inexpensive, off-the-shelf cameras [
12,
31,
33,
34,
35].
In line with aforementioned studies, the aim of the research presented herein was to verify the validity of such methodologies based on index calculation from the RGB bands of leaf images (obtained using a non-scientific-grade camera) for chlorophyll content evaluation in sugar beet leaves during the final growth stage. For this validation, in a first step, different size and color sugar beet leaves from commercial farms were collected and photographed in the laboratory, extracting RGB values. Twenty-five indices proposed in the literature were tested, and nine novel indices were proposed. In a second step, the best indices were tested in field conditions, comparing conventional chlorophyll-meter measurements with the values estimated using UAV-taken photographs of a sugar beet plot.
2. Materials and Methods
2.1. Sampling Site and Crop Management
The experiment was conducted during the sugar beet harvest season, from 18 October, 2018 until 21 October, 2018. Solar radiation conditions during the experiment (13.87, 11.29, 11.32 and 6.65 MJ·m
−2 on days 1 to 4, respectively) were obtained from an automatic weather station belonging to the SIAR (Agroclimatic Information System for Irrigation) network of the Spanish Ministry of Agriculture, Food and Environment. The leaves were taken from a commercial plantation located in Magaz de Pisuerga (41°58′1.2″N, 4°26′44.2″W, altitude 740 m.a.s.l.). A description of the experimental site is provided in previous work [
36]. The sowing, delayed due to an unusually long rainfall period, was conducted on 16 April, 2018. The average population density was 125,000 plants·ha
−1 (considering a sowing distance of 13.7 cm and a seedling emergence of 86%). Crop management practices strictly followed the Spanish Research Association for Sugar Beet Crop Improvement (AIMCRA)’s recommendations, described in [
37]. The irrigation dose was 460 mm.
Five plants of ‘Fernanda KWS’ variety were randomly selected and brought to the laboratory within 40 min.
Once at the laboratory, seven leaves were randomly selected from each plant in order to obtain a representative sample in terms of sizes and colorings.
Figure 1 (left) shows some of the leaves used in the experiment.
2.2. Chlorophyll Measurements
For each leaf, 10 instant and nondestructive measurements, with an area of 0.71 cm2, were taken using a CCM-200 plus (Opti-Sciences; Hudson, NH, USA) optic chlorophyll meter. This apparatus measures the absorbance at λ = 653 nm and λ = 931 nm and calculates a chlorophyll concentration index (CCI) value that is proportional to the amount of chlorophyll in the sample. The average of the 10 measurements was considered as the chlorophyll content in the leaf at each moment.
2.3. RGB Data Acquisition
Each leaf was vertically and independently placed on a 50 × 30 cm polyurethane sheet (
Figure 1, right), and pictures were taken during 4 days at different times (at 11:00, 13:30 and 17:00 on 18 October; at 10:00 on 19 October and 20 October; and at 10:00, 14:00 and 16:00 on 21 October) in a laboratory room, using only natural light coming from a single north-facing window covered with a thin translucent curtain. Photographs were taken in vertical position with a Sony α55 (SLT-A55V, Sony, Tokyo, Japan) camera with an APS-CCMOS sensor of 16.2 MPx resolution and a Sony SAL 55-200 mm lens. The camera was mounted on a tripod in a fixed position (1 m height and 2.5 m far from the leaf). Every picture was taken with 4912 × 3264 pixels resolution, 55 mm focal distance, ISO-1600 sensitivity setting, using the aperture-priority mode with white balance set to manual and punctual measurement of exposition. Before each photograph was taken, a grey/white Lastolite Ezybalance card (Manfrotto, Cassola, Italy) was placed in front of the leaf to be photographed and used for exposure and color correction. It should be clarified that Ezybalance works as a reflectance standard of 18% of brightness [
38].
2.4. Image Processing and Color Indices Calculation
Photographs were saved in ARW format, and were analyzed using Adobe Photoshop v.14.1 (Adobe Systems Inc., San Jose, CA, USA). Besides RGB levels, this software offers a channel for the brightness scale in the histogram tool. A routine was designed to remove the background of each photograph, using several functionalities of the program to keep only the leaf pixels. The average RGB value for the total surface of each leaf was calculated by means of the histogram tool (
Figure 2). Accordingly, each leaf presented specific RGB coordinates per shoot, based on which the different color indices were calculated.
Pictures were analyzed according to three main groups of vegetation indices: indices reported in the literature (subdivided into those previously studied by Kawashima et al. [
32] and those from other sources), new indices obtained by means of stepwise linear regression (SLR) and new indices obtained by principal component analysis (PCA) (
Table 1). The indices were calculated for different datasets: for each shoot, for all the shoots taken on the same day and for the mean of daily shoots. They were also calculated for three global datasets: for all the eight shoots together, for the means of all days and, finally, for the mean of the darkest shoot (18 October 2018, 17:00) and for the brightest one (21 October 2018, 14:00). Besides this, the entire dataset was divided into two subdatasets for validation purposes, as discussed below.
New Indices and Validation
In order to explore the possibility of obtaining new indices, PCA and SLR methodologies were used. In a previous step, to ensure an appropriate validation of the calculated indices, the dataset was divided into two subdatasets, the first one was used to generate the parameters of the new indices and the second one was used to estimate the statistical error. Subdatasets were defined by choosing shoots alternately. Thus, the first subdataset included shoots #1 and #3 from 18 October, shoot #1 from 20 October and shoot #2 from 21 October; the second subdataset consisted of photographs from shoot #2 from 18 October, shoot #1 from 19 October, and shoots #1 and #3 from 21 October.
PCA indices: PCA multivariate procedure was conducted with IBM SPSS Statics 21 (IBM Corp., Armonk, NY, USA).
PCA1 and
PCA2 indices (
Table 1) were obtained as the values of eigenvectors associated with the highest eigenvalues, i.e., the principal components (PC
1) of the datasets, which explained the largest amount of variance in the data, according to the procedure described by Saberioon et al. [
12]. Results from the PCA analysis for
PCA1 and
PCA2 are presented in the supporting information file.
I1 index was derived from the other PCA indices.
SLR indices: the SLR functionality in R software (v.2.15.3, R Development Core Team, Vienna, Austria [
46]) was used as an alternative procedure to obtain new indices. Some literature indices were considered as inputs, and resultant indices were calculated by adding and removing variables one by one [
44], using AIC (Akaike Information Criterion) method for variable selection [
47]. Indices
SLR1 to
SLR5 and
I2 (
Table 1) were obtained using this procedure.
2.5. Field Test
After the in-lab assessment of the 34 RGB indices, the 4 vegetation indices for which the best performance was found were tested in real field conditions, in a sugar beet plot located in Soto de Cerrato, Palencia, Spain (41°56′23″N, 4°26′78″W, altitude 730 m.a.s.l.) on 8 November, 2019. Images were captured with a low-cost Xiaomi Mi Drone (1080P version), with f/2.6 aperture and a 1/2.3 inches CMOS sensor (
Figure 3, left).
The soil was removed using Adobe Photoshop (in a similar fashion to the background removal in the in-lab photographs), and average RGB values of photographs taken at different heights were used for the calculation of the indices.
For comparison purposes, 500 chlorophyll measurements were randomly conducted in the same sugar beet plot with the CCM-200 plus chlorophyll meter (
Figure 3, right).
2.6. Statistical Analysis
STATISTICA v. 8.0 (TIBCO, Palo Alto, CA, USA) was used for descriptive analysis, including maximum, minimum, mean, standard error and coefficients of variation (CV) of the different chlorophyll concentrations and for all the indices. Correlation matrices—with Pearson correlation coefficients (R) and p-values for each index—were also calculated for each shoot, for all the shoots in a day, for the mean of daily shoots and, finally, for three global subdatasets (i.e., for all the eight shoots together, for the means of each day together, and for the means of one shoot of the first day and one shoot of the last day).
Graphs to allow comparisons between different indices and chlorophyll concentrations were also plotted, including regression equations and coefficients of determination (R2).
For validation of the new indices, the root-mean-square error (RMSE) was calculated as follows:
where M
i = measurement with CCM-200, P
i = predicted by index, and n = number of observations. Units are CCI.
4. Conclusions
Thirty-four RGB vegetation indices were evaluated for chlorophyll content estimation in spring-sown sugar beet leaves at the final stage of the cultivation period. Based on correlations with spectroscopically determined chlorophyll contents over four days, two indices from the scientific literature, i.e., (R−B)/(R+G+B) and IPCA, and two novel indices, I1 and I2, were selected as the most suitable choices. The four selected indices showed a good performance for estimating the chlorophyll content on the basis of RGB information from photographs taken in the laboratory under natural light in different conditions, allowing to compare measurements conducted at different hours in a day and on different days. However, the new proposed RGB indices (I1 and I2) were found to improve the already good performance of the two vegetation indices selected from the literature, especially for leaves with low chlorophyll contents. In a first field validation, conducted in nonoptimized conditions, I2 index obtained the best estimation of chlorophyll content from UAV-taken photographs, followed by (R−B)/(R+G+B). The feasibility of this kind of analysis for sugar beet, which involves inexpensive, conventional off-the-shelf cameras and commonly used image-processing software, was confirmed. The proposed method can pave the way to estimate in an indirect manner the N status and the sucrose content decrease in the senescence period in a facile manner, although further studies are needed to confirm this point.