Acessibilidade / Reportar erro

CORRELATIONS AMONG VEGETATION INDICES AND PEANUT TRAITS DURING DIFFERENT CROP DEVELOPMENT STAGES

ABSTRACT

Active optical sensors have been widely used for the spatial and temporal monitoring of peanut culture because they are accurate, non-destructive methods for rapidly obtaining data. The objective of this study was to determine the optimal stage of crop growth for collecting sensor readings based on correlations between quality indicators. In addition, we compared vegetation indices (Normalized Difference Vegetation Index [NDVI], Normalized Difference Red-Edge Index, [NDRE], and Inverse Ratio Vegetation Index, [IRVI]) by monitoring temporal variability in the peanut crop in order to determine which of them obtained the best reading quality throughout the process. The experiment was performed on the 2016/17 crop in the agricultural area of the municipality of Dumont in the state of São Paulo, Brazil. The experimental design was based on the basic assumptions of statistical quality control and contained 63 sample points in a 30 × 30 m grid. The parameters were evaluated at 30, 45, 60, 75, and 119 days after sowing (DAS) using proximal sensing with GreenSeeker and OptRX sensors. We found that 45 and 60 DAS were the optimal times for monitoring peanut crop variability. For spatiotemporal monitoring of the culture with control charts, NDRE showed the best readings throughout the process when compared to NDVI and IRVI.

KEYWORDS
remote sensing; Arachis hypogaea L.; control charts; precision agriculture

INTRODUCTION

The largest peanut-producing state in Brazil is São Paulo. Production is concentrated in the region of Ribeirão Preto, where cultivation is usually done in rented areas using rotation schemes and crop succession, often to facilitate the renewal of sugarcane plantations (Barbosa et al., 2014Barbosa RM, Homem BFM, Tarsitano MAA (2014) Custo de produção e lucratividade da cultura do amendoim no município de Jaboticabal, São Paulo. Revista Ceres 61(4):475-481. DOI: http://dx.doi.org/10.1590/0034-737X201461040005
http://dx.doi.org/10.1590/0034-737X20146...
).

Peanut is considered one of the most important legumes, not only for its economic value but also nutritionally. Additionally, it is widely used in crop rotation and succession, particularly in sugarcane and pasture reforestation areas, because it is a short-cycle crop and its operations are fully mechanized (Grotta et al., 2008Grotta DCC, Furlani CEA, Silva RP, Reis GN, Cortez JW, Alves PJ (2008) Influência da profundidade de semeadura e da compactação do solo sobre a semente na produtividade do amendoim. Ciência e Agrotecnologia 32(2):547-552. DOI: http://dx.doi.org/10.1590/S1413-70542008000200031
http://dx.doi.org/10.1590/S1413-70542008...
).

Considering the economic significance of this crop, it is important to increase its productivity. This can be achieved through the use of modern techniques and methods that allow greater knowledge of crop status by providing accurate temporal monitoring. According to Grohs et al. (2009)Grohs DS, Bredemeier C, Mundstock CM, Poletto N (2009) Modelo para estimativa do potencial produtivo em trigo e cevada por meio do sensor GreenSeeker. Engenharia Agrícola 29(1): 101-112. DOI: http://dx.doi.org/10.1590/S0100-69162009000100011
http://dx.doi.org/10.1590/S0100-69162009...
, within a given crop there are areas with different productivity potentials that need different types of management.

Owing to population growth, agricultural production is expected to double by 2050 in order to meet the food demand. Precision agriculture (PA) is the key to improve resource efficiency and productivity in order to help achieve this goal under the various constraints encountered in agriculture, such as soil degradation, rising costs, climate change, lack of labor, and limited availability of arable land. To overcome these challenges, PA develops and uses sensing methodologies to provide information on crop health indicators and stages of growth (Narvaez et al., 2017Narvaez FY, Reina G, Torres-Torriti M, Kantor G, Cheein FA (2017) A Survey of Ranging and Imaging Techniques for Precision Agriculture Phenotyping. IEEE/ASME Transactions on Mechatronics 22(6):2428-2439. DOI: http://dx.doi.org/10.1109/TMECH.2017.2760866
http://dx.doi.org/10.1109/TMECH.2017.276...
).

Remote-sensing techniques appear to have a high potential for both data processing and collection in agricultural areas. These data can be obtained by field radiometry, aerial photographs, and satellite images, and can accurately provide information on field variability (Motomiya et al., 2012Motomiya AVA, Molin JP, Motomiya WR, Baio FHR (2012) Mapeamento do índice de vegetação da diferença normalizada em lavoura de algodão. Pesquisa Agropecuária Tropical 42(1): 112-118. DOI: http://dx.doi.org/10.1590/S1983-40632012000100016
http://dx.doi.org/10.1590/S1983-40632012...
). According to Amaral et al. (2015a)Amaral LR, Molin JP, Portz G, Finazzi FB, Cortinove L (2015a) Comparison of crop canopy reflectance sensors used to identify sugarcane biomass and nitrogen status. Precision Agriculture 16(1): 15—28. DOI: https://doi.org/10.1007/s11119-014-9377-2
https://doi.org/10.1007/s11119-014-9377-...
, many active optical sensors can be found, but there is little research comparing the efficiency of these devices for the determination of the parameters of cultures. Using terrestrial remote-sensing with the Normalized Difference Vegetation Index (NDVI) associated with agronomic characteristics of the peanut crop, it was previously verified that results obtained with the GreenSeeker sensor can be used to determine the population characteristics of plant cover and yield estimates, and as an additional tool for farmers to detect crop potential (Zerbato et al., 2016Zerbato C, Rosalen DL, Furlani CEA, Deghaid J, Voltarelli MA (2016) Agronomic characteristics associated with the normalized difference vegetation index (NDVI) in the peanut crop. Australian Journal of Crop Science 10(5):758-764. DOI: http://dx.doi.org/10.21475/ajcs.2016.10.05.p7167
http://dx.doi.org/10.21475/ajcs.2016.10....
).

