ORIGINAL RESEARCH article

Front. Plant Sci., 15 March 2019

Sec. Technical Advances in Plant Science

Volume 10 - 2019 | https://doi.org/10.3389/fpls.2019.00279

Applications of Unmanned Aerial Vehicle Based Imagery in Turfgrass Field Trials

  • 1. Crop and Soil Sciences Department, University of Georgia, Tifton, GA, United States

  • 2. Agronomy Department, University of Florida, Gainesville, FL, United States

  • 3. Turf Scout LLC., Greensboro, NC, United States

Abstract

Recent advances in remote sensing technology, especially in the area of Unmanned Aerial Vehicles (UAV) and Unmanned Aerial Systems (UASs) provide opportunities for turfgrass breeders to collect more comprehensive data during early stages of selection as well as in advanced trials. The goal of this study was to assess the use of UAV-based aerial imagery on replicated turfgrass field trials. Both visual (RGB) images and multispectral images were acquired with a small UAV platform on field trials of bermudagrass (Cynodon spp.) and zoysiagrass (Zoysia spp.) with plot sizes of 1.8 by 1.8 m and 0.9 by 0.9 m, respectively. Color indices and vegetation indices were calculated from the data extracted from UAV-based RGB images and multispectral images, respectively. Ground truth measurements including visual turfgrass quality, percent green cover, and normalized difference vegetation index (NDVI) were taken immediately following each UAV flight. Results from the study showed that ground-based NDVI can be predicted using UAV-based NDVI (R2 = 0.90, RMSE = 0.03). Ground percent green cover can be predicted using both UAV-based NDVI (R2 = 0.86, RMSE = 8.29) and visible atmospherically resistant index (VARI, R2 = 0.87, RMSE = 7.77), warranting the use of the more affordable RGB camera to estimate ground percent green cover. Out of the top ten entries identified using ground measurements, 92% (12 out of 13 in bermudagrass) and 80% (9 out of 11 in zoysiagrass) overlapped with those using UAV-based imagery. These results suggest that UAV-based high-resolution imagery is a reliable and powerful tool for assessing turfgrass performance during variety trials.

Introduction

Constraints such as high demands of time and labor during field phenotyping limit the ability of turfgrass breeders to collect more comprehensive data during early stages of selection and later in advanced trials. Recent improvements in remote sensing technologies provide opportunities to mitigate this bottleneck for future breeding advances. Non-invasive remote sensing methods such as digital image analysis and spectral reflectance have been widely used for quantifying turfgrass cover and quality (Menges et al., 1985; Richardson et al., 2001; Jiang and Carrow, 2007; Xiong et al., 2007). Vegetation indices such as NDVI and ratio vegetation index (RVI), calculated from reflectance at red and near infrared (NIR) bands, have been validated to predict turfgrass health in previous studies (Fitz-Rodríguez and Choi, 2002; Jiang and Carrow, 2007; Bremer et al., 2011). UAV platforms are evolving rapidly and offer advantages over other vehicles when used for remote sensing. When compared with satellite-based and ground-based remote sensing, UAV-based imagery has higher spatial resolution (1–2 cm per pixel) (Xiang and Tian, 2011) and thus more powerful statistical analytics at an affordable price. Both spatial and temporal resolution are important for turfgrasses because they can experience intermittent drought stresses in sandy soils in only three or 4 days without precipitation (Zhang et al., 2018).

In agricultural studies, the use of UAV-based imagery is increasing. High-resolution digital images have been acquired from aerobatic model aircrafts for estimating the nutrient status and crop biomass of corn (Zea mays L.), alfalfa (Medicago sativa L.), and soybeans [Glycine max (L.) Merr.] (Hunt et al., 2005). Primicerio et al. (2012) developed a six-rotor aerial platform with a multi-spectral camera to map the vineyard vigor of wine grapes (Vitis vinifera L.). In more recent studies, UAV systems were applied to monitor crop growth in wheat (Triticum spp.) (Lelong et al., 2008; Torres-Sánchez et al., 2014), predict yield in rice (Oryza Sativa L.) (Zhou et al., 2017), detect disease in potato (Solanum tuberosum L.) (Sugiura et al., 2016), and detect weeds in sunflower (Helianthus annuus L.) fields (Torres-Sánchez et al., 2013).

Recent studies with UAVs on turfgrass have also shown promising results. Xiang and Tian (2011) used an unmanned helicopter to monitor turfgrass response after glyphosate application and found only a 1.5% difference in the estimation of herbicide damage between aerial images and ground surveys. Caturegli et al. (2016) used UAV-based multi-spectral imagery to estimate the nitrogen status of hybrid bermudagrass (C. dactylon L. × C. transvaalensis Burtt-Davy), zoysiagrass (Z. matrella L. Merr.), and seashore paspalum (Paspalum vaginatunt Swartz.), concluding that UAV imagery can adequately assess the spatial variability of nitrogen status for these turfgrass species in large areas such as golf courses and sod farms. Another study assessed the health of creeping bentgrass (Agrostis stolonifera L.) using a RGB camera and a modified NIR camera under different mowing heights (Sommer et al., 2017). It was reported that color indices calculated from RGB images were better correlated (R2∼0.8) with ground truth data compared to mNDVI from the modified NIR camera (R2∼0.6). To our knowledge, no investigation has been conducted regarding the use of UAV-based imagery on turfgrass variety trials. Prior to being used by turfgrass breeders in variety trials, data extracted based on UAV images need to be compared with ground measurements. Therefore, the overall goal of the present study was to assess the potential use of visual (RGB) and multispectral images collected with a UAV platform in replicated turfgrass field trials. The objectives were (1) to examine the correlation between UAV-based measurements and ground measurements; (2) to determine the feasibility of developing a general model to predict ground measurements for two turfgrass species in different sampling dates; and (3) to assess if the information extracted from UAV-based imagery can help turfgrass breeders to make better decisions (genotype rankings) during selection processes.

