Discrimination of Tomato Plants (Solanum lycopersicum) Grown under Anaerobic Baffled Reactor Effluent, Nitrified Urine Concentrates and Commercial Hydroponic Fertilizer Regimes Using Simulated Sensor Spectral Settings

We assess the discriminative strength of three different satellite spectral settings (HyspIRI, the forthcoming Landsat 9 and Sentinel 2-MSI), in mapping tomato (Solanum lycopersicum Linnaeus) plants grown under hydroponic system, using human-excreta derived materials (HEDM), namely, anaerobic baffled reactor (ABR) effluent and nitrified urine concentrate (NUC) and commercial hydroponic fertilizer mix (CHFM) as main sources of nutrients. Simulated spectral settings of HyspIRI, Landsat 9 and Sentinel 2-MSI were resampled from spectrometric proximally sensed data. Discriminant analysis (DA) was applied in discriminating tomatoes grown under these different nutrient sources. Results showed that the simulated spectral settings of HyspIRI sensor better discriminate tomatoes grown under different fertilizer regimes when compared to Landsat 9 OLI and Sentinel-2 MSI spectral configurations. Using the DA algorithm, HyspIRI exhibited high overall accuracy (OA) of 0.99 and a kappa statistic of 0.99 whereas Landsat OLI and Sentinel-2 MSI exhibited OA of 0.94 and 0.95 and 0.79 and 0.85 kappa statistics, respectively. Simulated HyspIRI wavebands 710, 720, 690, 840, 1370 and 2110 nm, Sentinel 2-MSI bands 7 (783 nm), 6 (740 nm), 5 (705 nm) and 8a (865 nm) as well as Landsat bands 5 (865 nm), 6 (1610 nm), 7 (2200 nm) and 8 (590 nm), in order of importance, were selected as the most suitable bands for discriminating tomatoes grown under different fertilizer regimes. Overall, the performance of simulated HyspIRI, Landsat 9 OLI-2 and Sentinel-2 MSI spectral bands seem to bring new opportunities for crop monitoring.


