Monitoring of Khorasan ( Triticum turgidum ssp. Turanicum ) and Modern Kabot Spring Wheat ( Triticum aestivum ) Varieties by UAV and Sensor Technologies under Different Soil Tillage

: Khorasan wheat ( Triticum turgidum ssp. turanicum ( Jakubz. )) is an ancient tetraploid spring wheat variety originating from northeast parts of Central Asia. This variety can serve as a full-ﬂedged alternative to modern wheat but has a lower yield than modern varieties. It is commonly known that wheat growth is inﬂuenced by soil tillage technology (among other things). However, it is not known how soil tillage technology affects ancient varieties. Therefore, the main objective of this study was to evaluate the inﬂuence of different soil tillage technologies on the growth of the ancient Khorasan wheat variety in comparison to the modern Kabot spring wheat ( Triticum aestivum ) variety. The trial was arranged in six small plots, one half of which was sown by the Khorasan wheat variety and the other half of which was sown by the Kabot wheat variety. Three soil tillage methods were used for each cultivar: conventional tillage (CT) (20–25 cm), minimum tillage (MTC) with a coulter cultivator (15 cm), and minimization tillage (MTD) with a disc cultivator (12 cm). The soil surface of all of the variants were leveled after tillage (harrows & levelling bars). An unmanned aerial vehicle with multispectral and thermal cameras was used to monitor growth during the vegetation season. The ﬂight missions were supplemented by measurements using the GreenSeeker hand-held sensor and plant and soil analysis. The results showed that the Khorasan ancient wheat was better suited the conditions of conventional tillage, with low values of bulk density and highvalues of total soil porosity, which generally increased the nutritional value of the yield in this experimental plot. At the same time, it was found that this ancient wheat does not deplete the soil. The results also showed that the trend of developmental growing curves derived from different sensors was very similar regardless of measurement method. The sensors used in this study can be good indicators of micronutrient content in the plant as well as in the grains. A low-cost RGB camera can provide relevant results, especially in cases where equipment that is more accurate is not available.


