Monitoring Onion Crop “Cipolla Rossa di Tropea Calabria IGP” Growth and Yield Response to Varying Nitrogen Fertilizer Application Rates Using UAV Imagery

: Remote sensing (RS) platforms such as unmanned aerial vehicles (UAVs) represent an essential source of information in precision agriculture (PA) as they are able to provide images on a daily basis and at a very high resolution. In this framework, this study aims to identify the optimal level of nitrogen (N)-based nutrients for improved productivity in an onion ﬁeld of “Cipolla Rossa di Tropea” (Tropea red onion). Following an experiment that involved the arrangement of nine plots in the onion ﬁeld in a randomized complete block design (RCBD), with three replications, three different levels of N fertilization were compared: N150 (150 kg N ha − 1 ), N180 (180 kg N ha − 1 ), and e N210 (210 kg N ha − 1 ). The crop cycle was monitored using multispectral (MS) UAV imagery, producing vigor maps and taking into account the yield of data. The soil-adjusted vegetation index (SAVI) was used to monitor the vigor of the crop. In addition, the coverage’s class onion was spatially identiﬁed using geographical object-based image classiﬁcation (GEOBIA), observing differences in SAVI values obtained in plots subjected to differentiated N fertilizer treatment. The information retrieved from the analysis of soil properties (electrical conductivity, ammonium and nitrate nitrogen), yield performance and mean SAVI index data from each ﬁeld plot showed signiﬁcant relationships between the different indicators investigated. A higher onion yield was evident in plot N180, in which SAVI values were higher based on the production data.


Introduction
The International Society of Precision Agriculture (ISPAG) defines precision agriculture (PA) as "a management that gathers, processes and analyzes temporal, spatial and individual data and combines it with other information to support management decisions according to estimated variability for improved resource use efficiency, productivity, quality, profitability and sustainability of agricultural production" [1]. PA allows environmental impact costs at the farm level to be reduced through the targeted use of inputs. In addition, using agrochemicals and fertilizers more efficiently brings benefits in terms of yield and yield quality [2]. PA can be considered a cyclic system with several steps divided into data collection, data analysis, management decisions, and management decision evaluation [3]. In these frameworks, remote sensing (RS) techniques for aerial monitoring provide several solutions to collect data and map production [4]. Among the platforms, considering the differences in terms of spatial and temporal resolution, both satellites and the latest unmanned aerial vehicles (UAVs) are commonly used in PA applications. As far as satellites are concerned, nano-satellites can provide images with a relatively high temporal resolution (3 days) and a spatial resolution of about 3 m, which is useful for PA

Study Site
The onion field is located in Campora S. Giovanni, in the municipality of Amantea (Cosenza, southern Italy, 39 • 02 52" N, 16 • 06 03" E, 6 m a.s.l) (Figure 2a). The farm is part of a consortium whose cultivated area is over 700 hectares [31], involving more than 100 farms and giving a total production of 21Mt for an economic value close to 11 million euro [32]. The onion represents a significant crop production in the Calabria region [6]. The success of this product on the national and international markets relies on its unique intrinsic characteristics, such as the organoleptic properties [33] and the presence of nutraceutical compounds that make the Tropea red onion a "functional food" [34,35]. These characteristics allowed obtaining the brand European Protected Geographical Indication "Cipolla Rossa di Tropea Calabria IGP".
Drones 2021, 5, x FOR PEER REVIEW 5 of 25 success of this product on the national and international markets relies on its unique intrinsic characteristics, such as the organoleptic properties [33] and the presence of nutraceutical compounds that make the Tropea red onion a "functional food" [34,35]. These characteristics allowed obtaining the brand European Protected Geographical Indication "Cipolla Rossa di Tropea Calabria IGP". The experimental site's soil is a Typic Xeropsamments [36], and the fields are flat (slope ≤ 3%). The main soil properties, determined at the beginning of the experiment by soil sampling (before the start of the cultivation cycle, at the end of September) and referring to the 0-20 cm top layer (Ap horizon), were as follows: pH 7.51 (in water), electrical conductivity (EC1:2) 431.11 µS cm −1 , sandy texture (USDA), soil organic matter (SOM) 1.49 %, total organic carbon (TOC) 14.9 g kg −1 , total nitrogen (TN) 0.50 %, total inorganic carbon (TIC, CaCO3-C) 35 g kg −1 , available P (Olsen) 36.18 mg kg −1 , nitrate (NO3 --N) 41.93 mg kg −1 , ammonium (NH4 + -N) 10.27 mg kg −1 .
The study site covers a surface of 0.63 hectares and is less than 300 m from the sea coastline. The field (Figure 2b) is crossed by two paths, 3-m wide, where the agricultural vehicles pass while executing phytosanitary treatments and fertilization. The experimental site's soil is a Typic Xeropsamments [36], and the fields are flat (slope ≤ 3%). The main soil properties, determined at the beginning of the experiment by soil sampling (before the start of the cultivation cycle, at the end of September) and referring to the 0-20 cm top layer (Ap horizon), were as follows: pH 7.51 (in water), electrical conductivity (EC1:2) 431.11 µS cm −1 , sandy texture (USDA), soil organic matter (SOM) 1.49 %, total organic carbon (TOC) 14.9 g kg −1 , total nitrogen (TN) 0.50 %, total  10.27 mg kg −1 . The study site covers a surface of 0.63 hectares and is less than 300 m from the sea coastline. The field (Figure 2b) is crossed by two paths, 3-m wide, where the agricultural vehicles pass while executing phytosanitary treatments and fertilization.
The plants are placed at 15 × 15 cm from each other on a 3-m-wide soil bed. Any weeds on the soil bed were removed manually. The transplant was completed at the beginning of October (autumn), while the harvest was in early February (mid-winter) (Figure 2c).