Introduction
Food insecurity is a large and growing challenge in sub-Saharan Africa [1,2]. It is estimated that at least one out of four people are hungry and undernourished in sub-Saharan Africa. The World Bank estimates that in 2030, nearly 9 in 10 extremely poor people will be living in Sub-Saharan Africa [3]. This is exacerbated by droughts and soil nutrients deficiencies resulting from limited fertilizer applications [1,2]. This is in turn associated with high fertilizer and food prices, amongst other factors. According to FAO [4] the annual food inflation increased from 5% in 2014 to 6% in 2018 whereas in Europe it remained stable and declined in Latin America, Asia, and Oceania. Subsequently, the improvement of crop production which leads to food security has been amongst the principal priorities required to fulfil the goals of sustainable human development as well as the African Union's Agenda 2063 [5]. Furthermore, food demand is anticipated to triple in sub-Saharan Africa after the projected 2.5-fold increase in population increase [6]. Specifically, a 60% increase in agricultural and horticultural production will be required by the increasing population in the light of diminishing water and soil nutrient resources [6]. The major concern is the current dietary transition, which is in favour of vegetables such as tomatoes amongst other crops, which is projected to increase, especially in urban areas, while water and soil nutrients are in decline [6][7][8].
Tomato fruits have a critical dietary role in providing folate, vitamins A, C and E; as well as antioxidants (lycopene, beta-carotene, gamma-carotene); trace elements of flavonoids; phytosterols and water-soluble vitamins important for human health [7]. To circumvent the challenge of decreasing soil nutrients and increase the production of vegetables (tomatoes) within a small land area, efforts have been exerted towards improving soil fertility and reducing expenses associated with commercial hydroponic fertilizer mix (CHFM) through the use of anaerobic baffled reactor (ABR) effluents and nitrified urine concentrate (NUC) both as a source of soil nutrients and water [8,9]. Smith and Smith [9], for instance, noted that nitrogen recovered from wastewater supported a high increase in tomato (Solanum lycopersicum) plant canopy volume, flower and fruit production when compared to plants treated with commercial hydroponic fertilizer mix (CHFM) which contained N, P, K, Ca, Mg and Si. Al-Hamdan, Cruise et al. (2014) in Jordan noted that treatment of tomato crop using waste water facilitated and increased their fruit size by up to 2 cm in diameter, and weight up to 78.7 g in relation to those administered with potable water in their field experiment. However, the challenge that has been lurking in the agricultural sector is the lack of comprehensive spatial explicit frameworks as well as objective criteria for crop growth and productivity monitoring. Spatially explicit data is important in effectively and precisely managing production both in the field and greenhouses to meet the increasing demand of high quality and safe agricultural products such as tomatoes. Information on vegetable crop type, growth, productivity or health status was previously measured in situ or done through routine field surveys which are often time consuming and date lagged [10,11]. Despite the fact that these in situ methods obtained plausible levels of accuracy in characterizing crops, they lacked spatial representativeness. Consequently, there is need for spatial explicit techniques that can be operationally used not only to characterize the crop's areal extent, but also their physiognomies due to lack or excess of soil nutrients. This information can help deduce and understand crop quality, growth and productivity patterns, which are critical in ensuring food security and coming up with well-informed intervention mechanisms or management strategies where necessary.
Meanwhile, earth observation technologies offer spatially explicit non-destructive synoptic views, innovative and economically feasible timely spatial scale means of generating farm scale crop monitoring. Literature shows that crop/plant physiology and structure such as leaf area index, water content plant pigment content, canopy architecture and canopy density are associated with specific key spectral wavebands [12][13][14][15]. Subsequently the variations in these biochemical and physical characteristics of crop plants caused by various crop management practices (i.e., different fertilizer regimes) facilitates unique variations in their spectral finger prints. This makes remote sensing and earth observation facilities to be critical reservoirs of spatially and crop explicit information required in ensuring good quality of crop produce. Remotely sensed data is robust and very sensitive to subtle vegetation traits such as those induced by different water and nutrient regimes. Rajah et al. showed that hyperspectral remotely sensed data could discriminate common dry beans that were rain-fed from those that were irrigated. Lu, et al. [16] discriminated tomato crops that were infected with multi-diseases at different phenological stages using hyperspectral data. Their results exhibited a high overall classification accuracy of 100% in discriminating multi-diseases during the early, asymptomatic and late stages of leaves growth. Above all, the advent and advancement of earth observation facilities has unveiled opportunities for assessing previously unresolved crop quality, growth and productivity related questions linked to plant physiognomies such as that induced by different fertilizer regimes on plant spectral characteristics [11]. Despite hyperspectral data's trade-off between cost and accuracy, it remains the most accurate spatial data for monitoring crop growth and productivity. Unlike broadband satellite data, hyperspectral data has numerous contiguous spectral channels with a potential ability to detect and characterize subtle differences in plant traits such as those induced by different fertilizer application regimes in relation. For example,Česonienė, et al. [17] demonstrated that hyperspectral data could discriminate between conventionally and organically grown carrots in Lithuania to Jeffries-Matusita distances ranging between 1.98 and 2.00. They attributed their results to the ability of hyperspectral data to discern on the variations in the cell structure conditions and canopy structure of the carrots grown under different farming methods. Meanwhile, multispectral sensors like Landsat, Satellite Pour l'Observation de la Terre (SPOT), MODerate Resolution Imaging Spectroradiometer (MODIS) are characterized with broad bands making it difficult to discern subtle plant traits as they tend to mask out critical plant information.
Although hyperspectral sensors provide accurate datasets, a number of sensors have been or are being developed with improved sensing capabilities because of the exorbitant acquisition costs associated with it [18]. For example, the earth observation community recently witnessed the launching of Sentinel 2 multispectral imager (MSI) and Landsat 8 OLI etc. Sentinel 2 MSI has been the first freely available sensor with a set of spectral wavebands covering the red-edge section of the electromagnetic spectrum (B5 (705 nm), 6 (740 nm), and 7 (783 nm)) at a relatively fine spatial resolution of 20m. The sensor has a wide swath-width of 290 km, coupled with a high spatial resolution of 10 m as well as a five-day temporal resolution making it a better facility for crop mapping and monitoring. Both sensors (i.e., Sentinel 2 MSI and Landsat 8 OLI) have been tested in various environmental application areas with plausible findings and conclusions [19][20][21][22]. However, in some instances they have been reported to experience challenges, especially when applied at farm level monitoring. This has been attributed to the presence of broad wavebands which are perceived to be concealing most important information. As a result, now new sensors such as the proposed Landsat 9 OLI-2 (with improved noise-to-signal ratio), Environmental Mapping and Analysis Program (EnMAP) and Hyperspectral Infrared Imager (HyspIRI) are being developed. The National Aeronautics and Space Administration agency (NASA) is looking forward to launching the state-of-the-art HyspIRI instruments covering the visible and near-infrared section (Vis/NIR) as well as the thermal infrared (TIR). There is a need to evaluate their performance in discriminating subtle plant properties resulting from crop management practices in relation to the available broadband multispectral sensors. These sensors could help by providing a spatially explicit non-destructive method of characterizing crop quality, growth and productivity patterns which are required in the food industry for pricing as well as in agricultural production for ensuring food security.
The upcoming hyperspectral instruments have a potential to supply the much-needed spatially explicit, accurate, consistent information on vegetable crops. Both of these instruments will be spectrometric covering the spectral ranges of 420-2450 nm and 380-2510 nm at different sampling distances of 6.5 nm for EnMAP's VNIR and 10 nm for EnMAP's SWIR section as well as HyspIRI's VSWIR [23,24]. The swath width of HyspIRI will be 185 km at 30 and 60 m spatial resolutions whereas EnMAP have 30-km-wide coverage across-track at a ground-sampling unit of 30m. The temporal resolution of EnMAP will be 4 days at the equator whereas that of HyspIRI will be 5 days. The fine spectral, spatial and temporal resolutions of these sensors make them relatively more suitable for agriculture applications. The major advantage with such instruments is that they will avail quality data at relatively low costs for data scarce regions such as the sub-Saharan Africa where resources are limited. In this regard, there is need to compare the performance of these hyperspectral sensors to the recently launched freely available and forthcoming multispectral sensors (i.e., Sentinel 2 MSI and Landsat 9 OLI-2 with improve spectral settings) so as to ascertain their full potential.
Considering the fact that literature states that hyperspectral data is generally characterized by high collinearity and that there is no specific algorithm that is suitable for discriminating vegetation characteristics at different places and times, this study also examined the performance of discriminant analysis (DA) and partial least squares discriminant analysis (PLS-DA) in discriminating tomato plants grown under different fertilizer regimes. PLS-DA and DA have been widely used in discriminating plants with different characteristics using both hyperspectral and broadband sensors data [25][26][27]. These Algorithms have been widely used because they offer an opportunity to evaluate and interpret the minute spectral pattern variations in plants, especially those grown under different management regimes. PLS-DA and DA classification ensembles construct a distinctive spectrum that represents the spectral signatures of the plants samples while simplifying the discrimination process when compared with other methods such as K-nearest neighbours [28][29][30]. Despite the optimal performance of these methods in vegetation discrimination, the robustness and accuracy of these two is yet to be established [31] particularly in discriminating fine spectral variations of vegetable crops (i.e., tomatoes) induced by different fertilizer application regimes. Furthermore, hyperspectral data could be widely renowned for its exceptional performance in literature in relation to broad bands sensors but the performance of these forthcoming sensors remains undocumented. This study therefore, sought to assess the discriminative performance of HyspIRI, Landsat 9 OLI-2 and Sentinel 2 MSI in characterizing tomato (Solanum lycopersicum) crops grown under commercial hydroponic fertilizer mix, anaerobic baffled reactor effluent and nitrified urine concentrate as nutrient sources. HyspIRI was selected as a representative of hyperspectral sensors in this study due to the fact that both EnMAP and HyspIRI will have similar spectral and spatial characteristics and that HyspIRI is going to cover a larger spectral portion of the electromagnetic spectrum particularly in the VSWIR sensor in relation to that of EnMAP. In addition, this study also assessed the performance of DA and PLS-DA in discriminating tomato crops grown under different fertilizer regimes using HyspIRI, Landsat 9 OLI-2 and Sentinel 2 MSI spectral data simulated from hyperspectral data.