Materials and Methods

UAV System and Cameras

A Solo quadcopter (3D Robotics, Berkeley, CA, United States) with vertical take-off and landing capabilities was used to collect aerial images for this study (Figure 1). The UAV system, equipped with four brushless motors, can fly by either remote control or autonomously with Global Positioning Systems (GPS) and waypoint navigation system. The whole system consists of the drone, controller, and a ground station with software for mission planning, flight control and telemetry (Torres-Sánchez et al., 2013). The ground station interface between the pilot-in-command and the UAV has supporting software that implements a flight plan and monitors the flight. The software also has a telemetry system that collects relevant flight data and information including GPS position data and flight time. This information is known as the telemetry log and is useful during image processing. Two persons were involved during the flight mission, one being the remote pilot who was in charge of taking-off and landing the UAV by activating the programmed flight plan during the flight operations, and the second who acted as a visual observer that kept watch of the UAV for potential collision threats with other air traffic. Two cameras mounted separately on two similar quadcopters were used in this study. The first was a GoPro Hero 4 (GoPro, Inc. San Mateo, CA, United States) visual camera which acquires 7-megapixel images in true color (Red, R; Green, G, and Blue, B, bands) with 8-bit radiometric resolution (Figure 1B). The second camera was a Parrot Sequoia (MicaSense, Seattle, WA, United States) multispectral camera that measures at four narrow spectral bands (green: 530–570 nm; red: 640–680 nm; red edge: 730–740 nm; NIR: 770–810 nm) (Figure 1C,D).

FIGURE 1

Study Site and UAV Flights

The study site was part of a multi-state breeding project funded by U.S. Department of Agriculture – Specialty Crop Research Initiative. Two advanced breeding trials of bermudagrass and zoysiagrass were planted using plugs in June 2016 on a loamy sand (Tifton-Urban land complex, pH 5.3) at the University of Georgia Tifton Campus. Field plots were arranged as a randomized complete block design with three replications. The plot size was 1.8 by 1.8 m for bermudagrass and 0.9 m by 0.9 m for zoysiagrass. Forty advanced lines from two breeding programs (University of Georgia and Oklahoma State University) and four commercial cultivars (“Celebration,” “Latitude 36,” “TifTuf,” and “Tifway”) were included in the bermudagrass trial. In the zoysiagrass study, 40 advanced lines from two breeding programs (Texas A&M and University of Florida) and 3 commercial cultivars (“Empire,” “Palisades,” and “Zeon”) were tested. The UAV-based RGB and multispectral images were taken on 28 September 2017, 5 April 2018, 8 May 2018, and 18 June 2018 at 0930 h with no cloud cover. The flight altitude was 30 m for the RGB camera, resulting in image resolution of 2.2 cm per pixel. For the multispectral camera, the flight altitude was 46 m (to aid in image stitching later because the multispectral camera has a narrower field of view) which resulted in an image resolution of 4.3 cm per pixel. The UAV speed was set to 4.5 ms-1 for all flights. All flights were controlled using a free software, “Tower Beta” (DroidPlanner, 2018), installed on a smart device (Tablet) and by creating a flight plan for UAV equipped with each camera prior to any flight. The software maintained the desired flight altitude and flight speed during the flight while capturing the images at 80% front and side overlaps.

Image Preprocess and Data Acquisition

Geotagging of the RGB images collected by the GoPro Hero 4 was completed using the information provided by the telemetry log saved on the smart device and using the open source software, ‘Mission Planner’ (Oborne, 2018). Multispectral images were geotagged during flight with an on-board GPS in the Parrot Sequoia camera. Geotagged images were then processed and stitched in Pix4Dmapper Pro 4.2.27 (Pix4D SA, Lausanne, Switzerland) to generate an orthomosaic consisting of information collected for individual bands (RGB, NIR, and Red). Standard templates of “Ag RGB” and “Ag Multispectral” in Pix4Dmapper were used for stitching RGB and multispectral images respectively. The georeferenced orthomosaic was exported in a TIFF format for further analysis in ArcGIS (Esri, Redlands, CA, United States). For data analysis, a shape file consisting of the individual field plot information was created in ArcMap version 10.4.1 (Esri, Redlands, CA, United States). Data were extracted within each polygon (each polygon represented a plot) using the ArcMap feature in zonal statistics (Figure 2). Three different vegetation indices including NDVI, Green NDVI (GNDVI), and Normalized Difference Red Edge Index (NDRI) were calculated from the multispectral images (Table 1). Six color indices were calculated with normalized values for R, G, and B bands from the digital image (Torres-Sánchez et al., 2014; Saberioon et al., 2014; Zhou et al., 2017). The following normalization scheme was applied to the color indices:

FIGURE 2

Table 1