Introduction
Tillage is an important process that affects soil properties such as the temperature, moisture content, bulk density, or stability of soil aggregates [1]. The method of processing and habitat conditions give soil a specific value that determines the growth, development, and yield of crops [2]. Conventional soil tillage, which is still the most commonly used method, is one of the most demanding tilling operations, as its depth is most often between 20 and 25 cm [3]. The disadvantage of this tillage method is that it can adversely affect the structure of the soil due to the increased disturbance of soil aggregates [4]. Conventional tillage with moldboard ploughing (25-28 cm depth) followed by spring tine harrowing compared to reduced tillage (two harrowings with spring tine harrows up to 4-6 cm depth and sowing) has a worsening effect on the soil. This means reduced organic matter, poorer soil structure, reduced aggregate content, and increased bulk density [5]. Bilalis et al. [6] evaluated a field experiment comparing the influence of three different tillage methods (conventional, minimization, and no-till) on the quality, yield, and growth of four varieties of wheat ('Siette', 'Panifor', 'Myrto', and 'Estero'). The minimum tillage and no-till methods had the best effects on the soil (porosity and total nitrogen) while the highest grain yield was shown through the use of conventional tillage (CT) with the 'Siette', 'Myrto', and 'Estero' varieties, while the 'Panifor' variety had the highest yield after the use of minimization tillage. Grandy et al. [7] compared ploughing tillage to a depth of 30 cm with a notill system. If the strategies and conditions for the habitat are chosen correctly, no-till can achieve the same results or even better results than conventional tillage. In the EU, ecological systems are becoming increasingly popular than conventional ones, with the share of ecologically managed land amounting to around 4% [8]. In the Czech Republic, conventional tillage predominates, mostly in areas rich in precipitation. In 2016, 66.5% of arable land was cultivated with conventional technology, 32.1% with conservation technology, and 1.4% with zero tillage [9]. The influence of different tillage on winter wheat yield at three locations in the Czech Republic during 6 seasons showed that the reduced tillage system combined with a high input level achieved higher yields for all wheat varieties than conventional tillage [10].
Wheat is among one of the biggest strategic agriculture crops around the world. Wheat has been the staple food of the population for more than 8000 years [11]. With a growing population, the demand for wheat will increase by up to 70% by 2050. This will mainly be due to a change in the eating habits of a growing population and climate change [12]. The original civilizations in Europe, West Asia, and Africa grew ancient wheat, which is now terminologically called einkorn, emmer, khorasan, and spelt. These are varieties bred from either wild or ancient wheat. Recently, the popularity of these natural varieties has been growing, mainly due to the allegedly better tolerance of people who are allergic to modern wheat [13,14] and protection against the development of type-2 diabetes [15]. Khorasan wheat (Triticum turgidum ssp. turanicum) has a registered trademark KAMUT ® . The KAMUT ® trademark must meet given specifications, especially regarding the nutritional properties of the grain and its growing conditions, which are focused on organic cultivation [16].
Differences in plant growth, which are also sometimes due to the response to the tillage technology used, can be assessed by modern monitoring systems. Platforms for crop monitoring provide either remote or proximal sensing. Remote sensing platforms can be distinguished according to the operative height from which the crop canopy is scanned. Optical remote sensing includes traditional satellite systems for monitoring natural resources such as Landsat or Sentinel 2. Using satellite images, it is only possible to yield predictions on larger plots [17][18][19]. However, these satellite images are still insufficient for objectives such as breeding experiments or plant health monitoring, where there is a very high demand for spatial accuracy.
With the technical progress of the past two decades, it has become possible to obtain images from unmanned aerial vehicles, which have a very high spatial resolution and can thus monitor the crop canopy with very high accuracy [20]. Another advantage in comparison to satellite images is the flexibility of the acquisition time and the ability to focus on selected phenological phases of the monitored stands. Another advantage of modern unmanned aerial vehicles (UAVs) is the ability to acquire images with several sensors simultaneously, for example, thermal images are supplemented by multispectral images [21].
Proximal sensors such as the GreenSeeker handheld crop sensor can be used for nitrogen uptake monitoring, yield prediction, and thus for determining site-specific nitrogen management strategies in wheat [22]. The GreenSeeker handheld crop sensor is not suitable for measuring the normalized difference vegetation index (NDVI) status of medium and large agricultural fields due to the demands of data collection and processing to obtain the Agronomy 2021, 11, 1348 3 of 16 resulting map [23]. On the other hand, Ali et al. [22] found that GreenSeeker can give better information regarding crop canopy than contact sensors, which measure greenness on one spot on a leaf. Proximal sensors can find a better application in the regular monitoring of smaller agricultural plots, as in the case of the Naser et al. [24] study.
Sensors mounted on the platforms described above can provide an image in a minimum of one spectral band up to several spectral bands in various wavelength ranges. A multispectral camera is usually used for crop status monitoring and can be used to determine aspects such as structure, vigor, variability, or the chlorophyll content in leaves [25][26][27][28][29]. Thermal images can provide information on the condition of crops damaged by drought [30]. The spectral resolution of a thermal camera normally ranges from 7 to 14 µm in one spectral band [31].
Crop canopy status is usually expressed by vegetation indices as the ratio of difference in the sums of the sensor measurements, usually in two or more bands based on the strong reflectance in the near infrared part and strong absorption in red part of the electromagnetic spectrum [32]. The NDVI is one of the longest used vegetation indices, with its origins dating back to 1974 [33]. This index is used as a reliable tool to determine and evaluate the structure, vitality, health, and other biophysical parameters of agricultural crops [34][35][36]. The main advantage of this index is its easy readability and orientation in the determined values. In general, values lower than 0.25 signal bare soil. Higher values indicate varying degrees of vitality and crop canopy coverage, which is important, for example, in the tillering growth stage [37]. The NDVI alternative for the estimation of chlorophyll in leaves is the green normalized vegetation index (GNDVI) developed by Gitelson et al. [38], which uses the near-infrared and green electromagnetic parts of the spectrum in various ratios for chlorophyll content estimation. The usability of the GNDVI index was proven using corn, garlic, oats, onions, potatoes, sunflowers, alfalfa, and grapevines [39]. The triangular greenness index (TGI) spectral index [40] can fully serve the same purposes as NDVI or GNDVI if only images in the visible part of the spectrum are available [41]. Some authors, Gitelson & Merzlyak [42] or Xie et al. [43], for example, found through their research that vegetation indices with the red-edge part of the electromagnetic spectrum can be sensitive to leaf chlorophyll, as can the traditional VIs. The advantage of the red edge bands is using the relatively narrow part of the electromagnetic spectrum that relates to a position between the red absorption maximum and near infra-red high reflectance. Xie et al. [43] mentioned that the chlorophyll index green (CIG; [44]) has a close relationship to both chlorophyll content and leaf area index (LAI; [45]), but the chlorophyll index red edge (CIR; [46]) is more sensitive to small changes in the canopy, gap fraction, and senescence [47]. Using red-edge bands, it is possible to capture its shift to longer wavelengths, which is usually caused by increasing chlorophyll content [48]. The CIG and CIR were used by the authors of Xie et al. [43] with the aim of analyzing how changes in chlorophyll content and LAI contribute to reflectance variability in the RED, NIR, and red-edge bands. They monitored the spectral properties of winter wheat, barley, alfalfa, and maize.
It is clear from the literature review that the influence of soil tillage technology to wheat growth was studied in the past. Nevertheless, there is no information from literature regarding this influence on the ancient Khorasan variety. Therefore, the main objectives of this study were to evaluate how different types of soil tillage technology can affect the growth and conditions of old Khorasan wheat compared to modern Kabot spring wheat and to assess the suitability of the methods used to measure the condition of the stand during the vegetation season.