This study addresses a subject of current relevance in the academic area of peanut crops, as it uses different vegetation indices with active or canopy sensors to monitor plants during their development and provide increased accuracy in data collection. Statistical Quality Control (SQC) was also used to study the growth stages of peanut, a statistic that has shown good results but is “not usually applied” in precision agriculture. This statistic is being used in the area of agricultural mechanization and has yielded good results in monitoring operational quality, as observed by Carneiro et al. (2017bCarneiro FM, Furlani CEA, Gírio LAS, Silva RP, Malheiros EB (2017b) Mechanized and individualized application of N-P-K in combined operations in sugarcane. Engenharia Agrícola 37(4):676-688. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v37n4p676-688/2017
http://dx.doi.org/10.1590/1809-4430-Eng....
), Cassia et al. (2015)Cassia MT, Voltarelli MA, Silva RP, Zerbato C, Lima PH (2015) Monitoramento da operação de colheita mecanizada de sementes de soja. Revista Brasileira de Engenharia Agrícola e Ambiental 19(12):-1209-1214. DOI: http://dx.doi.org/10.1590/1807-1929/agriambi.v19n12p1209-1214
http://dx.doi.org/10.1590/1807-1929/agri...
, Menezes et al. (2018)Menezes PC, Silva RP, Carneiro FM, Gírio LAS, Oliveira MF, Voltarelli MA (2018) Can combine headers and travel speeds affect the quality of soybean harvesting operations? Revista Brasileira de Engenharia Agrícola e Ambiental 22(10):732-738., Ormond et al. (2016)Ormond ATS, Voltarelli MA, Paixão CSS, Alcântara AS, Kazama EH, Furlani CEA, Silva RP (2016) Qualidade na colheita mecanizada de milho semeado em diferentes velocidades. Revista Brasileira de Milho e Sorgo 15(3):582-593. DOI: http://dx.doi.org/10.18512/1980-6477/rbms.v15n3p582-593
http://dx.doi.org/10.18512/1980-6477/rbm...
, Tavares et al. (2015)Tavares TO, Santinato F, Silva RP, Voltarelli MA, Paixão CSS, Santinato R (2015) Qualidade do recolhimento mecanizado do café. Coffee Science 10(4):455-463., Voltarelli et al. (2015)Voltarelli MA, Silva RP, Cassia MT, Ortiz DF, Torres LS (2015) Qualidade do corte basal de cana-de-açúcar utilizando-se de três modelos de facas. Engenharia Agrícola 35(3):528-541. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v35n3p528-541/2015
http://dx.doi.org/10.1590/1809-4430-Eng....
, and others.

The objective of this study was to investigate which stage of crop growth is optimal for collecting sensor readings, applying the best correlations between the evaluated quality indicators (vegetation indices, fresh biomass (in natura) and dry matter, canopy width, chlorophyll index, yield, and maturation). In addition, vegetation index readings (NDVI, NDRE, and IRVI) were compared by monitoring temporal variability in peanut crops and observing which of them obtained the best reading quality throughout the process.

MATERIAL AND METHODS

The study of the temporal variability of vegetation indices in the peanut crop was carried out on the 2016/17 crop, in the agricultural area of the municipality of Dumont located near the coordinates 21°15′22.05′'S and 47°58′15.00″W in the state of São Paulo, Brazil.

The soil in the experimental area has a claylike texture and is classified as Red Latosol. According to the Köppen climatic classification, the climate of this region is tropical with dry winter (Aw; Alvares et al., 2013Alvares CA, Stape JL, Sentelhas PC, Gonçalves JLM, Sparovek G (2013) Köppen's climate classification map for Brazil. Meteorologische Zeitschrift 22(6):711-728. DOI: https://doi.org/10.1127/0941-2948/2013/0507
https://doi.org/10.1127/0941-2948/2013/0...
).

The experimental design was based on SQC. Sixty-three sampling points in a 30 × 30 m mesh grid were used for greater representativity of data collected regarding the biophysical characteristics of peanuts associated with NDVI, NDRE, and IRVI (Fig. 1). The behavior of the quality indicators, which were fresh and dry biomass, canopy width, chlorophyll index, maturation, yield, and productivity, was evaluated at 30, 45, 60, 75, and 119 days after sowing (DAS) using NDVI, NDRE, and IRVI by means of control charts to monitor the quality of readings of the vegetation indices for the peanut crop; this correlation was obtained from the best correlations between the quality indicators or parameters evaluated.

FIGURE 1
Location of the experimental area containing in total 63 sampling points in the municipality of Dumont, state of São Paulo, Brazil. Also, the sample points used in the experiment, spaced with a 30 × 30 m mesh.

Each sampling point was composed of two rows of 5 m in length with 0.90 m spacing between rows, comprising 9 m2 of useful area per point. All evaluations were carried out at all points or plots to monitor the temporal variability of the peanut crop.