IndicesEquationReference
NDVI: normalized difference vegetation index(RNIR−RRed)/(RNIR+RRed)Lee et al., 2011
GNDVI: green normalized difference vegetation index(RNIR−RGreen)/(RNIR+RGreen)Gitelson et al., 1996
NDRI: normalized red edge index(RNIR−RRededge)/(RNIR+RRededge)Mutanga and Skidmore, 2004
NDI: normalized difference index(g−r)/(g+r)Woebbecke et al., 1995
ExG: excessive green index2g−r−bWoebbecke et al., 1995
ExR: excessive red1.4r−gMeyer and Neto, 2008
ExGR: excessive green index minus excess red index3g−2.4r−bMeyer and Neto, 2008
VARI: visible atmospherically resistant index(g−r)/(g+r−b)Gitelson et al., 2002
GLI: green leaf index(2g−b−r)/(2g+b+r)Louhaichi et al., 2001

Equations for color indices used in the study.

where R, G, and B are the values of the red, green, and blue bands, respectively. Those six color indices included normalized difference index (NDI), excessive green index (ExG), excessive red (ExR), excessive green index minus excess red index (ExGR), VARI, and green leaf index (GLI) (Table 1).

Ground Measurements

Ground truth measurements including turfgrass quality (TQ), percent green cover, and canopy spectral reflectance of the plots were collected on the same day within an hour or two of the UAV-based imagery data collection. Visual ratings of TQ were based on the National Turfgrass Evaluation Program using a 1 to 9 scale (9 = excellent performance and 1 = poor performance, 6 = minimum acceptable quality) (Morris and Shearman, 2008). Percent green cover was estimated from digital images collected using a digital camera (Powershot G5; Canon, Tokyo, Japan) mounted to an enclosed photo box (56 cm by 56 cm) with four 9-W compact fluorescent lamps (TCP; Lighthouse Supply, Bristol, VA, United States). Each image was analyzed using SigmaScan Pro (version 5.0; Systat Software, San Jose, CA, United States) for percent green cover (0–100%) using a hue range from 60 to 120 and saturation range from 10 to 100 as outlined by Richardson et al. (2001). Canopy spectral reflectance values were measured using a Crop Circle ACS470 sensor (Holland Scientific, Lincoln, NE, United States), equipped with a decimeter level GPS (Raven Industries, Sioux Falls, SD, United States). The spectral sensor, with active light source, measured light reflectance in three spectral bands centered on 550 nm (green), 650 nm (red), and 730 nm (NIR). The system was mounted to a mobile cart at 61 cm above ground with a target ground area of 35 by 6.4 cm. Data were collected and processed using TurfScout platform (TurfScout, Greensboro, NC), where NDVI values were calculated within the program.

Data Analysis

All data were subjected to analysis of variance using SAS 9.4 (SAS Institute Inc., Cary, NC, United States). For each parameter, Fisher’s protected LSD at 0.05 probability level was used to separate significant means and to mark the top statistical group in bermudagrass and zoysiagrass entries. Turf Performance Index was calculated by summing up the number of times an entry entered the top statistical group (Wherley et al., 2013). Spearman’s rank correlation and linear regression were performed between UAV-based measurements and ground measurements in SAS using the CORR and Generalized Linear Model (GLM) procedures, respectively. Graphs were generated using SigmaPlot 14 (Systat Software, Inc. Point Richmond, CA, United States). Box plots were generated using ggplot2 in Rstudio (Wickham, 2016).

Results and Discussion

Distribution of Ground Truth Measurements and UAV-Based Measurements

Ground percent green cover, TQ, ground and UAV-based NDVI in bermudagrass ranged from 10 to 100%, 2 to 8, 0.17 to 0.66, and 0.35 to 0.80, respectively (Figure 2). The aforementioned parameters in zoysiagrass ranged from 10 to 80%, 2 to 6, 0.10 to 0.53, and 0.28 to 0.74, respectively. Parameters such as percent green cover and NDVI were affected by different sampling dates due to the combination of temperature, irradiance, and precipitation. In general, sampling dates in September 2017 and April 2018 were drier than the other two dates. Weather data were obtained from a weather station located about 1 kilometer away from the study site (Figure 3). Although no soil moisture data and leaf water content was taken, virtual symptoms of leaf wilting under varying degrees was observed in the plots due to the combination of evapotranspiration demand and irrigation frequency (once every 2 weeks). The study site was originally established to evaluate long-term drought resistance in advance experimental lines, and the irrigation frequency was designed to encourage periodic drought. Zoysiagrasses responded more rapidly compared to bermudagrasses in these drier environments as indicated by percent green cover and NDVI (Figure 4). This was likely due to the relative poor drought resistance of zoysiagrass as compared to that of bermudagrass (Zhang et al., 2018). Seasonal effects and interaction between turfgrass species and environments are not the focus of this study. Data over different seasons and species were collected to better elucidate the relationship between ground truth and UAV-based measurements and to build a more robust model.

FIGURE 3

FIGURE 4

Correlation Between Ground Measurements and UAV-Based Measurements

For multispectral imagery, ground truth measurements including TQ, percent green cover, and ground NDVI were positively correlated with the UAV-based vegetation indices NDVI (0.64, 0.89, and 0.92, P < 0.0001), GNDVI (0.59, 0.82, and 0.88, P < 0.0001), and NDRI (0.60, 0.68, and 0.76, P < 0.0001), respectively (Table 2). For RGB imagery, the aforementioned ground truth measurements were positively correlated with the UAV-based color indices VARI (0.59, 0.88, and 0.88, P < 0.0001), NDI (0.48, 0.56, and 0.60, P < 0.0001), ExG (0.44, 0.66, and 0.67, P < 0.0001), ExR (−0.28, −0.70, and −0.68, P > 0.0001), ExGR (0.53, 0.79, and 0.79, P < 0.0001), and GLI (0.63, 0.79, and 0.81, P < 0.0001), respectively. Further results and discussion is focused on UAV-based NDVI and VARI due to their higher correlation with the ground truth measurements in the study.