Experimental Area and Soil Tillage
The study was conducted in the 2020 vegetation season and was located in six regular plots in the size of 4 × 50 m nearČenovice (49 • 47 42.78 N, 15 • 6 35.94 E), a central Bohemian region in the Czech Republic. The experimental area has a south-eastern aspect with a 2.63 • slope and is 481 m a.s.l. on average. The soil is Haplic Luvisol. The soil classification (sampled from a depth of 5-10 cm) according to WRB (The World Reference Base) in percentages was clay 15.1%; silt 30.9%; very fine sand 9.2%; fine sand 14.3%; medium sand 12.1%; coarse sand 9.4%; very coarse sand 9%. Soil classification was performed on a Horiba LA-960 device (dry dispersion). Temperature and precipitation were measured by an agrometeorological station located next to the experiment. Table  1 lists the precipitation rates and temperatures for the observed 2020 cropping season for wheat varieties. Meteorological data are divided according to the main phenological stages on the BBCH scale [49]. Temperature values for the dates of UAV scanning are also given there. The experimental area was divided into six plots. The size of one plot was 4 × 50 m (200 m 2 ) with a 4 m wide handling gap between the parcels. Khorasan was sown in one-half of the trial area (plots 1-2-3) and the Kabot variety was sown in the other half (plots 4-5-6) of the experimental field. Three different soil tillage technologies were used for each section of the cultivar ( Figure 1). Plots 1 and 4 were cultivated using conventional tillage (CT) with moldboard plowing at a depth of 20-25 cm. Plots 2 and 5 were cultivated by means of minimum tillage (MTC) using a coulter cultivator twice (at a depth of 15 cm). Plots 3 and 6 were cultivated using minimization tillage (MTD) with the use of a disc cultivator twice (at a depth of 12 cm). The soil surface for all variants was leveled with a harrows (4 cm) and leveling bars.
Sowing took place on 12.4.2020. After the levelling the soil surface using harrows, the NPK fertilizer was manually scattered on the individual plots at a rate of 200 kg·ha −1 (4 kg per parcel). Further fertilization using NPK took place after the first measurement on 18.5 on 7.5. Mustang was applied against the developed dicotyledonous weeds. LAV fertilizer at a dose of 200 kg·ha −1 was applied manually on 2.6. The harvest took place on 21 August 2020 and was completed using the FORTSCHRITT E516B combine harvester. The harvest took place on the individual plots. The combine harvester hopper was emptied onto an unfolded tarpaulin, and the contents of each plot were separated and weighed separately. Sowing took place on 12.4.2020. After the levelling the soil surface using harrows, the NPK fertilizer was manually scattered on the individual plots at a rate of 200 kg•ha −1 (4 kg per parcel). Further fertilization using NPK took place after the first measurement on 18.5 on 7.5. Mustang was applied against the developed dicotyledonous weeds. LAV fertilizer at a dose of 200 kg•ha −1 was applied manually on 2.6. The harvest took place on 21 August 2020 and was completed using the FORTSCHRITT E516B combine harvester. The harvest took place on the individual plots. The combine harvester hopper was emptied onto an unfolded tarpaulin, and the contents of each plot were separated and weighed separately.