Remote-sensing evaluations were performed using GreenSeeker® (Trimble, Sunnyvale, CA) and OptRX® (Ag Leader, 2202 South Riverside Drive, Ames, IA 50010, USA) active optical sensors. Fig. 2 shows the evaluated periods in relation to the peanut growth stages. The vegetation index readings were taken in the canopy of the plant (Grohs et al., 2011Grohs DS, Bredemeier C, Poletto N, Mundstock CM (2011) Validação de modelo para predição do potencial produtivo de trigo com sensor óptico ativo. Pesquisa Agropecuária Brasileira 46(4):446-449. DOI: http://dx.doi.org/10.1590/S0100-204X2011000400015
http://dx.doi.org/10.1590/S0100-204X2011...
) between 8:00 a.m. and 12:00 p.m. For the productivity and maturation parameters, only one evaluation was performed after peanut ripening. At 119 DAS, only NDVI, IRVI, and NDRE were collected because this evaluation was performed one day before peanut ripening (120 DAS).

FIGURE 2
Growth stage of peanuts with evaluation dates.

The GreenSeeker is an active optical sensor (Trimble brand model 500) with a reading time between 1 and 2 s. This sensor emits electromagnetic radiation with a wavelength in the red (660 nm) and near-infrared (NIR; 770 nm) bands, and the light reflected by the plant is captured by the sensor. GreenSeeker automatically calculates the NDVI (Motomiya et al., 2014Motomiya AVA, Valente IMQ, Molin JP, Motomiya WR, Biscaro GA, Jordan RA (2014) Índice de vegetação no algodoeiro sob diferentes doses de nitrogênio e regulador de crescimento. Semina: Ciências Agrárias, 35(1): 169-178. DOI: http://dx.doi.org/10.5433/1679-0359.2014v35n1p169
http://dx.doi.org/10.5433/1679-0359.2014...
; Amaral et al., 2015aAmaral LR, Molin JP, Portz G, Finazzi FB, Cortinove L (2015a) Comparison of crop canopy reflectance sensors used to identify sugarcane biomass and nitrogen status. Precision Agriculture 16(1): 15—28. DOI: https://doi.org/10.1007/s11119-014-9377-2
https://doi.org/10.1007/s11119-014-9377-...
). The other obtained index along with the NDVI was the IRVI, which is still rarely used in research.

For a more accurate reading, the manufacturer recommends the working height of this sensor to be between 0.6 and 1.2 m from the target, in this case, the plant. The sensor readings were maintained at a working height of 0.6 to 0.7 m.

The OptRX is an active optical sensor (Ag Leader, model ACS430) with a reading time of 5 Hz (five readings per second). The reading height above the canopy was 0.6 to 0.7 m, with an image reading range (0.36 to 0.42 m) of 60% of the reading height. Two sensors were used to perform the average reading in real time.

According to the manufacturer (Ag Leader), in the United States, where OptRx was created, readings of vegetation indices are usually performed in real time in conjunction with product applications, using a sensor installed in a sprayer to allow the simultaneous realization of these operations. With this sensor, it is possible to obtain NDRE and NDVI. However, for this experiment, NDRE chosen as NDVI was collected with GreenSeeker.

The Laboratory of Machines and Agricultural Mechanization (LAMMA) built the support for transportation of the OptRx sensor during the field experiments, enabling the monitoring of the entire phenological cycle of the crop, from sowing to harvest. This support took into account the model developed by Professor Brenda V. Ortiz, University of Auburn (Alabama, USA). The study used a bicycle as a support, providing transport of two active optical sensors (GreenSeeker and Crop Circle) for field evaluations (Carneiro et al., 2017aCarneiro FM, Zerbato C, Menezes PC, Girio LAS, Oliveira MF, Furlani CEA (2017a) Sensoriamento Terrestre. Revista Cultivar Máquinas 15(173):20-21.).

Table 1 shows the calculations related to each index. The evaluated vegetation indices in this study were the NDVI and IRVI, obtained by GreenSeeker, and the NDRE, obtained by OptRX.

TABLE 1
Vegetation indices.