Table 2

Ground measurements
UAV-based Measurements
PGCNDVI.GNDVI.UGNDVI.UNDRI.UVARINDIExGExRExGRGLI
TQ†0.650.660.640.590.600.590.480.44−0.280.530.63
PGC0.880.890.820.680.880.560.66−0.700.790.79
NDVI.G0.920.880.760.880.600.67−0.680.790.81
NDVI.U0.980.760.910.490.80−0.810.900.90
GNDVI.U0.730.850.380.82−0.820.900.89
NDRI.U0.700.410.67−0.590.720.76
VARI0.690.69−0.78−0.780.84
NDI0.12−0.190.300.38
ExG−0.890.950.93
ExR−0.92−0.82
ExGR0.97

Spearman’s rank correlation coefficients between ground measurements (TQ, PGC, and NDVI.G) and UAV-based measurements (NDVI.U, GNDVI.U, NDRI.U, VARI, NDI, ExG, ExR, ExGR, and GLI) in bermudagrass and zoysiagrass variety trials from September 2017 to June 2018 in Tifton, GA.

TQ, turfgrass quality; PGC, percent green cover; NDVI.G, normalized difference vegetation index at ground level; NDVI.U, UAV-based NDVI; RNDVI.U, UAV-based green NDVI; NDRI.U, UAV-based normalized difference red edge index; VARI, UAV-based visible atmospherically resistant index; NDI, normalized difference index; ExG, excess green; ExR, excess red; ExGR, excess green minus excess red; GLI, green leaf index. †All coefficients are significantly at 0.001 P level except when labeled with “ns.”

Predicting Ground Measurements Using UAV-Based Measurements

Regression analysis indicated that ground percent green cover in bermudagrass can be predicted using UAV-based NDVI (PGC = 217.11 × NDVI. U – 73.27, R2 = 0.88, RMSE = 7.37) and VARI (PGC = 567.35 × VARI + 33.86, R2 = 0.89, RMSE = 6.94) (Table 3 and Figure 5) with only ∼7% difference according to the root mean square errors. Similarly, models fitted to predict ground percent green cover in zoysiagrass had R2 of 0.73 (PGC = 140.39 × NDVI.U – 30.88, RMSE = 7.70) and 0.69 (PGC = 693.31 × VARI + 35.02, RMSE = 8.24) using UAV-based NDVI and VARI, respectively. Ballesteros et al. (2014) also found that the color index VARI was a good predictor for estimating percent green cover in maize (PGC = 293.8 × VARI + 67.13, R2 = 0.95, RMSE = 6.71) and onion (Allium cepa L.) (PGC = 279.1 × VARI + 68.13, R2 = 0.82, RMSE = 6.29). Models of row crops like maize and onion had steeper slopes in predicting percent green cover than turfgrass in our study, which can be caused by multiple factors. Firstly, maize and onion crops have wider leaves than turfgrass so percent green cover maybe more responsive to the change in VARI. Secondly, for row crops, percent green cover is more of an interest before canopy closure and this would likely to impact the threshold settings in obtaining percent green cover. In turfgrass, percent green cover is used throughout every stage of growth. Therefore, results on row crops may not be directly applicable for turfgrass. When data from both the bermudagrass and zoysiagrass experiments were combined, ground NDVI was still accurately predicted using UAV-based NDVI (R2 = 0.84–0.86, RMSE = 0.03), indicating that UAV-based sensor information is consistent across at least these two turfgrass species. Sommer et al. (2017) reported that color indices (R2 = 0.80) from a RGB camera were better correlated with ground measurement NDVI than mNDVI (R2 = 0.60) from a modified NIR camera. In our study, VARI was a good predictor for ground percent green cover in our study, supporting the use of a more affordable digital camera for data collection if ground percent green cover is the interest. Multi-spectral cameras may provide more information for the detection of weed pressure (Torres-Sánchez et al., 2013), disease incidence (Mahlein, 2016; Brodbeck et al., 2017), and drought stress (Zarco-Tejada et al., 2013).

Table 3

Adjusted R2p-valueRMSE†
Bermudagrass
PGC (%) = 217.11 × NDVI.U − 73.270.88∗∗∗7.37
PGC (%) = 567.35 × VARI + 33.860.89∗∗∗6.94
TQ = 9.19 × NDVI.U − 0.580.68∗∗∗0.58
TQ = 22.92 × VARI + 4.000.63∗∗∗0.62
NDVI.G = 0.85 × NDVI.U − 0.080.86∗∗∗0.03
Zoysiagrass
PGC (%) = 140.39 × NDVI.U − 30.880.73∗∗∗7.70
PGC (%) = 693.31 × VARI + 35.020.69∗∗∗8.24
TQ = 4.41 × NDVI.U + 1.780.13∗∗∗1.04
TQ = 14.17 × VARI +3.860.05∗∗1.09
NDVI.G = 0.73 × NDVI.U − 0.060.84∗∗∗0.03
Species Combined
PGC (%) = 177.49 × NDVI.U − 48.740.86∗∗∗8.29
PGC (%) = 571.61 × VARI + 34.380.87∗∗∗7.77
TQ = 7.47 × NDVI.U + 0.410.49∗∗∗0.87
TQ = 23.21 × VARI + 3.920.47∗∗∗0.89
NDVI.G = 0.91 × NDVI.U − 0.130.90∗∗∗0.03