UAV Canopy Sensing and Data Processing
Four flight were performed during the cropping season (see Table 1) using a fixedwing eBeeX drone with a MicaSense RedEdge MX camera (AgEagle Aerial Systems Inc., Wichita, KN, USA) and a DuetT dual camera (senseFly SA, Route de Genève 38, Cheseaux-sur-Lausanne, Switzerland) comprising of a 20 Mpx RGB sensor with a 28 mm focal lens-S.O.D.A. camera and an infrared camera developed by FLIR technology for thermal mapping (details in Table 2). Table 2. Center wavelengths for individual spectral bands and bandwidth of sensors used in this study.

UAV Canopy Sensing and Data Processing
Four flight were performed during the cropping season (see Table 1) using a fixedwing eBeeX drone with a MicaSense RedEdge MX camera (AgEagle Aerial Systems Inc., Wichita, KN, USA) and a DuetT dual camera (senseFly SA, Route de Genève 38, Cheseauxsur-Lausanne, Switzerland) comprising of a 20 Mpx RGB sensor with a 28 mm focal lens-S.O.D.A. camera and an infrared camera developed by FLIR technology for thermal mapping (details in Table 2).  (above elevation data) with a spatial resolution of 6.0 cm·px −1 for the resulting images. The MicaSense Red Edge MX camera required radiometric calibration before and after each flight using a calibrated reflectance panel (CRP). The CRP made it possible to convert the raw pixel values from the RedEdge images to absolute reflectance. The information regarding the calibration was stored in the EXIF metadata and was available for further image processing. Flights with the DuetT camera were performed with 80% longitudinal and 75% lateral overlap and at 91.8 m above the elevation data flight altitude, with a 12.0 cm·px −1 resulting spatial resolution for the thermal images. The S.O.D.A. camera had the following settings: 84% longitudinal and 83% lateral overlap and a resulting 2.1 cm·px −1 spatial resolution for the RGB images. Accuracy was ensured by VRS.MAX-CZEPOS (master auxiliary stations, RTCM 3.1. correction format) provided by the CZEPOS service from the Czech Office for Surveying, Mapping and Cadastre of the Czech Republic. The resulting average accuracy of these missions was 4.8 cm horizontally and 5.2 cm vertically.
The images were pre-processed and refined in the postflight tool eMotion SW were uploaded into Pix4D SW (Pix4D S.A., Route de Renens 24, Prilly, Switzerland) and processed by a common procedure that consisted of several steps that ensured high accuracy of the resulting orthomosaic (details available at: https://support.pix4d.com/hc/en-us/ articles/202557759-Menu-Process-Processing-Options, accessed on 15 May 2021). Images from Duet T thermal camera were processed using Duet T templates. The thermal maps were visualized in absolute temperature values. The resulting layers were in the GeoTIFF format georeferenced in the WGS84 UTM Zone 33N coordinate system.
The NDVI, GNDVI, CIR, and TGI (for images from S.O.D.A. and MicaSense camera) spectral indices and plant height, which derived from digital terrain model (DTM) and the digital surface model (DSM), were calculated for each of the data sets (details in Table 3). Plant height was computed from the RGB orthophoto as the difference between the DSM derived for each scanning term and the first scan date of the DTM, which was completed when the soil was best detectable.

Proximal Canopy Sensing, Soil and Plant Sampling
The GreenSeeker (version 1.00, Rev B, 2012, Trimble Inc., Sunnyvale, CA, USA) was used to obtain the NDVI values directly above the selected areas of plant cover. The NDVI values were able to convert the information from the RED and NIR reflectance immediately into values of the NDVI spectral index [50]. Data were collected at 72 sampling points (see Figure 1) during the vegetation season for wheat varieties on 18 May, 11 June, 14 July, and 24 July.
Soil samples were put into physical cylinders in the form of undisturbed soil samples (Kopecky's cylinder) with a capacity of 100 cm 3 from a depth of 5-10 cm. The evaluation was completed according to the Novak method [51]. The following characteristics were Agronomy 2021, 11, 1348 7 of 16 evaluated from the samples: initial moisture content, bulk density, total porosity, capillary, and non-capillary porosity. Samples were taken three days after the wheat harvest.
Other soil samples were taken from a depth of 10-15 cm from each parcel after harvesting the grain on 21 August 2020. The samples were given to an external laboratory at the Crop Research Institute in Prague, where the necessary analyses were carried out. The samples identified the levels of accessible nutrients (P, K, Mg, Ca)-determined in the Mehlich III extractant, as well as the microelement content (B, Cu, Mn, Zn, Fe)-also specified in the Mehlich III extractant [52]. This method is the most widely used method to determine the content of nutrients or risk elements in the soil. Methodology: The soil sample (weigh 10 g) was shaken on a rotary shaker for 10 min into a sealable plastic container with a capacity (200-400) mL with the Mehlich III extractant (modification by Trávník et al. [53]). The Mehlich III extractant was composed of the following: 0.2 mol·L −1 of concentrated acetic acid CH 3 COOH, 0.015 mol·L −1 of ammonium fluoride solution NH 4 F, 0.013 mol·L −1 HNO 3 , 0.25 mol·L −1 of ammonium nitrate NH 4 NO 3 , and 0.001 mol·L −1 of ethylendiaminotetraoctic acid EDTA in 1:10 ratio (w:v-10 g of soil and 100 ± 0.5 mL of extractant). After extraction, the suspension was immediately filtered through thick filter paper [54].
Whole plants were collected before the harvest date (21 August) from six selected locations of each parcel (more precisely, from sampling points where data was taken using the GreenSeeker) using electric scissors just above the ground. The area of each site was 0.25 m 2 (0.5 × 0.5 m). The locations were bound by a U-shaped rod. The number of spikes from the cut plants, the total weight of the spikes and the weight of a thousand seeds (WTS), and the content of macro and micronutrients in the grains were assessed. The plant samples were analyzed in an external laboratory at the Crop Research Institute in Prague.

Statistical Analysis
The statistical analysis was processed in Statistica 13.5.0.17 software by TIBCO Software Inc., Tulsa, OK, USA. It was a mainly graphical processing of the resulting values as 2D graphs of averages and confidence intervals. There was also a statistical analysis, which specifically used the ANOVA tool and Tukey's HSD (honestly significant difference) test at the level of significance of 0.05. For basic data processing, MS Excel software was used. The zonal statistics tool in QGIS SW was used for image analysis.

Evaluation of Soil Parameters under Different Tillage
Basic physical soil parameters (moisture content, bulk density, total porosity, capillary porosity, and non-capillary) are given in Table 4. The content of major macronutrients in the soil is reported in Table 5.  The plots with the ancient Khorasan wheat variety contained a generally higher content of macronutrients in the soil (see Table 5). This fact is supported by various studies [41,44,55]. In general, this ancient wheat does not have demanding growing conditions, but at the same time, the grain is rich in nutrients compared to the modern Kabot variety. The results in Table 5 showed that the different tillage methods played a significant role in the ability of each individual variety to absorb nutrients and use them for growth. Khorasan wheat prospered better in conventionally tilled soil, where the bulk density of the soil was the smallest and the total porosity was the largest. From the point of view of the tillage variants, there are interesting results for the MTC variant, where the highest values of all the tillage variants were found for bulk density and initial moisture content (see Table 4).

Evaluation of Growth Parameters under Different Tillage
It can be seen in Table 6 that Khorasan wheat grains generally have higher values for both macro and micronutrients. Only the macronutrient potassium (K) shows no difference, and the values are variable. According to Wijngaard and Arendt [56], Khorasan wheat is classified as a healthy food due to its higher content of some nutrients, which agrees with this study. Likewise, Slavin et al. [57] point to a higher content of beneficial substances in the whole grains of Khorasan wheat. In addition, the interaction between all of the components, such as macro and micronutrients, omega 3-fatty acids, fiber, probiotic oligosaccharides, etc., gives these grains valuable nutritional value, which according to Maki and Phillips [58] can positively affect human health, as they can even reduce type 2 diabetes or the risk of obesity [13].
The development of the height of the wheat varieties and the yield parameters (number and total weight of spikes and weight of thousands of seeds) resulting from the different tillage methods are shown in Figure 2. The growth parameters of the selected wheat varieties cultivated under different tillage methods during the entire vegetation season were evaluated by spectral indices (Figure 3). Table 6. Content of macronutrients (phosphorus (P), magnesium (Mg), potassium (K) and calcium (Ca)) and micronutrients (copper (Cu), zinc (Zn), iron (Fe), manganese (Mn) and boron (B)) in grains of the Khorasan and Kabot wheat varieties and plots with different soil tillage.  There are large differences in the weight of one thousand seeds (WTS) between the varieties. The Khorasan variety has a much higher average WTS than Kabot. Khorasan has large grains and a high WTS value (usually around 50 g). In the case of this study, WTS was 62 g on average, which is consistent with the study of Grausgruber et al. [59]. They stated that in some cases, the WTS could be higher than 60 g. Differences in WTS values were also evident in the tillage variants, especially for variants 3 and 6 (MTD), where there is an increase in average values. However, this is a statistically insignificant There are large differences in the weight of one thousand seeds (WTS) between the varieties. The Khorasan variety has a much higher average WTS than Kabot. Khorasan has large grains and a high WTS value (usually around 50 g). In the case of this study, WTS was 62 g on average, which is consistent with the study of Grausgruber et al. [59]. They stated that in some cases, the WTS could be higher than 60 g. Differences in WTS values were also evident in the tillage variants, especially for variants 3 and 6 (MTD), where there is an increase in average values. However, this is a statistically insignificant difference within the wheat varieties. Khorasan seemed to be much more sensitive to the soil conditions (compared to Table 4-soil parameters). On the other hand, the number of spikes and the total weight of the spikes have a diametrically different trend than WTS. With the increasing number of spikes on the plot, the WTS decreased, but the total weight of the spikes was higher because the spikes contained more seeds. Where there were fewer numbers of spikes, there were fewer number of shoots. The individual plants then had more space and more nutrients to help them grow and develop. Tari [60] found that wheat is sensitive to water access during the stem elongation and heading stages. A lack of water in these stages can lead to a decrease in yield. Our results indicate sufficient water during the stem elongation phase. However, soil tillage played a significant role in the availability of water and nutrients for the growth and the resulting yield in the case of both wheat varieties. The number and weight of the spikes are in accordance with the trend in the soil capillary porosity.

Variant
Height was calculated from UAV images as DSM-DTM (see Table 3). During the tillering stage (18 May), the Khorasan variety crops were higher compared to the Kabot crops, indicating good germination and relatively good adaptation to the habitat conditions of this variety. Khorasan plants generally reach a higher height, most often around 127 cm [61]. The evaluation of the images also showed significant differences in plant height between the experimental parcels with different soil tillage. The tallest plants on all measurement dates for both varieties were the variant that was cultivated with conventional plowing technology (CT). These variants also had the highest total porosity and the lowest bulk density values (details in Table 4), which are the soil conditions most suitable for growth. Initial moisture content at the time of soil sampling was medium to low. On the contrary, the lowest Khorasan wheat height was found in the plot with minimum tillage; in the case of the Kabot wheat, the lowest height was found in the minimum and minimization tillage variants. Crop height for both varieties are in accordance with bulk density and has the of the opposite trend that can be seen in the total porosity values. The growth trends of varieties on plots with different tillage is steady. In the stages after BBCH 60 (2 and 14 July), the crops had already formed spikes, and the development of the grains in comparison with the WTS values is relatively well captured in Figure 2. Even though the WTS reached the highest value, the crop height in the last scanning (14 July) was lower than in the previous one (2 July). UAV height measurement is a relatively common and reliable method of capturing the vitality of the crop, along with derived spectral indices [62]. In the case of our measurements, we had a spatial resolution of about 2 cm·px −1 , which is sufficient to derive a relatively accurate model.
The development of NDVI values calculated from the MicaSense camera related to individual plots and dates showed a very balanced trend, especially in the later crop scanning sessions (see details in Figure 3). During the first UAV scanning (18 May), the Kabot variety in a plot with conventional tillage was the most lively. The stand was the best involved here, it had the best tillering ability (higher capillary porosity and lower bulk density of soil). These results are in line with the number of spikes (see Figure 2). The Khorasan variety had relatively the same development on all three plots with different tillage. However, a plot with conventional tillage proved to be the most suitable for the Khorasan variety in terms of evaluation using the spectral indices. Conversely, Glamoclija et al. [63] stated that the Khorasan variety best benefits in worse soil conditions in mountainous areas. According to Rajičić [64], the quality and yield of grains are mainly affected by the climatic conditions of the habitat. Our findings are in accordance with the study by Jelínek et al. [41]. They found that the NDVI values of wheat cultivars were higher in places with better water availability and nutrients, regardless of wheat variety.
The development of the NDVI values from the GreenSeeker was consistent with the results from the MicaSense camera. A different measurement method and a different wavelength range, as was the case in the study of Kumhálová & Matějková [23], probably caused variations in the results. They concluded that another reason of variation could be early crop scanning with different stand density. While all pixels on the parcels were taken into account from the UAV, the GreenSeeker was used only on sampling points (details in Figure 1).
The GNDVI spectral index is a reliable indicator of the chlorophyll content of the leaves. It usually has the same trend of development as NDVI [65], and this was confirmed in our study. As it can be seen in Figure 3, the stand was gradually maturing during the last date of scanning (14 July), and the chlorophyll content in the leaves was decreasing. The Kabot variety on the CT variant matures in delay, probably due to the large number of tillers.
into account from the UAV, the GreenSeeker was used only on sampling points (details in Figure 1).
The GNDVI spectral index is a reliable indicator of the chlorophyll content of the leaves. It usually has the same trend of development as NDVI [65], and this was confirmed in our study. As it can be seen in Figure 3, the stand was gradually maturing during the last date of scanning (14 July), and the chlorophyll content in the leaves was decreasing. The Kabot variety on the CT variant matures in delay, probably due to the large number of tillers. The CIR spectral index uses the red edge spectral band, which is very sensitive to the state of the cell structures of the plant tissues. The index has a very similar course to G/NDVI. However, a higher chlorophyll content on the CT plot of the Kabot variety was detected on 11 June, when the stand corresponded to the beginning of the heading growth stage. As in the case of the previously evaluated spectral indices, the reason was probably due to more beneficial soil conditions for the growth and better tillering of the Kabot variety.
The TGI index belongs to the RGB index group [66]. The TGI index was calculated from two sources-the more expensive 5 bands MicaSense RedEdge MX and the low-cost S.O.D.A. cameras. As it can be seen from Figure 3, the development of the TGI values was very similar for both data sources. The differences were probably caused due to different spectral characteristics (detail in Table 2). Both of the centers and ranges of the wavelengths of the individual bands differed. The variance of the values between the individual terms of the TGI index derived from images from the MicaSense camera was much smaller. However, as shown in Figure 3, both models of the TGI index development described the condition of the stand in the scanned terms well and are consistent with the other indices used. They are thus a very good alternative to NDVI, for example, if only an RGB camera is available. Jelínek et al. [18] solved a similar problem with spectral indices derived from two satellite images-Sentinel 2 and Landsat 8-for the evaluation of wheat varieties. Both sources were usable, but they showed slight differences.
Trends of the selected indices were in line with the content of micronutrients, especially Fe, Mn, and Zn. Although these indices are more often an indicator for the nitrogen content in the leaves or canopy structure, as described by Walsh and Shafian [67], they could be a good indicator of micronutrient content in the plant as well as in the grains. Pandey et al. [30] confirmed that the sufficient intake of micronutrients at right time and dosage can significantly increase the yield parameters of the wheat in terms of the number of spikes, height grains per spike, height of crops, among others. In the case of the selected varieties, the micronutrient content in the grains was in line with the height of both wheat varieties and affected the development of spikes in the Kabot variety. Khorasan primarily used microelements for the growth of plants and not for the development of tillers and spikes, but these micronutrients were highly concentrated in the grains themselves, which generally increased the nutritional value of the yield.
The development of canopy temperature depends on the density of the stand as well as on the other factors such as the characteristics of the thermal camera, the meteorological conditions, and the intensity of thermal radiation at the time of taking the pictures [68]. A literature review by Pandey et. al. [30] that specifically focuses on the use of thermal images for the purposes of water stress detection in crops. However, the relatively high temperatures (see Figure 3) in the first campaign were affected by very sparse vegetation, such as when the bare soil was visible. The CT plot of the Kabot variety generally showed the lowest temperatures in all of the campaigns, which is in line with the results from the spectral indices. The development of temperatures has the opposite course than the used spectral indices. Thus, the results the structure of the stand. Table 7 provides an overview of the coefficient of determination among the selected spectral indices and the measured crop and soil parameters for the Khorasan and Kabot varieties regardless of soil tillage variants. The significance of the coefficients of determination is highlighted according to Chráska [69]. It can be seen that the varieties differ in terms of their growth requirements. A significant difference between the varieties is mainly recorded in the area of spike development and grain weight and the ability of the selected spectral indices and sensors to capture details.

Conclusions
The results showed that each tillage technology played a significant role in providing suitable conditions for the growth of Khorasan ancient wheat. Although the literature review describes that it is not a demanding crop for mountain conditions, the results of our study showed that Khorasan is better suited for the conditions of conventional tillage with moldboard plowing at 20-25 cm. A sufficient supply of both macro and micronutrients (especially Zn, Fe and Mn) can significantly positively influence the height of both wheat varieties and the development of spikes in the Kabot variety. Khorasan used micronutrients primarily for the growth of plants and not for the development of tillers and spikes. These micronutrients were then highly concentrated in the grains of the Khorasan variety in the plot with conventional tillage (with low values of bulk density and high values of total porosity), which generally increased the nutritional value of the yield in this experimental plot. An additional advantage of this variety is that it does not deplete the soil.
Our results from four types of sensors have shown the advantages of combining different data sources that can complement and emphasize details. The CIR index with the red edge spectral band highlighted the characteristics of the stand regarding the state of the canopy structure in the beginning of heading growth stage of the Kabot variety. However, the results also show that the trend of developmental growing curves derived from different sensors was very similar regardless of the measurement method (area measurement-UAV images; point measurement-GreenSeeker). Sensors used in this study can be good indicators of the micronutrient content in the plant as well as in the grains. Nevertheless, more research is needed for a more accurate confirmation of these findings. It can be concluded that a low-cost RGB camera can provide relevant results, especially in cases where equipment that is more accurate is not available.