The evaluated parameters were as follows:

  • Vegetation indices were acquired using GreenSeeker and OptRX sensors, which obtained NDVI, IRVI, and NDRE.

  • Fresh and dry biomass were collected only from the aerial part of the plant. The samples were obtained utilizing a 0.25 m2 (0.5 × 0.5 m) frame. Sampling was performed per row (5 m in length), 1 m per point, in the two adjacent rows at the same time. After placing the frames in each of the two rows on both sides of the sample point, the plants were cut to the ground with the aid of a knife. The aerial part of the peanut plant was taken and the roots were discarded. The collected plants were placed in paper bags and then weighed to obtain the wet mass (fresh biomass, Fig. 3a). After that, the samples were placed in a greenhouse with air circulation at 65 °C for 72 h (Gobbi et al., 2009Gobbi KF, Garcia R, Garcez Neto AF, Pereira OG, Ventrella MC, Rocha GC (2009) Características morfológicas, estruturais e produtividade do capim-braquiária e do amendoim forrageiro submetidos ao sombreamento. Revista Brasileira de Zootecnia 38(9): 1645-1654. DOI: http://dx.doi.org/10.1590/S1516-35982009000900002
    http://dx.doi.org/10.1590/S1516-35982009...
    ; Grohs et al., 2009Grohs DS, Bredemeier C, Mundstock CM, Poletto N (2009) Modelo para estimativa do potencial produtivo em trigo e cevada por meio do sensor GreenSeeker. Engenharia Agrícola 29(1): 101-112. DOI: http://dx.doi.org/10.1590/S0100-69162009000100011
    http://dx.doi.org/10.1590/S0100-69162009...
    ) to dry them (dry biomass, Fig. 3b). Finally, the dried samples were weighed with a semi-analytical balance (model BL 3200H) to obtain the dry mass.

  • Chlorophyll content was obtained using a Marcone chlorophyll Meter, model CCM-200 Plus. The readings were performed by randomly collecting three leaves per plot, and three readings were taken from each leaf to achieve greater accuracy.

  • Canopy width and the height of the plants were measured using a scale graded in centimeters. The height of the plants was measured from the ground level to the highest point of the plant, and the width of the canopy was measured from one end of the plant to the other.

  • The productivity of the plants was measured with a frame of 2 m2 The plants were collected in raffia bags and later taken to be trodden off manually to obtain the pods, which were then weighed. The data were transformed into kg ha- 1.

  • The maturation of the pods was determined by the hull-scrape method, in which the exocarp of the pods was scrubbed, exposing the staining of the mesocarp. The classification of the pods was performed by color according to the Peanut Maturation Classification Table, which was developed by the University of Georgia, USA (Williams & Drexler, 1981Williams EJ, Drexler JS (1981) A non-destructive method for determining peanut pod maturity. Peanut Science 8(2): 134-141. DOI: https://doi.org/10.3146/i0095-3679-8-2-15
    https://doi.org/10.3146/i0095-3679-8-2-1...
    ). In this experiment, 100 random pods were collected per sampling point to determine their maturation, and the area of maturation percentage was determined using black, brown, and orange classes (Santos et al., 2013Santos EP, Silva RP, Bertonha RS, Noronha RHF, Zerbato C (2013) Produtividade e perdas de amendoim em cinco diferentes datas de arranquio. Revista Ciência Agronômica 44(4):695-702.).

FIGURE 3
Fresh (a) and dry (b) biomass per sample.

The quality indicators were analyzed through SQC using individual value control charts, which are one of the most commonly used statistical process control tools among researchers because they allow the monitoring of data behavior over time.

The control charts for individual values are composed of the arithmetic mean of the sample and the Upper (UCL) and Lower Control (LCL) Limits, obtained from the mean and standard deviation of the values analyzed per parameter. The LCL equals the mean minus three times the standard deviation, while the UCL equals the mean plus three times the standard deviation (Toledo et al., 2008Toledo A, Tabile RA, Silva RP, Furlani CEA, Magalhães SC, Costa BO (2008) Caracterização das perdas e distribuição de cobertura vegetal em colheita mecanizada de soja. Engenharia Agrícola 28(4):710-719. DOI: http://dx.doi.org/10.1590/S0100-69162008000400011
http://dx.doi.org/10.1590/S0100-69162008...
).

Correlation analysis of quality indicators (vegetation indices, fresh and dry biomass, chlorophyll index, canopy width, maturation, and yield) was also performed in this study. These analyses were used to verify the relationship between the variables and the indicators, as well as whether there were any differences between them.

Galarça et al. (2010)Galarça SP, Lima CSM, Silveira G, Rufato AR (2010) Correlação de Pearson e análise de trilha identificando variáveis para caracterizar porta-enxerto de Pyrus communis L. Ciência e Agrotecnologia 34(4):860-869. DOI: http://dx.doi.org/10.1590/S1413-70542010000400010
http://dx.doi.org/10.1590/S1413-70542010...
and Figueiredo Filho & Silva Júnior (2009)Figueiredo Filho DB, Silva Júnior JA (2009) Desvendando os Mistérios do Coeficiente de Correlação de Pearson (r). Revista Política Hoje 18(1): 115-146. stated that the Pearson correlation coefficient (r) has values ranging from -1 to 1 and that its value suggests the strength of the relationship between the variables. A perfectly positive correlation between two variables is represented by values close to 1 (r = 1), whereas a perfectly negative correlation is represented by values close to -1 (r = -1), indicating that when one variable increases, the other decreases.

Dancey & Reidy (2006)Dancey C, Reidy J (2006) Estatística sem matemática para psicologia: usando SPSS para Windows. Porto Alegre, 608p. proposed a classification for the Pearson correlation coefficient as r = 0.10 to 0.30 (weak), r = 0.40 to 0.6 (moderate), and r = 0.70 to 1 (strong). According to Figueiredo Filho & Silva Júnior (2009)Figueiredo Filho DB, Silva Júnior JA (2009) Desvendando os Mistérios do Coeficiente de Correlação de Pearson (r). Revista Política Hoje 18(1): 115-146., the degree of linear statistical dependence between the variables is greater when the value is closer to 1 (independent of the signal), but the strength of this relationship is lower when the value is closer to zero.

All the analyses presented in this study were performed using the Minitab 16 statistical software.

RESULTS AND DISCUSSION

The correlation graphs were obtained using the Pearson coefficient in order to determine the behavior and relationships of indicators and the evaluated qualities (vegetation indices, fresh and dry biomass, chlorophyll index, canopy width, maturation, and productivity). In Table 2, the highest values of the correlation coefficient for the majority of the evaluated parameters are shown at 45 and 60 DAS, which correspond to development stages R1 and R2 of the peanut, respectively (Fig. 2). This correlation is especially strong at 45 DAS (R1 stage). As observed in Table 2, according to the Dancey & Reidy (2006)Dancey C, Reidy J (2006) Estatística sem matemática para psicologia: usando SPSS para Windows. Porto Alegre, 608p. classification, the highest correlation values were moderate (r = 0.40 to 0.60) at 45 and 60 DAS (stages R1 and R2).

TABLE 2
Pearson correlation coefficients (r) for NDVI, NDRE, and IRVI in relation to fresh and dry biomass, canopy width, yield, chlorophyll index, and maturity at 30, 45, 60, 75, and 119 DAS (n = 63).

The fresh and dry biomass provided the highest correlation values with vegetation indices when compared with the other parameters, especially the correlation between NDRE and biomass (Table 2). The biomass was the parameter that most influenced the readings of the active optical sensors, corroborating the results obtained by Amaral et al. (2015b)Amaral LR, Molin JP, Schepers JS (2015b) Algorithm for Variable-Rate Nitrogen Application in Sugarcane Based on Active Crop Canopy Sensor. Agronomy Journal 107(4) :1513-1523..

Regarding the vegetation indices, NDRE presented the highest correlations when compared to the other indices for the parameters under analysis. On sugarcane, Amaral et al. (2015a)Amaral LR, Molin JP, Portz G, Finazzi FB, Cortinove L (2015a) Comparison of crop canopy reflectance sensors used to identify sugarcane biomass and nitrogen status. Precision Agriculture 16(1): 15—28. DOI: https://doi.org/10.1007/s11119-014-9377-2
https://doi.org/10.1007/s11119-014-9377-...
and Taubinger et al. (2012)Taubinger L, Amaral LR, Molin JP (2012) Vegetation indexes from active crop canopy sensor and their potential interference factors on sugarcane. In: International Conference on Precision Agriculture. Indianapolis, Hyatt Regency. found that NDRE had less influence on the plant canopy and was more efficient at predicting biomass when compared to NDVI.

One limitation of NDVI that could explain these results is the saturation of the red electromagnetic wave band. This is due to the high absorption of chlorophyll pigments, which causes saturation in this band by the canopy of the culture. Furthermore, NDVI also presents a nonlinear relationship between the biophysical parameters (Baret & Guyot, 1991Baret F, Guyot G (1991) Potentials and limits of vegetation indexes for LAI and APAR assessment. Remote Sensing of Environment 35(2-3):161-173. DOI: https://doi.org/10.1016/0034-4257(91)90009-U
https://doi.org/10.1016/0034-4257(91)900...
). In most of the evaluated parameters at 60 and 75 DAS (R2 and R4 stages, respectively), the NDVI values were probably lower due to saturation.

Regarding the width of the canopy and the chlorophyll index, the highest correlations between the indices were at 45 DAS (stage R1), mainly in the NDRE readings. Steven et al (1990)Steven MD, Malthus TJ, Demetriades-Shah TH, Danson FM, Clark JA (1990) High spectral resolution indices for crop stress. In: Steven MD, Clark JA. Applications of remote sensing in agriculture. Butterworths, p209-228. explained that the red-edge spectrum is one of the best descriptors for the remote sensing of chlorophyll concentration because this wavelength represents the maximum slope of the vegetation reflectance spectra. The same authors mentioned that this occurs in the 680 to 750 nm bands, where the reflectance changes from very high to very low (in the region of red absorption by chlorophyll) NIR because of canopy and leaf growth.

Regarding the productivity and maturation parameters, the highest correlations were at 45 and 60 DAS (stage R1 and R2, respectively). As NDRE at 60 DAS presented a higher correlation between productivity values, it should be preferred for predicting productivity compared to other indices such as NDVI (Amaral et al., 2015bAmaral LR, Molin JP, Schepers JS (2015b) Algorithm for Variable-Rate Nitrogen Application in Sugarcane Based on Active Crop Canopy Sensor. Agronomy Journal 107(4) :1513-1523.).

Thus, given these results, the most feasible period to begin monitoring of productivity and maturation using NDRE would be 45 DAS (stage R1), because of the physiological potential of the plant. This stage is at the beginning of flowering, with high translocation of peanut fluids and nutrients to flower formation. At 30 DAS (emergence stage until the beginning of flowering), most of the parameters showed low correlations, particularly the canopy width as the plant was still small and had not yet covered the area between the rows of the crop, leaving the soil exposed. Feng et al. (2016)Feng W, Zhanga Hai-Yan, Zhanga Yuan-Shuai, Qi Shuang-Li, Heng Ya-Rong, Guo Bin-Bin, Ma Dong-Yun, Guo Tian-Cai (2016) Remote detection of canopy leaf nitrogen concentration in winter wheat by using water resistance vegetation indexes from in-situ hyperspectral data. Field Crops Research 198:238-246. DOI: https://doi.org/10.1016/jfcr.2016.08.023
https://doi.org/10.1016/jfcr.2016.08.023...
noted that the spectral reflectance of the plant canopy is affected by the crop canopy, soil, biophysical properties of the vegetation, and factors that affect the accuracy of agronomic index readings.

The data in Table 2 show that at 45 and 60 DAS, higher correlations were obtained between the parameters. The values presented lower variabilities during these periods, which can be explained by the developmental stage of the plants. In the R1 and R2 reproductive stages, the plant demands greater physiological potential and greater translocation of fluids and nutrients in its interior. Moreover, the control charts verified that lower dispersion values were found in the NDRE, with lower variability indicating higher process quality, as seen in Fig. 4.

FIGURE 4
Control chart of individual values for vegetation indices (NDVI, NDRE, and IRVI) at 30, 45, 60, 75, and 119 DAS.