Experimental Set-Up
A pot experiment was conducted in a hydroponic system that was set up in a polyethylene tunnel located at Newlands-Mashu Research Station under eThekwini Municipality, Durban, South Africa (29 • 46 25.648" E 30 • 58 28.329" S). The hydroponic system was designed to run three nutrients streams namely, anaerobic baffle reactor (ABR) effluents, nitrified urine concentrate (NUC) and commercial hydroponic fertilizer mix (CHFM) as a control. Each hydroponic system consisted of 150 L tank and the nutrient solution for each system was enclosed in a 100 L container stacked on the ground at the foot of each system.
Six-weeks-old, seedlings of 'Monica', a determinate tomato cultivar purchased from a local nursery (Sunshine seedlings, Pietermaritzburg, South Africa) were transplanted to 30 cm polyethylene pots filled with pine sawdust as a growing medium. The nutrient solution for each nutrient source was supplied to the plants using a pressure pump (DAB Model K30/70M, DAB Pumps, MarcoPolo, Mestrino, Italy) via a 20 m irrigation line. A 20 cm drip irrigation emitters (2 L) were placed and irrigation was performed at six intervals of 5 min/duration daily using a timer. The study was arranged using a complete randomised design with three replications of five plants each, giving fifteen experimental units per nutrient source.
Tomato plants of the control treatment were irrigated with a commercial hydroponic fertilizer mix (Hygroponic ® and Solu-cal ® ) at the rate of 800 g + 620 g/1000 L of water as recommended for hydroponic tomato production; NUC, commercial fertilizer application rate was used as a standard as recommended by Jonsson et al. (2004) and ABR effluents with no specified application rate. For the CHFM and NUC treatment, the fertilizer was mixed using municipal tap water whereas for the ABR treatment only effluent from the anaerobic baffle reactor component was used as nutrient and irrigation source. The experiment was allowed to run for 12 weeks before the crop was harvested.
Remotely sensed proximal spectra data was collected at vegetative stage (i.e., four weeks after transplanting) and two weeks after flowering. In this study all requirements for the photosynthetic activity of the tomato plants were the same except for the fertilisation treatments. This was done so as to make sure that any spectral differences between the CHFM, ABR, NUC treated tomato plants were mainly due to the fertiliser applications since all conditions were the same.