Experimental Design and Crop Management
The experimental setup consisted of field plots (5 m × 5 m each), delimited using steel stakes placed at the corners of each plot (Figure 3a,b). They were arranged in a randomized complete block design (RCBD), with three replications, in order to compare three different levels of N fertilization ( Figure 3c): N150 (150 kg N ha −1 ), the typical N dose applied by the local farmers; N180 (180 kg N ha −1 ); and N210 (210 kg N ha −1 ). The plants are placed at 15 × 15 cm from each other on a 3-m-wide soil bed. Any weeds on the soil bed were removed manually. The transplant was completed at the beginning of October (autumn), while the harvest was in early February (mid-winter) (Figure 2c).

Experimental Design and Crop Management
The experimental setup consisted of field plots (5 m × 5 m each), delimited using steel stakes placed at the corners of each plot (Figure 3a,b). They were arranged in a randomized complete block design (RCBD), with three replications, in order to compare three different levels of N fertilization ( Figure 3c): N150 (150 kg N ha −1 ), the typical N dose applied by the local farmers; N180 (180 kg N ha −1 ); and N210 (210 kg N ha −1 ). To better visually compare the three N fertilization levels, the nine test plots are highlighted in light blue for plots N150 (150 kg N ha −1 ), blue for plots N180 (180 kg N ha −1 ), and white for plots N210 (210 kg N ha −1 ).
According to traditional practices, basal fertilization was carried out before onion transplanting with 400 kg ha −1 of 20N-10P2O5-10K2O chemical fertilizer, involving all field plots. At the same time, N was provided with two top-dressing fertilizations during the cropping season, with different doses according to the treatments with ammonium nitrate (34% N) and ammonium sulfate (20.6% N).
Finally, at the end of the growing season, onions from each experimental plot were manually harvested, the relative fresh yield was determined, and soil sampling was carried out. To better visually compare the three N fertilization levels, the nine test plots are highlighted in light blue for plots N150 (150 kg N ha −1 ), blue for plots N180 (180 kg N ha −1 ), and white for plots N210 (210 kg N ha −1 ).

Soil Sampling and Statistical Analysis of Soil Parameters
According to traditional practices, basal fertilization was carried out before onion transplanting with 400 kg ha −1 of 20N-10P 2 O 5 -10K 2 O chemical fertilizer, involving all field plots. At the same time, N was provided with two top-dressing fertilizations during the cropping season, with different doses according to the treatments with ammonium nitrate (34% N) and ammonium sulfate (20.6% N). Finally, at the end of the growing season, onions from each experimental plot were manually harvested, the relative fresh yield was determined, and soil sampling was carried out.