Regressions between ground measurements including percent green cover (PGC), turfgrass quality (TQ), and ground normalized difference vegetation index (NDVI.G) and UAV-based measurements including NDVI.G and visible atmospherically resistant index (VARI).

∗∗ and ∗∗∗ indicate model significance level at 0.01, 0.001 P level. †RMSE is the standard deviation of the error.

FIGURE 5

All models for the sensor-based prediction of TQ were more poorly fitted (R2 = 0.63–0.68 for bermudagrass and 0.05–0.13 for zoysiagrass) than those for percent green cover. This agrees with previous findings that TQ is more difficult to predict based on percent green cover or NDVI (Sullivan et al., 2017), which was largely due to its subjective nature and greater coefficients of variation than optical sensor measurements (Bell et al., 2002; Lee et al., 2011; Bremer et al., 2011). Although TQ, based on 1–9 scale, is still an important parameter to be used when evaluating turf performance, objective measurements such as percent green cover and NDVI are being adopted more often for their reliable and repeatable estimation of quality (Richardson et al., 2001, 2008, Bremer et al., 2011). Results from our study demonstrated great potential for UAV-based imagery to be used in assessing turf performance.

Bremer et al. (2011) reported that separate models for different grass species may be required to predict TQ with ground NDVI. To date we have not found a report on whether ground-based NDVI or percent green cover of turfgrass can be predicted using UAV-based measurements, or whether these relationships are species dependent. In our study, the model based on combined data from both species had an increased R2 (0.90) and equal root mean square of error (RMSE = 0.03) in the prediction of ground NDVI using UAV-based NDVI compared to models fitted from individual species (Table 3). Using a species combined model to predict percent green cover in bermudagrass resulted in an increased error of 0.92% in UAV-based NDVI and 0.83% in VARI. The increase in errors were smaller when predicting zoysiagrass percent green cover (0.59% for NDVI.U and 0.47% for VARI). These results suggest that developing a general model to predict ground NDVI using a UAV platform for both bermudagrass and zoysiagrass, and possibly more species, without a significant sacrifice of accuracy is possible. However, before building a general model to increase the efficiency of data collection, certain issues need to be addressed and the limitation of the current model should be taken into account. Several UAV-based platforms and sensors need to be tested for the same purpose in order to determine if the model is platform/sensor dependent. The model may be limited when it comes to a different turfgrass species with broader leaf such as St. Augustinegrass [Stenotaphrum secundatum (Walt.) Kuntze] and under less ideal weather conditions, such as overcast and cloudy days. Compared to different species and weather conditions, repeated time scale is less of a concern for fitting the model. Further investigation, addressing those limitations, should be completed to improve the fit of general models to predict ground measurements from UAV-based systems.

Entry Ranking Using Ground Measurements and UAV-Based Measurements

Depending on the trait or genes targeted for improvement, typically the top performing 10–20% (Hanna et al., 2013) of genotypes are advanced to the next cycle of selection during early generation turfgrass breeding variety trials, which in the context of this study represents the best 5–10 performers. There was a 100% overlap in the identification of the top 5 bermudagrass entries using ground measurements (percent green cover and NDVI) and UAV-based measurements (NDVI and VARI) (Table 4), including Tif16118 (1st vs. 1st), Tif16116 (1st vs. 5th), Tif16115 (3rd vs. 2nd), TifTuf (3rd vs. 3rd), Tif16110 (3rd vs. 5th), and Tif16112 (3rd vs. 5th). An overlap of 93% was observed in identifying the top 10 bermudagrass entries between two groups of measurements, meaning 12 out of 13 identified by ground measurements were the same entries identified using UAV-based measurements. Less overlap (55%) was found in identifying the top 5 zoysiagrass entries using both types of measurements, but 80% of the top 10 zoysiagrass entries corresponded to ground and UAV-based rankings. The 55% overlap of the top 5 zoysiagrasses was due to the lack of statistical separation among 5th ranking entries using ground measurements, including DALZ1614, Empire, FZ1440, FAES1335, FAES1322, DALZ1625, and DALZ1410.

Table 4

Ground
UAV
Ground
UAV
BermudagrassTPI†RankTPI‡RankZoysiagrassTPIRankTPIRank
Tif161188181DALZ14097181
Tif161168155DALZ16047181
Tif161157372DALZ16067155
TifTuf7363Palisades71410
Tif161107355DALZ16145573
Tif161127355Empire5564
OSU13376763FZ14405555
Tif1610667310FAES133555410
OSU143967310FAES132255410
OSU140651055DALZ162555320
Tif1610251049DALZ141055232
Tif16104510310FZ142741255
Tif16117510027DALZ1605412410
OSU1408414310FZ1436412320
Tif16105414310FZ1327412232
Tif16114414118DALZ1611316410
Tif16108316310DALZ1311316320
OSU1433316118DALZ1603316320
OSU1418316118FAES1313316320
Tif16103316027DALZ1609316320
Tif16107316027FZ1309316320
Tif16101316027DALZ1408316232
Tif16113223216DALZ1313316232
Tifway223216FZ1252316232
OSU1412223118FZ1410316232
Latitude36223118FZ1407226410
OSU1402223118FAES1336226410
OSU1403223118FZ1429226410
OSU1417223027DALZ1310226320
OSU1257223027Zeon226232
Tif16119223027DALZ1626226140
OSU1435223027FZ1337226140
OSU1414133118FZ142213355
OSU1425133027DALZ1613133410
Tif16109133027FAES1337133320
Tif16120133027DALZ1601133320
OSU1423133027DALZ1607133320
Tif16111133027FZ1333133140
OSU1318038118DALZ131403955
Celebration038027DALZ1615039410
OSU1310038027FAES1307039320
OSU1409038027FAES1319039232
OSU1415038027FAES1312039140
OSU1420038027