The analysis regarding the temporal variability of the peanut crop was monitored using the control charts of individual values of the vegetation indices (NDVI, NDRE, and IRVI) at 30, 45, 60, 75, and 119 DAS, as shown in Fig. 4. For all the evaluated periods, NDRE presented higher process quality because of its lower variability.

Among the vegetation index readings, NDVI had the highest mean values because of the saturation problem related to biomass increase. Grohs et al (2009)Grohs DS, Bredemeier C, Mundstock CM, Poletto N (2009) Modelo para estimativa do potencial produtivo em trigo e cevada por meio do sensor GreenSeeker. Engenharia Agrícola 29(1): 101-112. DOI: http://dx.doi.org/10.1590/S0100-69162009000100011
http://dx.doi.org/10.1590/S0100-69162009...
also observed that NDVI increased with increasing biomass up to the saturation value. In the present study, the mean saturation of NDVI was 0.89 between stages R2 and R4 (at 60 and 75 DAS), as shown in Fig. 4c and 4d.

Analyzing each vegetation index separately, NDVI readings at 60 and 75 DAS (stage R2 and R4) demonstrated saturation because data were constant and showed the same behavior, as shown in Fig. 4c and 4d. The control charts represented excellent tools for the temporal monitoring of the spectral behavior of the indices because they facilitated the visualization of the moment when the saturation of the NDVI reading occurred.

Similarly, Zanzarini et al. (2013)Zanzarini FV, Pissarra TCT, Brandão FJC, Teixeira DDB (2013) Correlação espacial do índice de vegetação (NDVI) de imagem Landsat/ETM+ com atributos do solo. Revista Brasileira de Engenharia Agrícola e Ambiental 17(6):608-614. DOI: http://dx.doi.org/10.1590/S1415-43662013000600006
http://dx.doi.org/10.1590/S1415-43662013...
observed that this vegetation index has a limiting factor due to its rapid saturation at a certain stage of development due to the increase in biomass, and the consequent stability in the readings. However, at 119 DAS, one day before harvest, NDVI accurately collected the indices, probably owing to the decrease in biomass at this stage.

Regarding IRVI, Table 2 shows that its values were inversely proportional to those of NDVI. Kapp Júnior et al. (2016)Kapp Júnior C, Guimarães AM, Caires EF (2016) Use of active canopy sensors to discriminate wheat response to nitrogen fertilization under no-tillage. Engenharia Agrícola 36(5):886-894. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v36n5p886-894/2016
http://dx.doi.org/10.1590/1809-4430-Eng....
verified that the concentration of chlorophyll in the leaf tissue affects the reflection of the wavelength in the visible range. This pattern is expected because the lower the chlorophyll level and application of nitrogen in the crop, the lower the absorption of solar radiation in the region, resulting in a decrease in NDVI and an increase in the IRVI.

There was temporal variability among the indices because of the stage of peanut development, bands of readings used in the indices, and size of the plant canopy. In the initial stages, the crop had an incomplete canopy and this influenced the sensor readings due to soil reflectance. The outlying points can be explained by factors influencing the vegetation indices, such as canopy width, leaf geometry, stage of development or age of the plant, leaf color, and sensor imaging range.

Our results show that the optimum time for the use of active optical sensing from a terrestrial platform was identified. Besides, the employed method allows the monitoring of the growing stages of the culture because of its non-destructive nature and provides knowledge regarding the spatial variability of the crop through vegetation indexes, specifically showing if plants are vigorous. This will facilitate the faster detection of the causes that could affect the development of crops and consequently the yield, such as pests, diseases, and nematodes, among others. This work is very relevant because it shows the optimal time for the use of these sensors, which could give farmers greater time-saving.

CONCLUSIONS

Considering that the best correlations among the evaluated quality indicators were obtained at, 45 and 60 DAS (stage R1 and R2), particularly 45 DAS, these times should be considered the optimal time for the monitoring of peanut crop variability.

For the spatiotemporal monitoring of the culture through control charts, NDRE showed the best reading qualities throughout the process when compared with NDVI and IRVI.

The use of active optical sensors shows great potential for the detection of temporal variability in peanut crops. Furthermore, it was also possible to associate the biophysical characteristics of the crop with the vegetation indices. Through the use of Pearson correlation, it was possible to identify the parameters that best correlated with the vegetation indices.

The control charts applied in this study show promise for the temporal monitoring of vegetation index readings during the different stages of peanut growth.

This study will serve as a reference for future research regarding the use of remote sensing techniques at the terrestrial level. These approaches will also help both researchers and rural producers by identifying the optimal time for the use of active optical sensors in order to monitor peanut crop variability and verify plant vigor using a vegetation index. Finally, it also provides new knowledge in the area of physiology, irrigation with varied rate, and soil fertility application of fertilizers before and after the optimal time.

REFERENCES

  • Alvares CA, Stape JL, Sentelhas PC, Gonçalves JLM, Sparovek G (2013) Köppen's climate classification map for Brazil. Meteorologische Zeitschrift 22(6):711-728. DOI: https://doi.org/10.1127/0941-2948/2013/0507
    » https://doi.org/10.1127/0941-2948/2013/0507
  • Amaral LR, Molin JP, Portz G, Finazzi FB, Cortinove L (2015a) Comparison of crop canopy reflectance sensors used to identify sugarcane biomass and nitrogen status. Precision Agriculture 16(1): 15—28. DOI: https://doi.org/10.1007/s11119-014-9377-2
    » https://doi.org/10.1007/s11119-014-9377-2
  • Amaral LR, Molin JP, Schepers JS (2015b) Algorithm for Variable-Rate Nitrogen Application in Sugarcane Based on Active Crop Canopy Sensor. Agronomy Journal 107(4) :1513-1523.
  • Barbosa RM, Homem BFM, Tarsitano MAA (2014) Custo de produção e lucratividade da cultura do amendoim no município de Jaboticabal, São Paulo. Revista Ceres 61(4):475-481. DOI: http://dx.doi.org/10.1590/0034-737X201461040005
    » http://dx.doi.org/10.1590/0034-737X201461040005
  • Baret F, Guyot G (1991) Potentials and limits of vegetation indexes for LAI and APAR assessment. Remote Sensing of Environment 35(2-3):161-173. DOI: https://doi.org/10.1016/0034-4257(91)90009-U
    » https://doi.org/10.1016/0034-4257(91)90009-U
  • Buschmann C, Nagel E (1993) In vivo spectroscopy and internal optics of leaves as basis for remote sensing of vegetation. International Journal of Remote Sensing 14(4):711-722. DOI: https://doi.org/10.1080/01431169308904370
    » https://doi.org/10.1080/01431169308904370
  • Carneiro FM, Zerbato C, Menezes PC, Girio LAS, Oliveira MF, Furlani CEA (2017a) Sensoriamento Terrestre. Revista Cultivar Máquinas 15(173):20-21.
  • Carneiro FM, Furlani CEA, Gírio LAS, Silva RP, Malheiros EB (2017b) Mechanized and individualized application of N-P-K in combined operations in sugarcane. Engenharia Agrícola 37(4):676-688. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v37n4p676-688/2017
    » http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v37n4p676-688/2017
  • Cassia MT, Voltarelli MA, Silva RP, Zerbato C, Lima PH (2015) Monitoramento da operação de colheita mecanizada de sementes de soja. Revista Brasileira de Engenharia Agrícola e Ambiental 19(12):-1209-1214. DOI: http://dx.doi.org/10.1590/1807-1929/agriambi.v19n12p1209-1214
    » http://dx.doi.org/10.1590/1807-1929/agriambi.v19n12p1209-1214
  • Dancey C, Reidy J (2006) Estatística sem matemática para psicologia: usando SPSS para Windows. Porto Alegre, 608p.
  • Feng W, Zhanga Hai-Yan, Zhanga Yuan-Shuai, Qi Shuang-Li, Heng Ya-Rong, Guo Bin-Bin, Ma Dong-Yun, Guo Tian-Cai (2016) Remote detection of canopy leaf nitrogen concentration in winter wheat by using water resistance vegetation indexes from in-situ hyperspectral data. Field Crops Research 198:238-246. DOI: https://doi.org/10.1016/jfcr.2016.08.023
    » https://doi.org/10.1016/jfcr.2016.08.023
  • Figueiredo Filho DB, Silva Júnior JA (2009) Desvendando os Mistérios do Coeficiente de Correlação de Pearson (r). Revista Política Hoje 18(1): 115-146.
  • Galarça SP, Lima CSM, Silveira G, Rufato AR (2010) Correlação de Pearson e análise de trilha identificando variáveis para caracterizar porta-enxerto de Pyrus communis L. Ciência e Agrotecnologia 34(4):860-869. DOI: http://dx.doi.org/10.1590/S1413-70542010000400010
    » http://dx.doi.org/10.1590/S1413-70542010000400010
  • Gobbi KF, Garcia R, Garcez Neto AF, Pereira OG, Ventrella MC, Rocha GC (2009) Características morfológicas, estruturais e produtividade do capim-braquiária e do amendoim forrageiro submetidos ao sombreamento. Revista Brasileira de Zootecnia 38(9): 1645-1654. DOI: http://dx.doi.org/10.1590/S1516-35982009000900002
    » http://dx.doi.org/10.1590/S1516-35982009000900002
  • Grohs DS, Bredemeier C, Mundstock CM, Poletto N (2009) Modelo para estimativa do potencial produtivo em trigo e cevada por meio do sensor GreenSeeker. Engenharia Agrícola 29(1): 101-112. DOI: http://dx.doi.org/10.1590/S0100-69162009000100011
    » http://dx.doi.org/10.1590/S0100-69162009000100011
  • Grohs DS, Bredemeier C, Poletto N, Mundstock CM (2011) Validação de modelo para predição do potencial produtivo de trigo com sensor óptico ativo. Pesquisa Agropecuária Brasileira 46(4):446-449. DOI: http://dx.doi.org/10.1590/S0100-204X2011000400015
    » http://dx.doi.org/10.1590/S0100-204X2011000400015
  • Grotta DCC, Furlani CEA, Silva RP, Reis GN, Cortez JW, Alves PJ (2008) Influência da profundidade de semeadura e da compactação do solo sobre a semente na produtividade do amendoim. Ciência e Agrotecnologia 32(2):547-552. DOI: http://dx.doi.org/10.1590/S1413-70542008000200031
    » http://dx.doi.org/10.1590/S1413-70542008000200031
  • Kapp Júnior C, Guimarães AM, Caires EF (2016) Use of active canopy sensors to discriminate wheat response to nitrogen fertilization under no-tillage. Engenharia Agrícola 36(5):886-894. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v36n5p886-894/2016
    » http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v36n5p886-894/2016
  • Menezes PC, Silva RP, Carneiro FM, Gírio LAS, Oliveira MF, Voltarelli MA (2018) Can combine headers and travel speeds affect the quality of soybean harvesting operations? Revista Brasileira de Engenharia Agrícola e Ambiental 22(10):732-738.
  • Motomiya AVA, Molin JP, Motomiya WR, Baio FHR (2012) Mapeamento do índice de vegetação da diferença normalizada em lavoura de algodão. Pesquisa Agropecuária Tropical 42(1): 112-118. DOI: http://dx.doi.org/10.1590/S1983-40632012000100016
    » http://dx.doi.org/10.1590/S1983-40632012000100016
  • Motomiya AVA, Valente IMQ, Molin JP, Motomiya WR, Biscaro GA, Jordan RA (2014) Índice de vegetação no algodoeiro sob diferentes doses de nitrogênio e regulador de crescimento. Semina: Ciências Agrárias, 35(1): 169-178. DOI: http://dx.doi.org/10.5433/1679-0359.2014v35n1p169
    » http://dx.doi.org/10.5433/1679-0359.2014v35n1p169
  • Narvaez FY, Reina G, Torres-Torriti M, Kantor G, Cheein FA (2017) A Survey of Ranging and Imaging Techniques for Precision Agriculture Phenotyping. IEEE/ASME Transactions on Mechatronics 22(6):2428-2439. DOI: http://dx.doi.org/10.1109/TMECH.2017.2760866
    » http://dx.doi.org/10.1109/TMECH.2017.2760866
  • Ormond ATS, Voltarelli MA, Paixão CSS, Alcântara AS, Kazama EH, Furlani CEA, Silva RP (2016) Qualidade na colheita mecanizada de milho semeado em diferentes velocidades. Revista Brasileira de Milho e Sorgo 15(3):582-593. DOI: http://dx.doi.org/10.18512/1980-6477/rbms.v15n3p582-593
    » http://dx.doi.org/10.18512/1980-6477/rbms.v15n3p582-593
  • Rouse JW, Haas RH, Schell JA, Deering DW (1973) Monitoring vegetation systems in the great plains with ERTS. In: Third earth resources technology satellite-1 symposium. National Aeronautics and Space Administration v1, p309-330.
  • Santos EP, Silva RP, Bertonha RS, Noronha RHF, Zerbato C (2013) Produtividade e perdas de amendoim em cinco diferentes datas de arranquio. Revista Ciência Agronômica 44(4):695-702.
  • Steven MD, Malthus TJ, Demetriades-Shah TH, Danson FM, Clark JA (1990) High spectral resolution indices for crop stress. In: Steven MD, Clark JA. Applications of remote sensing in agriculture. Butterworths, p209-228.
  • Taubinger L, Amaral LR, Molin JP (2012) Vegetation indexes from active crop canopy sensor and their potential interference factors on sugarcane. In: International Conference on Precision Agriculture. Indianapolis, Hyatt Regency.
  • Tavares TO, Santinato F, Silva RP, Voltarelli MA, Paixão CSS, Santinato R (2015) Qualidade do recolhimento mecanizado do café. Coffee Science 10(4):455-463.
  • Toledo A, Tabile RA, Silva RP, Furlani CEA, Magalhães SC, Costa BO (2008) Caracterização das perdas e distribuição de cobertura vegetal em colheita mecanizada de soja. Engenharia Agrícola 28(4):710-719. DOI: http://dx.doi.org/10.1590/S0100-69162008000400011
    » http://dx.doi.org/10.1590/S0100-69162008000400011
  • Voltarelli MA, Silva RP, Cassia MT, Ortiz DF, Torres LS (2015) Qualidade do corte basal de cana-de-açúcar utilizando-se de três modelos de facas. Engenharia Agrícola 35(3):528-541. DOI: http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v35n3p528-541/2015
    » http://dx.doi.org/10.1590/1809-4430-Eng.Agric.v35n3p528-541/2015
  • Williams EJ, Drexler JS (1981) A non-destructive method for determining peanut pod maturity. Peanut Science 8(2): 134-141. DOI: https://doi.org/10.3146/i0095-3679-8-2-15
    » https://doi.org/10.3146/i0095-3679-8-2-15
  • Zanzarini FV, Pissarra TCT, Brandão FJC, Teixeira DDB (2013) Correlação espacial do índice de vegetação (NDVI) de imagem Landsat/ETM+ com atributos do solo. Revista Brasileira de Engenharia Agrícola e Ambiental 17(6):608-614. DOI: http://dx.doi.org/10.1590/S1415-43662013000600006
    » http://dx.doi.org/10.1590/S1415-43662013000600006
  • Zerbato C, Rosalen DL, Furlani CEA, Deghaid J, Voltarelli MA (2016) Agronomic characteristics associated with the normalized difference vegetation index (NDVI) in the peanut crop. Australian Journal of Crop Science 10(5):758-764. DOI: http://dx.doi.org/10.21475/ajcs.2016.10.05.p7167
    » http://dx.doi.org/10.21475/ajcs.2016.10.05.p7167

Edited by

Area Editor: Fabio Henrique Rojo Baio

Publication Dates

  • Publication in this collection
    09 Sept 2019
  • Date of issue
    Sept 2019

History

  • Received
    11 Mar 2019
  • Accepted
    29 May 2019
Associação Brasileira de Engenharia Agrícola SBEA - Associação Brasileira de Engenharia Agrícola, Departamento de Engenharia e Ciências Exatas FCAV/UNESP, Prof. Paulo Donato Castellane, km 5, 14884.900 | Jaboticabal - SP, Tel./Fax: +55 16 3209 7619 - Jaboticabal - SP - Brazil
E-mail: revistasbea@sbea.org.br