Remotely Sensed Data
The FieldSpec-3 (ASD Inc., Boulder, CO, USA) Analytic Spectral Device (ASD) FieldSpec instrument was used to acquire the spectral reflectance of tomatoes plants receiving ABR effluents, NUC and CHFM. The ASD measured the radiation at 1.4 nm intervals for the 350-1000 nm and 2 nm intervals for the 1000-2500 nm spectral regions. The reflectance measurements were conducted using bare fibre-optic held at nadir position~0.5m above the tomato canopies resulting in a field of view with a diameter of 0.225m. This diameter was found to be adequate to capture the reflectance of the tomato canopies. The normalization of spectral measurements was conducted after every 5 to 10 spectral measurements, using a standard white spectralon. A standard spectralon is optically flat to ±4% over the range of 250-2500 nm and ±1% over the photopic region of the spectrum [32,33]. Its diffuse reflectance standards are highly Lambertian [33]. This was done to circumvent the possible changes in weather conditions as well as irradiance from the sun [34]. The spectral measurements were conducted under clear skies during the day between 10:00 and 14:00 since this is the time with maximum net radiation. A total of 900 spectral samples were measured on canopies of tomatoes treated with CHFM (n = 300), ABR (n = 300) and NUC (n = 300). Noise regions between 350 and 399, 1355 and 1420 nm, 1810 and 1940 nm and 2470-2500 nm were cleaned prior to any analysis of the data [35]. The spectrometrically sensed proximal spectra were then resampled into the spectral settings of HyspIRI, Sentinel 2_MSI and Landsat 9 OLI-2 (  [37]. To simulate the ASD data into the spectral settings of broad spectral sensors, the PRISM in IDL-ENVI was used.

Discriminating Tomato Plants Grown under Different Fertilizer Regimes
Exploratory data analysis was conducted to determine if the data followed a normal distribution curve. Normality test was performed using the Kolmogorov-Smirnov test. Further, we assessed spectral separability as well as administered a pre-filter [38] after hyperspectral data resampling it to HyspIRI, Sentinel-2 MSI and Landsat 8/9 OLI spectral configurations. Resampling was based on the Analysis of Variance test (ANOVA). We then conducted post hoc test to establish the channels that exhibited significant differences between the spectral data of the tomato crops receiving nutrients from ABR effluent, NUC and CHFM.
The other objective of this study was to assess the accuracies of partial least squares discriminant analysis and discriminant analysis algorithms in characterizing tomato crops grown under the three fertilizer regimes. In that regard, we used the PLSDA and the DA to classify the spectral reflectance of tomato crops growing under UNC, UNF and ABR fertilizer treatments. Details about DA and PLS DA are provided in Zhang, et al. [39] and Boulesteix [29]. Prior to conducting PLS-DA and DA, the spectral samples were partitioned into training (70%) and testing (30%) data.

Classification Accuracy Assessment
The 30 percent of the samples were used for classification accuracy assessment (90 per treatment). Confusion matrices were computed and used to evaluate the classification accuracies of the DA and PLS DA models. We further computed the overall, producer and user accuracies, as well as the kappa statistics for each set of spectral settings as classified by the two algorithms based on the confusion matrices. To compare the performance of the two algorithms, a McNemar's test was conducted as detailed by Manandhar, et al. [40] and de Leeuw, et al. [41].

Results
Analysis of variance tests results showed significant differences (α = 0.05) between tomato plants treated with different fertilizer combinations based on the spectral settings of HyspIRI, Sentinel-2MSI and Landsat 9 OLI. Figure 1 illustrates spectral signatures (mean spectral reflectance) of tomato plants grown under ABR, CHFM and NUC treatments. Specifically, the spectral signatures of tomatoes grown under NUC exhibited higher reflectance curves when compared to those growing under CHFM and ABR in all the sections of the electromagnetic spectrum across all the sensors. The CHFM spectral signature was the lowest in comparison to the other to fertilizer treatments across all the spectral signatures and sensors. For HyspIRI resampled data, significant differences were observed in the visible, NIR as well as the SWIR portions of the electromagnetic spectrum ( Figure 1). Meanwhile the most glaring differences in the reflectance of tomato plants grown under ABR, CHFM and NUC treatments were observed in the NIR portion of the electromagnetic spectrum based on the simulated Sentinel 2 MSI and Landsat 9 OLI-2 data (Figure 1). HyspIRI spectral settings exhibited more potential spectral windows of separability between tomato plants grown under ABR, CHFM and NUC treatments when compared with the broadband sensors. Potential spectral windows of separability exhibited by Sentinel 2 MSI and Landsat 9 OLI-2 were only in the near infrared regions (Figure 1).

Discriminating Tomatoes under ABR Effluent, NUC and CHFM
All sensors better characterized tomato crops administered with ABR when compared to those administered with CHFM and NUC. Specifically, high producer and user accuracies ranging between 91% and 100% were observed in characterizing tomatoes treated with ABR (Table 2). Meanwhile tomatoes administered with ABR and CHFM were characterized with slightly lower accuracies (Table 2). Moreover, HyspIRI, produced high accuracies characterized by kappa statistics of 0.99, whereas the spectral settings of Sentinel-2 MSI and Landsat 9 OLI-2′s spectral settings exhibited kappa statistics of 0.85 and 0.79, respectively. HyspIRI spectral settings exhibited high producer accuracies of 100%, 92% and 100% for ABR, CHFM and NUC respectively. Meanwhile Sentinel-2 MSI exhibited slightly lower producer accuracies (ABR = 91%, CHFM = 86% and NUC = 100%). Landsat OLI-2 exhibited producer accuracies that were comparable to those of Sentinel 2 MSI which were 91% for ABR, 83% for CHFM and 100% for NUC. The same trend could be observed on the user accuracies ( Table 2). In general, the CHFM treatments had slightly lower-class accuracies across all sensors in relation to the HEDM treatments.

Performance of DA and PLS-DA Algorithms in Discriminating Tomatoes under ABR Effluent, NUC and CHFM
When comparing the performance of algorithms, DA exhibited very high accuracies. For instance, DA's producer accuracies derived using DA across all fertilizer treatments and sensor

Discriminating Tomatoes under ABR Effluent, NUC and CHFM
All sensors better characterized tomato crops administered with ABR when compared to those administered with CHFM and NUC. Specifically, high producer and user accuracies ranging between 91% and 100% were observed in characterizing tomatoes treated with ABR (Table 2). Meanwhile tomatoes administered with ABR and CHFM were characterized with slightly lower accuracies (Table 2). Moreover, HyspIRI, produced high accuracies characterized by kappa statistics of 0.99, whereas the spectral settings of Sentinel-2 MSI and Landsat 9 OLI-2 s spectral settings exhibited kappa statistics of 0.85 and 0.79, respectively. HyspIRI spectral settings exhibited high producer accuracies of 100%, 92% and 100% for ABR, CHFM and NUC respectively. Meanwhile Sentinel-2 MSI exhibited slightly lower producer accuracies (ABR = 91%, CHFM = 86% and NUC = 100%). Landsat OLI-2 exhibited producer accuracies that were comparable to those of Sentinel 2 MSI which were 91% for ABR, 83% for CHFM and 100% for NUC. The same trend could be observed on the user accuracies (Table 2). In general, the CHFM treatments had slightly lower-class accuracies across all sensors in relation to the HEDM treatments. When comparing the performance of algorithms, DA exhibited very high accuracies. For instance, DA's producer accuracies derived using DA across all fertilizer treatments and sensor simulations ranged from a minimum of 89% whereas PLS-DA had a slightly lower minimum of 83% ( Table 2). The user accuracies derived using DA ranged from a minimum of 76% whereas those derived using PLS-DA ranged from 65% to 100%. The overall accuracies derived using DA were higher (i.e., 0.94-0.99) when compared to those derived using PLS-DA (i.e., 0.3-0.97). Uniformly, the kappa statistics derived using DA were higher, ranging between 0.79 and 0.99 whereas those derived using PLS-DA were lower ranging between 0.63 and 0.90 (Table 2). Table 3 summarizes potential spectral variables for discriminating tomato plants grown under ABR, CHFM and NUC treatments using HyspIRI, Sentinel 2 MSI and Landsat 9 OLI-2. HyspIRI wavebands 710, 720, 690, 840, 1370 and 2110 nm, Sentinel 2-MSI bands 7, 6, 5 and 8a as well as Landsat bands 5, 6, 7 and 8, in order of importance, were identified as the most optimal spectral variables for discriminating tomato plants grown under different fertilizer regimes. NIR was the most prominent optimal section with a potential of discriminating tomatoes plants grown under different fertilizer regimes across the three sensors, although it did not outperform the red-edge section of the electromagnetic spectrum which was the most influential region.

Discussion
We sought to compare the strength of HyspIRI's spectral configuration in relation to Landsat 9 OLI-2 and Sentinel 2 MSI spectral settings in characterizing tomato (Solanum lycopersicum) crops grown under CHFM, ABR and NUC treatment regimes. Results of this study showed that tomatoes that were fertilized using ABR could be optimally discriminated (i.e., Kappa statics ranging 0.79 to 0.99) from those that were administered with CHFM and NUC. The spectral reflectance of HDME fertilized crops was higher than that of those treated with CHFM. This indicates the potential of HEDM fertilizer for sustaining the health of vegetables in a manner comparable to that of conventional chemical fertilizers. This could be attributed to the fact that ABR effluents have nutrient properties favorable for tomato plants, which facilitate excessive vegetative growth, high biomass accumulation, delayed or uneven maturity [42][43][44]. Tomato crops with excessive vegetative growth, high biomass accumulation, delayed or uneven maturity tend to be easily detected and discriminated by satellite sensors compared to those which are not. Subsequently, tomatoes growing under the ABR treatments have a different spectral signature from those that are fertilized using CHFM and NUC. Literature illustrates that ABR tends to facilitate high biomass accumulation (i.e., increased leaf area index) hence the high classification accuracies exhibited by all remotely sensed data in this study in characterizing the HEDM fertiliser treatments [11,28,29]. For example, Al-Lahham, El Assi and Fayyad [45] illustrated that tomato crops that were administered with high quantities of waste water had big fruit sizes hence high biomass accumulation in relation to those that were administered with potable water. In a related study, Zavadil (2009) noted that primary treated waste water which contained an average of 14-fold nitrogen amounts (70.6 mg/L, which was 89% ammonia form) and also a 3-fold of the total phosphorus, resulted in high biomass accumulation and yields of lettuce salad, radishes and carrots in their study assessing the influence of sewage in relation to potable water on vegetable growth. Subsequently, the increases in biomass accumulation associated with wastewater treated vegetables could explain the discrimination of ABR treated tomato crops in comparison to those administered with CHFM and NUC, in this study. Furthermore, the most influential wavebands that facilitated high classification accuracies in discriminating tomato plants grown under human-excreta derived materials (HEDM) from those growing under CHFM conditions were from the red-edge section of the electromagnetic spectrum which is generally associated with healthy plants. The red-edge section of the electromagnetic spectrum is very sensitive to high chlorophyll, leaf angle distribution, leaf area index levels associated with plants grown in more nutritious environments such as the those exhibited by ABR. As aforementioned, HEDM have a high turnover of N P K plant nutrients which is comparable to that of commercial soluble fertilizers [46][47][48]. In this regard, the high turnover of plant nutrients such as nitrogen facilitates an increase in biomass, LAD, LAI which in turn makes the canopy spectral signature of those plants treated with HEDM to be discriminable from other plants especially in the red-edge and the NIR sections of the electromagnetic spectrum [49,50]. Subsequently, the red-edge and the NIR regions of the electromagnetic spectrum illustrate that the spectral signatures and the health of plants grown under HEDM are comparable to those grown under the CHFM.
When assessing the performance of sensors, HyspIRI outperformed the two multispectral sensors namely, Sentinel-2 MSI and Landsat OLI spectral settings in discriminating tomato crops grown under different fertilizer regimes. This could be explained by the fact that HyspIRI is a hyperspectral sensor characterized by narrow spectral wavebands that are more sensitive to the spectral reflectance of tomato crops grown under different fertilizer regimes than broadband sensor settings such as those of Landsat which could be masking out those minute tomato crops spectral variations. There is a consistently growing body of literature that supports the claim that hyperspectral sensors are more sensitive to minute vegetation spectral variabilities compared to broadband sensors due to the narrowed bandwidths configuration [51][52][53][54][55][56]. Specifically, Thenkabail, Smith and De Pauw [55,57] illustrated that narrow bands characterized different crop traits such as yield as well as spectral variations when compared to broadband spectral data. Also, they were able to better characterize wheat from barley using hyperspectral data in relation to the data from broadband sensors. They attributed this to the variation in spectral settings (bandwidths) of the sensors they used. These variations in spectral settings affected the detail that determined the accuracy of their models for plant trait characterization. Meanwhile, Clark [57] noted that there was no significant variation in the performance of HyspIRI and Sentinel-2 MSI as well as Landsat OLI in landcover classification of the San Francisco Bay Area in northern California, USA. However, their results confirmed that HyspIRI exhibited higher classification accuracies in their study.
Results of this study also illustrated that Sentinel 2 MSI and Landsat 9 performed satisfactorily in discriminating tomatoes grown under different fertilizer regimes, although Sentinel-2 MSI outperformed Landsat OLI. This could be explained by the fact that Sentinel-2 MSI spectral settings cover the red-edge portion of the electromagnetic spectrum which is critical in mapping and detecting various vegetation traits. Moreover, there is a large and growing body of literature illustrating that Sentinel-2 MSI performs better than Landsat OLI in vegetation mapping [22,57,58]. The study by Colkesen and Kavzoglu [58] illustrated that Sentinel-2 MSI outperformed Landsat OLI in discriminating alfalfa, sugar beet and bean in the agricultural lands of the Ferizli district, Turkey. Shoko and Mutanga [22] also illustrated the robustness of Sentinel-2 MSI remotely sensed data in better discriminating C3 Festuca costata Nees from the C4 Themeda triandra Forssk grasses in a mountainous area in South Africa. They also attributed the optimal performance of Sentinel-2 MSI to the presence of red-edge bands in discriminating Festuca costata from Themeda triandra grasses.
Even though this study successfully illustrated the robustness of the spectral settings of these sensors in discriminating tomato crops grown under different fertilizer regimes using simulated data, there is still need to assess their satellite remotely sensed data. This study only sought to evaluate the spectral configuration of these sensors, hence there is still need to evaluate their radiometric and spatial resolutions in a similar application setup. Generally, the reflectance of remotely sensed satellite data tends to be affected by numerous factors associated with the sensor's platform, atmospheric related influences etc. [59][60][61]. Factors that affect satellite remotely sensed data are different from those that affect spectrometric proximally sensed data. Spectrometric proximally sensed data tends to be less affected by atmospheric related influences; hence, it offers a robust dataset suitable for testing the spectral settings of sensors, particularly the forth coming ones. The performance of the same sensors' spectral settings derived using satellite platforms, therefore, still needs to be evaluated.
Although this was not the major objective of the study, DA outperformed PLS-DA in discriminating tomato crops grown under different fertilizer regimes. In this study PLS-DA failed to derive unnecessary variables for characterizing tomato crops grown under different fertilizer regimes. On the other hand, there are numerous studies that have illustrated the optimal performance of DA in dimension reduction as well as feature extraction [62][63][64]. Our results also indicated that DA selected more spectral bands in relation to PLS-DA as optimal spectral variables for discriminating tomatoes grown under different fertilizer treatment regimes. The optimal performance of DA in relation to PLS-DA could be attributed to the conservative nature of PLS-DA in classifying tomato plants [31].

Implications of the Study'S Findings for Horticultural Crop Production
The optimal performance of HyspIRI's spectral settings in characterizing tomato plants grown under different fertilizer treatments illustrates its great potential in providing additional invaluable spatially explicit information urgently required in horticulture management and decision-making processes particularly for successful site and crop specific management practices. For instance, this study shows that based on remotely sensed data HEDM fertilizers have a great potential of improving vegetable crops health in manner comparable to that of conventional fertilisers. The fact that this study showed that HyspIRI, Sentinel 2 MSI and Landsat 9 OLI-2 could detect the differences between HEDM fertilized plants and those that are grown under conventional fertilizers underscores the potential of these sensors in providing information that could also be used in other applications such as crop inventory, condition, production forecast, assessment of nutrient deficiencies as well as growth and health of tomato plants and other horticultural crops. Furthermore, these sensors offer cheap, fast and reliable spatial explicit information on crop conditions required in the administration of precise and proper fertilizer applications while optimizing resources and increasing net returns especially in sub-Saharan countries with limited data and financial resources for facilitating high production of horticultural crops. However, more and extensive research efforts still need to be exerted in fully ascertaining and exploiting the potential of remotely sensed data in other horticultural crop-management applications.

Conclusions
The prime objective of this study was to compare the strength of the forthcoming hyperspectral sensor HyspIRI's spectral settings in the context of characterizing the effects of different fertilizer-treatment regimes on tomato crops. Furthermore, the study assessed the performance of PLS-DA in relation to DA in discriminating tomatoes treated with ABR, NUC and CHFM. Based on the results exhibited by this study we conclude that: • The forthcoming HyspIRI sensor has the potential to accurately map tomato crops under various fertilizer regimes. Landsat and Sentinel, performed comparably to HyspIRI spectral settings.

•
Overall, all the sensors were able to characterise the comparable impact of HEDM in relation to that of CHFM fertilizers on the spectral characteristics of tomato plants as a proxy of their health. • DA offers optimal accuracies in characterizing tomatoes grown under different fertilizer regimes when compared to PLS-DA.
These findings are a substantial foundation upon which comprehensive precision agricultural assessments initiatives could be formed. These initiatives are required in order to attain sustainable agriculture as well as food security in regions such as sub-Saharan Africa where agricultural crop monitoring is currently hindered by the limited access to robust spatial data sets.