Turf performance index (TPI) and genotypic ranking using ground measurements (percent green cover and ground NDVI) and UAV-based measurements (NDVI and VARI) in bermudagrass and zoysiagrass from September 2017 to June 2018.

†Ground TPI summed up the number of times an entry entered the top statistical group in percent green cover and ground NDVI. ‡UAV TPI summed up the number of times an entry entered the top statistical group in UAV-based NDVI and VARI.

Stress tolerance is a common interest in turfgrass variety trials (Chandra, 2015), with performance during drought stress being of the utmost importance. In our study, both ground and UAV-based measurements were used to rank the entries from all dates combined and also only during the two drier dates (September 2017 and May 2018). There was an 83% overlap in the identification of the top 5 entries in both species over all dates, and when the two drier dates were analyzed separately (Table 5). A few entries had an increased ranking under conditions of reduced soil moisture, indicating some mechanism of drought tolerance or avoidance. These entries were Tif16115 (from 2nd to 1st) and TifTuf (from 2nd to 1st) bermudagrass, DALZ1614 (from 4th to 2nd), FA1427 (from 7th to 5th), and FAES1322 (from 9th to 5th) zoysiagrass.

Table 5

All dates
Drought
All dates
Drought
BermudagrassTPI†RankTPI‡RankZoysiagrassTPIRankTPIRank
Tif16118191101DALZ1409191101
Tif16115172101DALZ160419192
TifTuf172101DALZ160616385
OSU133716469DALZ161415492
Tif1611016485Empire14592
Tif1611616494Palisades14578
Tif1611215769FAES133513778
OSU140614876FZ142713785
Tif1610612976DALZ1605129610
OSU14391110315FAES132212985
Tif16102111076FZ1440129610
Tif161041110512DALZ13111012419
OSU1408913413DALZ16111012610
Tif1610591369DALZ16251012513
OSU1412615315FZ14221012419
OSU1433615315DALZ1410916513
Tif16108615121FAES1313916513
Tif16113615121FAES1336916513
Tif16114615121FZ1436916330
Tif16117615413DALZ1314820513
Tifway521218DALZ1408820419
OSU1417422218DALZ1603820330
OSU1418422121DALZ1609820419
Tif16103422121FZ1327820419
Tif16107422033FZ1407820419
Latitude36326121FZ1429820419
OSU1257326121DALZ1310727330
OSU1402326121DALZ1313727419
OSU1403326121DALZ1613727419
Tif16101326121FAES1337727237
Tif16119326218FZ1252727419
OSU1414232033FZ1309727513
OSU1425232121Zeon633419
OSU1435232121DALZ1615534330
Tif16109232033FAES1307534330
Tif16120232033FAES1319534330
OSU1318137033FZ1410534237
OSU1423137033DALZ1601438237
Tif16111137033DALZ1607438330
Celebration040033DALZ1626438237
OSU1310040033FZ1337438237
OSU1409040033FAES1312342142
OSU1415040033FZ1333342142
OSU1420040033

Turf performance index (TPI) based on TQ, percent green cover, ground NDVI, UAV-based NDVI and VARI and genotypic ranking based on TPI in bermudagrass and zoysiagrass from September 2017 to June 2018.

†All dates TPI summed up the number of times an entry entered the top statistical group for all five parameters on all four dates. ‡Drought dates TPI summed up the number of times an entry entered the top statistical group for all five parameters on September 2017 and May 2018.

In summary, predicting objective variables such as ground percent green cover and NDVI is more reliable than subjective visual TQ. Data (VARI and NDVI) extracted based on high-resolution UAV imagery provide accurate estimates of ground measurement such as percent green cover and NDVI. Given the advantages in field of view, spatial and temporal resolution, UAV-based imagery would equip turfgrass breeder with powerful tools for assessing turfgrass performance, which would greatly increases the efficiency of data collection in relatively large trials (>1,000 plots). Before developing a general model for predicting ground measurement for turfgrass, further investigation of more UAV platforms/sensors on different turfgrass species under varying weather conditions and locations is needed. Further research using more advanced thermal and hyperspectral sensors should be conducted during periods of drought stress, as well as for disease detection, to determine their value against more readily available UAV-based imagery platforms.

Statements

Author contributions

JZ and BS designed the experiments. JZ, SV, and WP collected UAV images. JZ and BS took the ground measurements. JZ processed the images and analyzed the data. DS contributed the data analysis on the ground measurement of NDVI. KK provided critical reviews and supported the costs of the study.

Funding