Soil Sampling and Statistical Analysis of Soil Parameters
Soil sampling was carried out within each plot in early February at the harvest. Soil samples were collected (Ap horizon, 0-20 cm soil layer) using a manual auger from the middle of each plant's inter-row space to minimize any border and plant effect and then thoroughly mixed to form a unique composite sample representative of the replicate. Three composite samples (replicates) were taken per each treatment, thereby obtaining nine composite soil samples. In the laboratory, each sample was split into two aliquots: a representative amount of moist field soil (100 g) was promptly (within 24 h) processed for NH 4 + -N and NO 3 − -N analysis, whereas the remaining aliquot (100 g) was air-dried, sieved to pass through a 2-mm sieve, and then stored at room temperature before electrical conductivity (EC) determination. Soil EC was measured in 2-mm-sieved and air-dried soil at 25 • C in a 1:2 (w/v) soil-to-water ratio mixture (EC 1:2 at 25 • C). The soil NH 4 + -N and NO 3 − -N concentration in 2 M KCl soil extracts (1:10, w/v) was determined colorimetrically by the Berthelot reaction and Griess-Ilosvay method, respectively, by using a flow injection analysis system (FIAS 400 PerkinElmer, Inc., Richmond, CT, USA) equipped with an AS90 autosampler (PerkinElmer) and linked to a UV/Vis spectrophotometer Lambda 25 (PerkinElmer).
Soil chemical data are reported as mean values (n = 3) and were expressed on a dry-weight (DW) basis (105 • C, 24 h). Data were first tested for deviation from normality (Kolmogorov-Smirnov test) and homogeneity of within-group variances (Levene's test). A one-way analysis of variance (ANOVA) was performed to highlight the effect of soil treatments on tested soil chemical variables and onion yield. Tukey's HSD test at p < 0.05 was used for multiple pairwise comparisons of treatment means. Statistical analyzes were performed by using R software [37].

UAV and Image Data Acquisition
The UAV surveys were carried out between the early middle (October) and the end of the cultivation cycle (early February). During this period, the UAV flights were performed three times; the first on 30 October 2019, the second on 29 November 2019, and the last on 31 January 2020 (Table 1), in correspondence with the following onion growth stages: leaf development, 109 BBCH-scale phase (30 October 2019); 30% of bulb development, 403 BBCH-scale phase (29 November 2019); complete bulb growth, 409 BBCH-scale phase (31 January 2020) [38].
The first two dates correspond to the phase of the crop cycle in which there is the full development of harvestable vegetative plant parts. On the last date of the survey, at the end of January, the harvest is close.
Surveys were carried out at 30 m of flight height using a multirotor UAV (Multirotor G4 Surveying Robot-Service Drone GmbH) (Figure 4a). Multirotor G4 is equipped with six high-efficiency brushless electric motors, gimbals, and flight control with a dual 32-bit processor. The UAV was equipped with an MS camera Tetracam µ-MCA06 snap (Tetracam Inc.-Chatsworth, GA, USA). This camera has six different 1.3-megapixel CMOS sensors (4:3 format, 1280 × 1024 pixels) and an interchangeable bandpass filter housed in the unit called "master," as it can synchronize the unit called "slave" (Table 1).
In this way, each sensor shoots simultaneously, and all images of the bands are synchronized. These bands are blue (480-500 nm), green (540-560 nm), red (670-690 nm), red edge (710-730 nm), NIR1 ("master" channel, 790-810 nm), and NIR2 (880-920 nm). The MS camera was equipped with its global navigation satellite system (GNSS), the FirePoint 100 GPS. All flights were conducted at an altitude of 30 m a.g.l. under constant scene illumination and cloud-free conditions. The first two dates correspond to the phase of the crop cycle in which there is the full development of harvestable vegetative plant parts. On the last date of the survey, at the end of January, the harvest is close.
Surveys were carried out at 30 m of flight height using a multirotor UAV (Multirotor G4 Surveying Robot-Service Drone GmbH) ( Figure 4a). Multirotor G4 is equipped with six high-efficiency brushless electric motors, gimbals, and flight control with a dual 32-bit processor. The UAV was equipped with an MS camera Tetracam µ-MCA06 snap (Tetracam Inc.-Chatsworth, GA, USA). This camera has six different 1.3-megapixel CMOS sensors (4:3 format, 1280 × 1024 pixels) and an interchangeable bandpass filter housed in the unit called "master," as it can synchronize the unit called "slave" (Table 1). In this way, each sensor shoots simultaneously, and all images of the bands are synchronized. These bands are blue (480-500 nm), green (540-560 nm), red (670-690 nm), red edge (710-730 nm), NIR1 ("master" channel, 790-810 nm), and NIR2 (880-920 nm). The MS camera was equipped with its global navigation satellite system (GNSS), the FirePoint 100 GPS. All flights were conducted at an altitude of 30 m a.g.l. under constant scene illumination and cloud-free conditions.
The ground sample distance (GSD) was 1.5 cm. Ground control points (GCPs) were placed in the field to obtain a high-quality geolocated output. GCPs were made using 50 cm × 50 cm black and white polypropylene panels (Figure 4b). Each GCP's position coor- The ground sample distance (GSD) was 1.5 cm. Ground control points (GCPs) were placed in the field to obtain a high-quality geolocated output. GCPs were made using 50 cm × 50 cm black and white polypropylene panels (Figure 4b). Each GCP's position co-

Image Data Pre-Processing
To perform the conversion from digital number (DN) to reflectance value, the reflectance of grey on the field calibration target (50 cm × 50 cm polypropylene panel) was measured through the spectroradiometer Apogee Ps-300. All the images were extracted from the camera and subsequently converted from the *.RAW format to 10-bit TIFF format through PixelWrench II (version 1.2.4, Tetracam Inc-Chatsworth, GA, USA). Images were then aligned, stacked, and radiometrically corrected using Pix4D mapper Pro (version 4.3-Pix4D SA, Prilly, Switzerland). Reflectance values measured on a grey calibration panel placed in the field were extracted in correspondence with each central band and used to calibrate images. Finally, reflectance orthomosaics were generated for each band and stacked into a single multiband orthomosaic for each of the six bands (blue, green, red, red edge, NIR1, and NIR2).

Image Data Processing: Segmentation, Classification, and SAVI Maps
A GEOBIA process was developed on the orthomosaics resulting from UAV image processing to discriminate onion crops from the soil. Image data processing was performed using the eCognition Developer 10 software (Trimble GeoSpatial, Munich, Germany).
Considering the characteristics of the field that includes portions of soil well visible between individual plants and in the paths of the passage of agricultural vehicles, classification was performed to separate onions and soil. The classification was developed by exploiting only the vegetation's spectral response in the different bands, as in Modica et al. [39] and Messina et al. [6]. Firstly, the procedure of segmentation of the image was performed. Segmentation is a fundamental pre-requisite for classification/feature extraction [40] as the resulting segments (i.e., image objects containing soil or vegetation/onion plants in our case) form the basis for the subsequent classification step [41]. Furthermore, object-based approaches provide accurate results when using very-high-resolution imagery [42,43]. This step foresees the image segmentation into separate, non-overlapping regions [44], subsequently extracted in the form of vectorial objects. We implemented the segmentation step on the Multiresolution Segmentation (MRS) algorithm [45], which works by identifying single-pixel objects and merging them with closer objects following the relative homogeneity criterion [46]. This criterion depends on combining the original image's spectral and shape properties and new objects generated by the merging process. Therefore, two parameters regulated homogeneity criteria: shape and compactness. The shape parameter is set by a weight attributed to the shape of objects with respect to color. The shape parameter value can range from 0 to 0.9. Color homogeneity depends on the standard deviation of spectral colors. When choosing the parameter's value, the higher the selected value, the greater the influence of the shape with respect to the color during segmentation and vice versa [47]. Compactness is the second parameter; the importance of the shape with respect to the smoothness depends on it. It derives from the product of width and length calculated on numbers of pixels [48].
In addition to the shape and compactness parameters, the scale parameter determines the final size and dimension of the objects deriving from segmentation [40,45,46,49]. So smaller values of the scale parameter produce smaller objects and vice versa with higher values. Scale parameter indirectly defines the maximum allowed heterogeneity for the resulting image objects. In addition to the parameters indicated so far, several weights can be attributed to the input data (in this case, band layers). These parameters were chosen to perform segmentation: 0.1 for shape, 0.5 for compactness, 1 for scale parameter, applied to bands/layers provided by Tetracam µ-MCA06: blue, green, red, red edge, NIR1, and NIR2. Before choosing these parameters, some trial-and-error tests were performed, based on visual interpretation, by testing attributing different values to the segmentation parameters until the segmentation considered better was obtained. Specifically, the segmentation Drones 2021, 5, 61 9 of 24 assumed as the most helpful with the research objectives was the one generating objects capable of isolating individual plants.
For monitoring onion crops in UAV imagery, the soil-adjusted vegetation index (SAVI) was chosen to analyze the vegetative vigor [4,6,50,51]. SAVI was developed by Huete [52] to minimize the effects of soil background on the vegetation signal by using a constant soil adjustment factor, L, which can take a value between 0 and 1 depending on the level of vegetation cover inserted in the normalized difference vegetation index (NDVI) formula, according to this formula Equation (1): where ρ is the reflectance at the given wavelength.
The SAVI, being one of the soil-corrected vegetation indices [53], allows further reduction of the background contribution reflectance by making it easier to detect onions, discriminating them from background soil as shown in [6]. Initially, the SAVI was calculated on the images of the entire field and histograms were used for a preliminary analysis over the months investigated. Then, descriptive statistics were used to analyze SAVI values referring only to onions identified within the single plots. For this reason, after the segmentation step, the classification of onions and bare soil was performed. The onion crops were classified by exploiting a SAVI (L = 1) threshold value ≥ −0.03 in the October image and a SAVI (L = 0.5) threshold value > 0 in November and January images.
These threshold values were chosen as a result of some trial-and-error tests. It was judged better based on visual interpretation considering its ability to detect onions, following the same methodology shown in Modica et al. [39]. The data obtained concerning only the class "onions" was used to create a mask to be applied to each of the three SAVI maps derived from UAV images.
A geometrical mask of 0.3 m × 5 m was used to exclude any weeds from each plot in the space between the soil beds. The masks were applied to the UAV images to obtain only parts of SAVI maps concerning the onion crop.

Chemical Soil Parameters and Fresh Onion Yield
N fertilization managements significantly affected (p < 0.05) all soil chemical properties investigated (Table 2). In particular, the increase of N dose applied determined a consistent raise of soil EC, NH 4 + -N and NO 3 − -N concentrations. Regarding soil EC, at the lowest N dose applied (N150), the value recorded was equal to 286.8 µS cm −1 , while greater N doses (N180 and N210) showed higher values of +38% and +250%, respectively, for N180 and N210. Soil EC determination represents an indirect measure of soil soluble salts concentration in its liquid phase. EC expresses the availability of soluble species present in the soil as an indicator of soil salinity status, including the chemical elements essential for the plant's mineral nutrition and the saline compounds that are potentially harmful to crops, especially for the susceptibility to salinity. This parameter's value can depend on several factors related to the soil's intrinsic properties and management. In agricultural soils, factors such as the climatic trend, irrigation water quality and fertilization could significantly affect soluble salts concentration in the soil [54,55]. Among vegetable crops, the onion is a species that is highly susceptible to salinity [21,56]. For this species, a salt tolerance threshold equal to 1.2 dS m −1 was estimated [57][58][59][60].
Soil NH 4 + -N and NO 3 − -N concentrations observed in N150 treatment were the lowest retrieved and equal to 7.3 and 38.8 mg kg −1 , respectively. Increased N dose distribution augmented both mineral N forms (NH 4 + -N and NO 3 − -N) by +133% and +51% for N180 and by +414% and +334% for N210. Mineral N fertilization strategy is one of the main factors that sustain onion growth and production [61,62]. Indeed, N had a critical role during photosynthesis, is fundamental to form chlorophyll, controls numerous processes of plant metabolism, and is also an indispensable part of co-enzymes and proteins. Therefore, N has been reported as the most critical yield-limiting nutrient for many crops [20]. In any case, the fertilization strategy should be carefully thought through. In fact, excesses of N can determine accelerated growth, delay of bulb maturity, increased susceptibility to pests and diseases, reduced dry-matter content, limited shelf life, and result in low onion yield and quality [63][64][65]. Moreover, a prudent strategy for N fertilization management is required in order to avoid N losses in the environment that can cause air (N gaseous emission) and water (N leaching) pollution and ensure crop production sustainability. Therefore, soil availability for the crop could meet crop requirements in order to achieve a higher fertilizer utilization efficiency, avoid pollution and, as a result, reduce costs and obtain the highest crop profit.
Fresh onion yield was significantly affected by the fertilization management applied (p < 0.05) ( Figure 5). The highest yield was registered for N180 treatments, equal to 890 q ha −1 , while lower productions were observed for N150, equivalent to 692 q ha −1 (−21%), and N210 with 619 q ha −1 (−30%). The yield registered in the current experiment is near the productivity observed in other studies carried out within and abroad of Europe [66][67][68]. In particular, concerning fertilization management, as in our experiment, the highest yield with similar production was also obtained with a dose of 180 kg N ha −1 by Díaz-Pérez et al. [69], Rodrigues et al. [70], and Gonçalves et al. [71].
As argued by García et al. [72] and Machado and Serralheiro [73], the increase of the level of salts in the soil solution can affect water and nutrient absorption, resulting in a metabolic imbalance that can lead to the reduction of plant growth and production, as observed in the current research and also reported by Medeiros Pessoa et al. [74] and Lima and Bull [75]. Indeed, soil salinity stress can determine a stunted development of plant foliage, with a reduction of leaves per plant and leaf height necessary to achieve a high production [76][77][78]. Shannon and Grieve [79] observed that due to salinity stress, leaves change their color from brilliant green to dull blue-green caused by inhibition in the synthesis of chlorophyll, and leaf tips express burn symptoms with a reduction of photosynthesis. Among plant growth phases, in particular, bulb growth is the period most sensitive to water and salinity stress, according to Kadayifci et al. [80] and Sta-Baba et al. [78], with negative consequences on both fresh weight and product quality. Moreover, beyond its direct effect on plant physiology, salinity weakening the plant can increase susceptibility to the diseases [81,82]. Among product quality indices, salinity negatively affects the mean bulb morphology, weight, diameter, and consistency [82,83].
(p < 0.05) ( Figure 5). The highest yield was registered for N180 treatments, equal to 890 q ha −1 , while lower productions were observed for N150, equivalent to 692 q ha −1 (−21%), and N210 with 619 q ha −1 (−30%). The yield registered in the current experiment is near the productivity observed in other studies carried out within and abroad of Europe [66][67][68]. In particular, concerning fertilization management, as in our experiment, the highest yield with similar production was also obtained with a dose of 180 kg N ha −1 by Díaz-Pérez et al. [69], Rodrigues et al. [70], and Gonçalves et al. [71].  (180 kg N ha −1 ) and N210 (210 kg N ha −1 )) on fresh onion yield (mean ± St. Dev, n = 3). Significant effects due to fertilization management (Trt) are presented as F-and P-value estimated by oneway ANOVA. Different letters indicate significant differences among treatments (Tukey's HSD test, p < 0.05).
Onion yield is strictly linked to N fertilization and is directly proportional to the amount of N supplied to the crop. However, onion growth is stimulated by N fertilization up to a specific limit, beyond which no benefit is observed for the crop [65,84,85]. Among the studies regarding onion N fertilization, Jilani et al. [86] observed the highest onion yield with 200 kg N ha −1 , while other studies [65,68,[87][88][89][90] reported a range between 100 and 180 kg N ha −1 to achieve the best yield performance. According to several authors, the increment of yield related to N fertilization augments is attributable to an increased plant height, leaf area, bulk weight, and diameter [65,89,90]. Moreover, it was also reported that optimal N fertilization could increase vitamin C concentration in bulbs [90]. Furthermore, it was observed that higher than optimal N fertilizer doses reduce the crop productivity level due to the appearance of phytotoxicity phenomena on plants [91,92] and plant susceptibility to diseases [93,94], also decreasing the commercial production [69,95]. However, differences in optimal N dose could be ascribable to the different varieties, environmental conditions, soil characteristics, and management. Even so, our results agree with the trend observed from the other studies, as argued above, highlighting a positive yield trend from N150 to N180 and then decrease as the administered dose increases (N210). This evidence suggests that increasing the N dose up to 180 kg N ha −1 produces increasing benefits for plant metabolism and growth, thus stimulating crop yield of Tropea red onion; increased N dose beyond this threshold harmed production, triggering a vicious circle wherein the progressive increase of N application determined a reduction of nutrient assimilation from the crop due to the progressive increase of soil EC, and thus soil salinity, to which onion is particularly susceptible.

Soil-Adjusted Vegetation Index Maps
SAVI vigor maps of the entire field are shown in Figure 6. Considering the image of October, at the beginning of the transplant, SAVI's values are shallow, ranging from 0 to 0.15. This is mainly due to the soil's predominant presence, with a value of 0, compared to the still tiny plants with thin leaves, especially in the early stages of growth [26]. The histograms reported in Figure 7 show the frequency distribution of SAVI values as a percentage of each imagery's total values. In October, more distribution of values is at −0.1. A more homogeneous distribution in the range between −0.2 and 0 can be observed in November. In January, a large part of the values' distribution is above 0, between 0.1 and 0.5. These SAVI maps were inherent to the whole field. In November, SAVI's value ranged from values below and equal to 0 up to 0.7. The crop's presence is beginning to become evident, and its coverage relative to the soil is more significant than in October. The paths followed by agricultural vehicles are also easily identifiable.
Lastly, in January, SAVI's value range is between 0.15 and 0.9. The lowest values, close to 0 (black in the maps), are primarily found in agricultural vehicles' side paths. The medium-high values, between 0.4 and 0.9, belong to the crop's vegetation, which is prevalent concerning the soil, also considering the upcoming harvest.
The histograms reported in Figure 7 show the frequency distribution of SAVI values as a percentage of each imagery's total values. In October, more distribution of values is at −0.1. A more homogeneous distribution in the range between −0.2 and 0 can be observed in November. In January, a large part of the values' distribution is above 0, between 0.1 and 0.5. These SAVI maps were inherent to the whole field. To analyze only the onion, SAVI computed on the onion crop within the nine plots was considered. In particular, other maps were generated using the masks produced on eCognition, following the same methodology adopted in Messina et al. [6]. The masks obtained by exporting a vector file containing only the class "onions" were applied to the UAV images with the aim to focus only onions present within the plots subjected to the different fertilization treatments. Consequently, parts of the scene of the image identified as bare soil were excluded from the map (Figures 8-10).
In Figure 8, showing SAVI calculated on the onions mask in the October image, what is shown in Figure 6 is confirmed, as the crop takes on values between 0 and 0.15. To analyze only the onion, SAVI computed on the onion crop within the nine plots was considered. In particular, other maps were generated using the masks produced on eCognition, following the same methodology adopted in Messina et al. [6]. The masks obtained by exporting a vector file containing only the class "onions" were applied to the UAV images with the aim to focus only onions present within the plots subjected to the different fertilization treatments. Consequently, parts of the scene of the image identified as bare soil were excluded from the map (Figures 8-10).
In Figure 8, showing SAVI calculated on the onions mask in the October image, what is shown in Figure 6 is confirmed, as the crop takes on values between 0 and 0.15.
The SAVI mean in the single plots, as can be seen in Table 3, is close to 0. This can be explained by taking it into account during the segmentation phase. The classification algorithm was not always able to correctly separate the vegetation from the soil in the background, as shown in Jeong et al. [26] and Messina et al. [6]. For this reason, during the onion classification phase, pixels with a value equal to or very close to 0 were selected, including both vegetation, consisting of thin leaves, and portions of the underlying soil.     This affected the average SAVI value of the nine plots in the October image. This is observable in the boxplots' trend in Figure 11, in which numerous outliers represent the crop's SAVI values. The SAVI mean is similar in all the plots, considering that no differentiated fertilizer treatment has been performed yet.
In Figure 9, showing SAVI calculated on onions mask in November image, the crop takes on values between 0.15 and 0.6. The SAVI mean in the single plots, as can be seen in Table 3, is higher than that seen in October, with values between 0.11 and 0.23. In particular, the lowest value was found in plot 7 (N150), while the higher value belongs to plots 4 (N210), 5 and 9 (N180). Furthermore, it is possible to notice also, as in plots 2, 5, and 9 (N180), that the average of the SAVI is always equal or superior to 0.20. Small portions of the crop where there are higher SAVI values (progressively colored from light green to dark green in Figure 9), between 0.5 and 0.8 (except for plot 7 (N150)), are also highlighted by the outliers in the boxplots in Figure 11. The SAVI calculated on the onions mask in the January image, the crop, close to harvest, takes on values between 0.3 and 0.9, as shown in Figure 10 (between yellow and dark green), confirming a progressive increase compared to November. The SAVI mean in the single plots, shown in Table 3, ranges from 0.36 to 0.53. In detail, the lowest values, 0.36 and 0.38, were found in plots 1 (N150), 3 (N210) and 7 (N150). On the other hand, the higher value belongs to plots 8 (N210) and 9 (N180) with SAVI means equal to 0.49 and 0.53, respectively, made evident by looking at the predominant dark green in Figure 10. In particular, these very high SAVI values, starting at 0.75 and exceeding 0.9, can be inferred from the outliers and the trend of the boxplots in Figure 11.

Relationship between Nitrogen Fertilization, Onion Yield, and SAVI Mapschemical Parameters and Fresh Onion Yield
The information retrieved from the analysis of soil properties (EC, NH 4 + -N, and NO 3 − -N), yield performance, and mean SAVI index data from each field plot showed significant relationships between the different indicators investigated. In particular, the observation of the soil analysis related to the soil EC and the concentration of NH 4 + -N and NO 3 − -N confirmed that the increasing levels were going from plots N150 to plots N210, while the SAVI index, interpreting the effect from the interaction between soil properties and onion plants, was positively correlated with yields, as can be observed in Table 3.
Along the cropping cycle, in the initial phase of the crop cycle, in October (leaf development, 109 BBCH phase), before the differentiation of the fertilizer treatment, no significant differences in vegetative vigor were observed in the different plots. This is shown by a similar mean SAVI value in plots with different N fertilizer doses: 0.002 (N150), 0.005 (N180), and 0.004 (N210). In the following months, the differences became more evident (November, 30% of bulb development, 403 BBCH phase; January, complete bulb growth, 409 BBCH phase). In November, the mean SAVI value in plots N180 is 0.22 and is higher than both plots N150 (0.16) and N210 (0. 19), and in January, the trend is the same with the mean SAVI value in plots N180 being higher (0.45) than both plots N150 (0.38) and N210 (0.42). Therefore, these changes were a consequence of the fertilizer dose effects in terms of soil NH 4 + -N and NO 3 − -N concentrations, as confirmed by the increase of EC, with further repercussions on onion yields, as argued above.
In more detail, in November, a higher average SAVI was observed in plots 2, 5, and 9 (N180) compared to the other treatments. The SAVI value in plot 3 (N210) was lower than plot 1 (N150), while plots 4 and 8 were equal and lower, respectively, than plots 5 and 9 (N180) (Figure 9).
In the SAVI map of the last surveys, those of January ( Figure 10), carried out two weeks before the harvest, in two out of three cases, the average SAVI in plots 1 and 6 (N150) was equal to that of plots 3 and 4 (N210). In contrast, for plots N180, the average SAVI was the highest. This evidence was corroborated at harvest by higher yields in plots 2, 5, and 9 (N180), with a mean production of 890 q ha −1 . The second-highest mean yield was obtained in plots 1, 6, and 7 (N150), equal to 692 q ha −1 , while the lowest yield, 619 q ha −1 , was obtained in plots 3, 4, and 8 (N210). Therefore, based on the production data, a higher mean Tropea red onion yield was retrieved in plots N180 which corresponded to higher SAVI values. Concerning the N210 plots, as shown in Table 3, the high content of NH 4 + -N and NO 3 − -N initially brought benefits in terms of growth (greater area occupied by the onions in plots N210 compared to N150) but up to a certain point, beyond which no benefit is observed for the crop with adverse effects (also in terms of yield) shown by lower SAVI values in most cases (and in all cases in January) compared to plots N180 and equal to those of plots N150 (in January in two out of three cases). Among the UAV survey epochs, our study highlighted the reliability of SAVI information from November (at 30% of bulb development, 403 BBCH phase) to obtain valuable data to detect plant status and predict the future yield.
This evidence, specific for the onion, highlights the importance of monitoring onion plants during the first stages of bulb growth, considering, as argued by Shannon and Grieve [79] and Pasternak et al. [96], the higher salinity sensitivity of this species during the early growth phases due to its small and shallow rooting system.
Finally, according to Al-Gaadi et al. [97], Venancio et al. [98], and Nagy et al. [99], who found positive relations between SAVI and plant growth in potato, corn, and wheat, respectively, our study confirms the reliability of the SAVI index to give helpful information regarding plant development [6,27,28] and yield expectation also for onion crop typically characterized by sparse vegetation and a low level of soil covering.

Conclusions
This study dealt with the identification of the optimal level of nitrogen (N)-based nutrients for an improved productivity in an onion field by monitoring the entire crop cycle, producing vigor maps using MS UAV imagery. To our best knowledge, only a few studies have dealt with onion crop monitoring using RS techniques.
The combined results obtained from analysis of soil chemical properties, vigor maps, and yields obtained in differentiated N fertilizer treatment plots showed a comprehensive overview of monitoring onion growth and yield response to varying N fertilizer application rates. In particular, the analysis of some soil chemical parameters, soil EC, NH 4 + -N and NO 3 − -N concentrations showed the rise of these parameters from the plots N150, treated with the lowest N fertilizer, to the plots N210, treated with the highest dose. Changes in soil chemical parameters due to fertilization management affected the yield with the highest retrieved in N180 plots followed by N150 plots, while the lowest yield was observed in the N210 plots. Our results agree with the other studies' trend, highlighting a positive yield trend from N150 to N180 and decreasing as the administered dose increases (N210). This evidence suggests that increasing the N dose up to 180 kg N ha −1 produces increasing benefits for plant metabolism and growth, thus stimulating the crop yield of Tropea red onion. On the contrary, the dose of N210 initially brought benefits in terms of growth but up to a certain point, beyond which no benefit is observed for the crop with adverse effects also for the farm, in terms of lower yield obtained and higher dose of fertilizer used, and nevertheless for the environment.
Regarding the SAVI index used for producing vigor maps, generally, its highest average was found in plots whose yield was the highest, i.e., plots N180. Reductions in vegetation vigor in plots where the experiment included applying the highest dose are probably attributable to adverse effects caused by excess salinity, proven by high EC. In order to be able to deepen the aspects dealt with in this case study, it would also be interesting to test the proposed approach not only on onions but also on other herbaceous crops.
Our study's evidence can be helpful in the onion farms to increase their real sustainability. With this regard, one aspect concerns economic sustainability and, therefore, the improvement of efficiency in the use of fertilizers based on spatial and temporal variability in the field and obtaining the best production results both in terms of quantity and quality. In contrast, another aspect is related to environmental sustainability, and given the type of experiment subject of this work, in addition to the identification of the most appropriate dose in the fertilizer treatment for better production, the harmful effects on the environment that excessive doses of N fertilizer can cause were also considered, primarily due to the phenomenon of leaching. Therefore, our scientific contribution aims to support Tropea red onion production by increasing profitability and maintaining agroecosystem quality.