The authors would like to acknowledge the funding source from United States Department of Agriculture National Institute of Food and Agriculture Specialty Crop Research Initiative. Grant number 2015-51181-24291.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    BallesterosR.OrtegaJ. F.HernándezD.MorenoM. A. (2014). Applications of georeferenced high-resolution images obtained with unmanned aerial vehicles. Part II: application to maize and onion crops of a semi-arid region in Spain.Precis. Agric.15593614. 10.1007/s11119-014-9357-6

  • 2

    BellG. E.MartinD. L.WieseS. G.DobsonD. D.SmithM. W.StoneM. L.et al (2002). Vehicle-mounted optical sensing.Crop Sci.42197201. 10.2135/cropsci2002.1970

  • 3

    BremerD. J.LeeH.SuK.KeeleyS. J. (2011). Relationships between normalized difference vegetation index and visual quality in cool-season turfgrass: II. factors affecting NDVI and its component reflectances.Crop Sci.5122192227. 10.2135/cropsci2010.12.0729

  • 4

    BrodbeckC.SikoraE.DelaneyD.PateG.JohnsonJ. (2017). Using unmanned aircraft systems for early detection of soybean diseases.Adv. Anim. Biosci.8802806. 10.1017/S2040470017001315

  • 5

    CaturegliL.CornigliaM.GaetaniM.GrossiN.MagniS.MigliazziM.et al (2016). Unmanned aerial vehicle to estimate nitrogen status of turfgrasses (X Hu, Ed.).PLoS One11:e0158268. 10.1371/journal.pone.0158268

  • 6

    ChandraA. (2015). “Plant genetics and genomics to improve drought and salinity tolerance for sustainable turfgrass production in the Southern United States,” in Proceedings of the ASA, CSSA and SSSA International Annual Meetings, Minneapolis, MN.

  • 7

    DroidPlanner (2018). Tower Beta. [Mobile Application Software]. Available at: https://play.google.com

  • 8

    Fitz-RodríguezE.ChoiC. Y. (2002). Monitoring turfgrass quality using multispctral radiometry.Trans. ASAE45865871. 10.13031/2013.8839

  • 9

    GitelsonA. A.KaufmanY. J.MerzlyakM. N. (1996). Use of a green channel in remote sensing of global vegetation from EOS-MODIS.Remote Sens. Environ.58289298. 10.1016/S0034-4257(96)00072-7

  • 10

    GitelsonA. A.KaufmanY. J.StarkR.RundquistD. (2002). Novel algorithms for remote estimation of vegetation fraction.Remote Sens. Environ.807687. 10.1016/S0034-4257(01)00289-9

  • 11

    HannaW.RaymerP.SchwartzB. M. (2013). “Warm-season grasses: biology and breeding,” in Turfgrass: Biology, Use and Management, edsCisarJ.HorganB.BonosS. (Madison, WI: ASA), 543590.

  • 12

    HuntE. R.CavigelliM.DaughtryC. S. T.McmurtreyJ. E.WalthallC. L. (2005). Evaluation of digital photography from model aircraft for remote sensing of crop biomass and nitrogen status.Precis. Agric.6359378. 10.1007/s11119-005-2324-5

  • 13

    JiangY.CarrowR. (2007). Broadband spectral reflectance models of turfgrass species and cultivars to drought stress.Crop Sci.4716111618. 10.2135/cropsci2006.09.0617

  • 14

    LeeH.BremerD. J.SuK.KeeleyS. J. (2011). Relationships between normalized difference vegetation index and visual quality in turfgrasses: effects of mowing height.Crop Sci.51323332. 10.2135/cropsci2010.05.0296

  • 15

    LelongC.BurgerP.JubelinG.RouxB.LabbéS.BaretF.et al (2008). Assessment of unmanned aerial vehicles imagery for quantitative monitoring of wheat crop in small plots.Sensors835573585. 10.3390/s8053557

  • 16

    LouhaichiM.BormanM. M.JohnsonD. E. (2001). Spatially located platform and aerial photography for documentation of grazing impacts on wheat.Geocarto Int.166570. 10.1080/10106040108542184

  • 17

    MahleinA. K. (2016). Plant disease detection by imaging sensors – parallels and specific demands for precision agriculture and plant phenotyping.Plant Dis.100241251. 10.1094/PDIS-03-15-0340-FE

  • 18

    MengesR. M.NixonP. R.RichardsonA. J. (1985). Light reflectance and remote sensing of weeds in agronomic and horticultural crops.Weed Sci.33569581. 10.1017/s0043174500082862

  • 19

    MeyerG. E.NetoJ. C. (2008). Verification of color vegetation indices for automated crop imaging applications.Comput. Electron. Agric.63282293. 10.1016/J.COMPAG.2008.03.009

  • 20

    MorrisK. N.ShearmanR. C. (2008). NTEP Turfgrass Evaluation Guidelines.Beltsville, MD: National Turfgrass Evaluation Program.

  • 21

    MutangaO.SkidmoreA. K. (2004). Narrow band vegetation indices overcome the saturation problem in biomass estimation.Int. J. Remote Sens.2539994014. 10.1080/01431160310001654923

  • 22

    OborneM. (2018). Mission Planner. Available at: http://ardupilot.org/planner/index.html

  • 23

    PrimicerioJ.DiGennaroF.FiorilloE.GenesioL.LugatoE.MateseA.et al (2012). A flexible unmanned aerial vehicle for precision agriculture.Precis. Agric.13517523. 10.1007/s11119-012-9257-6

  • 24

    RichardsonM. D.KarcherD. E.HignightK.RushD. (2008). Drought tolerance and rooting capacity of kentucky bluegrass cultivars.Crop Sci.4824292436. 10.2135/cropsci2008.01.0034

  • 25

    RichardsonM. D.KarcherD. E.PurcellL. C. (2001). Quantifying turfgrass cover using digital image analysis.Crop Sci.4118841888. 10.2135/cropsci2001.1884

  • 26

    SaberioonM. M.AminM. S. M.AnuarA. R.GholizadehA.WayayokA.Khairunniza-BejoS. (2014). Assessment of rice leaf chlorophyll content using visible bands at different growth stages at both the leaf and canopy scale.Int. J. Appl. Earth Obs. Geoinf.323545. 10.1016/J.JAG.2014.03.018

  • 27

    SommerH.LulisT. T.KaminskiJ. E. (2017). “Assessment of turfgrass health using aerial images from small unmanned aircraft,” in Proceedings of the ASA, CSSA and SSSA Annual Meeting, Tampa, FL.

  • 28

    SugiuraR.TsudaS.TamiyaS.ItohA.NishiwakiK.MurakamiN.et al (2016). Field phenotyping system for the assessment of potato late blight resistance using RGB imagery from an unmanned aerial vehicle.Biosyst. Eng.148110. 10.1016/J.BIOSYSTEMSENG.2016.04.010

  • 29

    SullivanD.ZhangJ.KowalewskiA. R.PeakeJ. B.AndersonW. F.WaltzF. C.et al (2017). Evaluating hybrid bermudagrass using spectral reflectance under different mowing heights and trinexapac-ethyl applications.Horttechnology274553. 10.21273/HORTTECH03436-16

  • 30

    Torres-SánchezJ.López-GranadosF.De CastroA. I.Peña-BarragánJ. M. (2013). Configuration and specifications of an unmanned aerial vehicle (UAV) for early site specific weed management (D Abbott, Ed.).PLoS One8:e58210. 10.1371/journal.pone.0058210

  • 31

    Torres-SánchezJ.PeñaJ. M.de CastroA. I.López-GranadosF. (2014). Multi-temporal mapping of the vegetation fraction in early-season wheat fields using images from UAV.Comput. Electron. Agric.103104113. 10.1016/J.COMPAG.2014.02.009

  • 32

    WherleyB.ChandraA.GenovesiA.KearnsM.PepperT.ThomasJ.et al (2013). Developmental response of St. Augustinegrass cultivars and experimental lines in moderate and heavy shade.HortScience4810471051. 10.21273/HORTSCI.48.8.1047

  • 33

    WickhamH. (2016). ggplot2: Elegant Graphics for Data Analysis.New York, NY: Springer-Verlag. 10.1007/978-3-319-24277-4

  • 34

    WoebbeckeD. M.MeyerG. E.VonBargenK.MortensenD. A. (1995). Color indices for weed identification under various soil, residue, and lighting conditions.Trans. ASAE38259269. 10.13031/2013.27838

  • 35

    XiangH.TianL. (2011). Development of a low-cost agricultural remote sensing system based on an autonomous unmanned aerial vehicle (UAV).Biosyst. Eng.108174190. 10.1016/J.BIOSYSTEMSENG.2010.11.010

  • 36

    XiongX.BellG. E.SolieJ. B.SmithM. W.MartinB. (2007). Bermudagrass seasonal responses to nitrogen fertilization and irrigation detected using optical sensing.Crop Sci.4716031610. 10.2135/cropsci2006.06.0400

  • 37

    Zarco-TejadaP. J.González-DugoV.WilliamsL. E.SuárezL.BerniJ. A. J.GoldhamerD.et al (2013). A PRI-based water stress index combining structural and chlorophyll effects: assessment using diurnal narrow-band airborne imagery and the CWSI thermal index.Remote Sens. Environ.1383850. 10.1016/J.RSE.2013.07.024

  • 38

    ZhangJ.PoudelB.KenworthyK. E.UnruhJ. B.RowlandD. L.EricksonJ. E.et al (2018). Drought responses of above-ground and below-ground characteristics in warm-season turfgrass.J. Agron. Crop Sci.205112. 10.1111/jac.12301

  • 39

    ZhouX.ZhengH. B.XuX. Q.HeJ. Y.GeX. K.YaoX.et al (2017). Predicting grain yield in rice using multi-temporal vegetation indices from UAV-based multispectral and digital imagery.ISPRS J. Photogramm. Remote Sens.130246255. 10.1016/J.ISPRSJPRS.2017.05.003

Summary

Keywords

bermudagrass, zoysiagrass, drought, normalized difference vegetation index, visible atmospherically resistant index, remote sensing

Citation

Zhang J, Virk S, Porter W, Kenworthy K, Sullivan D and Schwartz B (2019) Applications of Unmanned Aerial Vehicle Based Imagery in Turfgrass Field Trials. Front. Plant Sci. 10:279. doi: 10.3389/fpls.2019.00279

Received

08 October 2018

Accepted

20 February 2019

Published

15 March 2019

Volume

10 - 2019

Edited by

Yann Guédon, Centre de Coopération Internationale en Recherche Agronomique pour le Développement (CIRAD), France

Reviewed by

Jorge Torres Sánchez, Instituto de Agricultura Sostenible (IAS), Spain; Grégory Beurier, Centre de Coopération Internationale en Recherche Agronomique pour le Développement (CIRAD), France; Dale Bremer, Kansas State University, United States

Updates

Copyright

*Correspondence: Brian Schwartz,

This article was submitted to Technical Advances in Plant Science, a section of the journal Frontiers in Plant Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics