# Estimates of Maize Plant Density from UAV RGB Images Using Faster-RCNN Detection Model: Impact of the Spatial Resolution

K. Velumani, R. Lopez-Lozano, S. Madec, W. Guo, J. Gillet, A. Comar, F.

Baret

► To cite this version:

K. Velumani, R. Lopez-Lozano, S. Madec, W. Guo, J. Gillet, et al.

HAL Id: hal-03769738

https://hal.inrae.fr/hal-03769738

Submitted on 5 Sep 2022

**HAL** is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Research Article

# Estimates of Maize Plant Density from UAV RGB Images Using Faster-RCNN Detection Model: Impact of the Spatial Resolution

K. Velumani ID<fnref n="1" />,<fnref n="1" />,<fnref n="2" /> R. Lopez-Lozano ID<fnref n="2" />,<fnref n="2" /> S. Madec ID<fnref n="3" />,<fnref n="3" /> W. Guo ID<fnref n="4" />,<fnref n="4" /> J. Gillet,<fnref n="1" /> A. Comar ID<fnref n="1" />,<fnref n="1" /> and F. Baret ID<fnref n="2" />

Correspondence should be addressed to K. Velumani; kaaviya.velumani@inrae.fr and R. Lopez-Lozano; raul.lopez-lozano@inrae.fr

Received 16 April 2021; Accepted 2 July 2021; Published 21 August 2021

Copyright © 2021 K. Velumani et al. Exclusive Licensee Nanjing Agricultural University. Distributed under a Creative Commons Attribution License (CC BY 4.0).

Early-stage plant density is an essential trait that determines the fate of a genotype under given environmental conditions and management practices. The use of RGB images taken from UAVs may replace the traditional visual counting in fields with improved throughput, accuracy, and access to plant localization. However, high-resolution images are required to detect the small plants present at the early stages. This study explores the impact of image ground sampling distance (GSD) on the performances of maize plant detection at three-to-five leaves stage using Faster-RCNN object detection algorithm. Data collected at high resolution (GSD ≈ 0.3 cm) over six contrasted sites were used for model training. Two additional sites with images acquired both at high and low (GSD ≈ 0.6 cm) resolutions were used to evaluate the model performances. Results show that Faster-RCNN achieved very good plant detection and counting (rRMSE = 0.08) performances when native high-resolution images are used both for training and validation. Similarly, good performances were observed (rRMSE = 0.11) when the model is trained over synthetic low-resolution images obtained by downsampling the native training high-resolution images and applied to the synthetic low-resolution validation images. Conversely, poor performances are obtained when the model is trained on a given spatial resolution and applied to another spatial resolution. Training on a mix of high- and low-resolution images allows to get very good performances on the native high-resolution (rRMSE = 0.06) and synthetic low-resolution (rRMSE = 0.10) images. However, very low performances are still observed over the native low-resolution images (rRMSE = 0.48), mainly due to the poor quality of the native low-resolution images. Finally, an advanced super resolution method based on GAN (generative adversarial network) that introduces additional textural information derived from the native high-resolution images was applied to the native low-resolution validation images. Results show some significant improvement (rRMSE = 0.22) compared to bicubic upsampling approach, while still far below the performances achieved over the native high-resolution images.

## 1. Introduction

Plant density at emergence is an essential trait for crops since it is the first yield component that determines the fate of a genotype under given environmental conditions and management practices [1-5]. Competition between plants within the canopy depends on the sowing pattern and its understanding requires reliable observations of the plant localization and density [6-9]. An accurate estimation of actual plant density

is also necessary to evaluate the seed vigor by linking the emergence rate to the environmental factors [10-13].

Maize plant density is measured by visual counting in the field. However, this method is labor intensive, time consuming, and prone to sampling errors. Several higher throughput methods based on optical imagery have been developed in the last twenty years. This was permitted by the technological advances with the increasing availability of small, light, and affordable high spatial resolution cameras and autonomous

<footnote n="1">¹Hiphen SAS, 120 Rue Jean Dausset, Agroparc, Bâtiment Technicité, 84140 Avignon, France</footnote>

<footnote n="2">²INRAE, UMR EMMAH, UMT CAPTE, 228 Route de l'Aérodrome, Domaine Saint Paul-Site Agroparc CS 40509, 84914 Avignon Cedex 9, France</footnote>

<footnote n="3">³Arvalis, 228, Route de l'Aérodrome-CS 40509, 84914 Avignon Cedex 9, France</footnote>

<footnote n="4">⁴International Field Phenomics Research Laboratory, Institute for Sustainable Agro-Ecosystem Services, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan</footnote>
vehicles. Unmanned ground vehicles (UGV) provide access to detailed phenotypic traits [14–16] while being generally expensive and associated with throughputs of the order of few hundreds of microplots per hour. Conversely, unmanned aerial vehicles (UAV) are very affordable with higher acquisition throughput than UGVs. When carrying very high-resolution cameras, they can access potentially several traits [17, 18] including plant density [19, 20].

Image interpretation methods used to estimate plant density can be classified into three main categories. The first one is based on machine learning where the plant density measured over a small set of sampling area is related to other canopy level descriptors including vegetation indices derived from RGB and multispectral data [21–23]. However, this type of method may lead to significant errors due to the lack of representativeness of the training dataset as well as the effect of possible confounding factors including changes in background properties or plant architecture under genetic control. The second category of methods is based on standard computer vision techniques, where the image is first binarized to identify the green objects that are then classified into plants according to the geometrical features defined by the operator (e.g. [24, 25]). The last category of methods widely used now is based on deep learning algorithms for automatic object detection [26–28].

The main advantage of deep learning methods is their ability to automatically extract low-level features from the images to identify the targeted objects. Although deep learning methods appear very promising, their generalization capacity is determined by the volume and diversity of the training dataset [29]. While large collections of images can now be easily acquired, labeling the images used to train the deep models represents a significant effort that is the main limiting factor to build very large training datasets. Few international initiatives have been proposed to share massive labeled datasets that will contribute to maximize the performances of deep learning models [30–34], with however questions regarding the consistency of the acquisition conditions and particularly the ground sampling distance (GSD).

The use of UAV images for plant detection at early stages introduces important requirements on image resolution, as deep learning algorithms are sensitive to object scales with the identification of small objects being very challenging [35, 36]. For a given camera, low-altitude flights are therefore preferred to get the desired GSD. However, low-altitude flights decrease the acquisition throughput because of a reduced camera swath forcing to complete more tracks to cover the same experiment and require additionally to slow down the flying speed to reduce motion blur. An optimal altitude should therefore be selected to compromise between the acquisition throughput and the image GSD. Previous studies reporting early-stage maize plant detection from UAVs from deep learning methods did not addressed specifically this important scaling issue [20, 26, 27]. One way to address this scaling issue is to transform the low-resolution images into higher resolution ones using super resolution techniques. Dai et al. [37] have demonstrated the efficiency of super resolution techniques to enhance segmentation and edge detection.

Later, Fromm et al. [38] and Magoulianitis et al. [39] showed improvements in object detection performances when using the super resolution methods. The more advanced super resolution techniques use deep convolutional networks trained over paired high- and low-resolution images [40–42]. Since the construction of a real-world paired high- and low-resolution dataset is a complicated task, the high-resolution images are often degraded using a bicubic kernel or less frequently using Gaussian noise to constitute the low-resolution images [43]. However, more recent studies have shown the drawbacks of the bicubic down-sampling approaches as it smoothens sensor noise and other compression artifacts, thus failing to generalize while applied to real world images [41]. More recent studies propose the use of unsupervised domain translation techniques to generate realistic paired datasets for training the super resolution networks [44].

We propose here to explore the impact of image GSD on the performances of maize plant detection at stages from three to five leaves using deep learning methods. More specifically, three specific objectives are targeted: (1) to assess the accuracy and robustness of deep learning algorithms for detecting maize plants with high-resolution images used both, in the training and validation datasets; (2) to study the ability of these algorithms to generalize in the resolution domain, i.e. when applied to images with higher and lower resolution compared to the training dataset; and (3) to evaluate the efficiency of data augmentation and preparation techniques in the resolution domain to improve the detection performances. Special emphasis was put here on assessing the contribution of two contrasting methods to upsample low-resolution images: a simple bicubic upsampling algorithm and a more advanced super resolution model based on GAN (generative adversarial network) that introduces additional textural information. Data collected over several sites across France with UAV flights completed at several altitudes providing a range of GSDs were used.

## 2. Materials and Methods

2.1. *Study Sites.* This study was conducted over 8 sites corresponding to different field phenotyping platforms distributed across the west of France and sampled from 2016 to 2019 (Figure 1). The list of sites and their geographic coordinates are given in Table 1. Each platform included different maize microplots with size 20 to 40 square meters. Depending on the experimental design of the platform, the microplots were sown with two to seven rows of maize of different cultivars and row spacing varying from 30 to 110 cm. The sowing dates were always between mid-April and mid-May.

2.2. *Data Acquisition and Processing.* UAV flights were carried out on the eight sites approximately one month after the sowing date, between mid-May and mid-June (Table 1). Maize plants were in between three and five leaf stage, ensuring that there is almost no overlap among individual plants from near nadir viewing. The microplots were weeded and consisted of only maize plants. Three different RGB cameras were used for the data acquisition: Sony Alpha (ILCE-6000)
<figure><img src="image_2.png" /><figcaption>FIGURE 1: Location of the study sites with example extracts of the maize microplots acquired from UAV. (a) A map displaying the location of the eight maize phenotyping platforms located in the west of France used in this study. (b) An illustration of the bounding boxes drawn around the maize plants. The examples shown are from the Tartas site (GSD = 0.27 cm) (A) and Tartas site (GSD = 0.63 cm) (B).</figcaption></figure>

with a focal length of 30 mm, DJI X7 (FC6540) with focal lengths of 24 mm and 30 mm, and the default camera with DJI Mavic 2 pro (L1D-20c) with a focal length of 10.26 mm mounted on AltiGator Mikrokopter (Belgium) and DJI Mavic 2 pro (China). To georeference the images, ground control points (GCPs) were evenly distributed around the sites and their geographic coordinates were registered using a Real-Time Kinematic GPS.

The flights were conducted at an altitude above the ground ranging between 15 and 22 meters, providing a ground sampling distance (GSD) between 0.27 and 0.35 cm (Table 1). For the Tartas and Selommes sites, an additional flight was done at a higher altitude on the same day providing a GSD between 0.63 and 0.66 cm.

The flights were planned with a lateral and front overlap of 60/80% between individual images. Each dataset was processed using PhotoScan Professional (Agisoft LLC, Russia) to align the overlapping images by automatic tie point matching, optimize the aligned camera positions, and finally georeference the results using the GCPs. The steps followed are similar to the data processing detailed by Madec et al. [15]. Once orthorectified, the multiple instances of the microplot present in the overlapping images were extracted using Phenoscript, software developed within the CAPTE research unit. Phenoscript allows to select, among the individual images available for each microplots, those with full coverage of the microplot, minimum blur, and view direction closer to the nadir one. Only these images were used in this study.

## 2.3. Manual Labeling of Individual Plants.
From each site, the microplots were labeled with an offline tool, LabelImg [45]: bounding boxes around each maize plant were interactively drawn (Figure 1(b)) and saved in the Pascal VOC format as XML files. The available sites (Table 1) were divided into three groups: (1) the first group ($T^h$) composed of six sites was used to train the plant detection models. It includes a total of 202 microplots corresponding to 19,841 plants. (2) The second group ($V^h$) corresponding to the Tartas and Selommes with low-altitude flights was used to evaluate the model performance at high resolution. It includes a total of 36 microplots corresponding to 3256 plants. (3) The third group ($V^l$) corresponds to the high-altitude flights in Tartas and Selommes was used to evaluate the model performance at low resolution. It includes a total of 36 microplots corresponding to 3256 plants. An example of images extracted from the three groups is shown in Figure 2.

## 2.4. The Faster-RCNN Object Detection Model.
Faster-RCNN [46], a convolutional neural network designed for object detection, was selected to identify maize plants in the image. Besides its wide popularity outside the plant phenotyping community, Faster-RCNN has also been proved to be suitable for various plant and plant-organ detection tasks [47–49]. We used the implementation of Faster-RCNN in the open-source MMDetection Toolbox [50], written in PyTorch, with pretrained weights on ImageNet. The Faster-
<table><caption>TABLE 1: The dataset used for the training and validation of the object detection models are listed here.</caption><thead><tr><th>Dataset #</th><th>Site name</th><th>Latitude (°)</th><th>Longitude (°)</th><th>Acquisition date</th><th>Camera</th><th>Flight altitude (m)</th><th>Focal length (mm)</th><th>GSD (cm)</th><th>No. microplots labeled</th><th>No. plants labeled</th><th>Average plant size (pixels)</th></tr></thead><tbody><tr><td>1</td><td>Castetis</td><td>43.46</td><td>-0.71</td><td>06-06-2019</td><td>FC6540</td><td>26</td><td>24</td><td>0.35</td><td>20</td><td>2239</td><td>939</td></tr><tr><td>2</td><td>St. Hermine</td><td>46.54</td><td>-1.06</td><td>23-05-2019</td><td>FC6540</td><td>26</td><td>35</td><td>0.33*</td><td>20</td><td>2674</td><td>806</td></tr><tr><td>3</td><td>Nerac</td><td>44.16</td><td>0.3</td><td>01-06-2017</td><td>ILCE-6000</td><td>25</td><td>30</td><td>0.32</td><td>44</td><td>3338</td><td>2253</td></tr><tr><td>4</td><td>Thenay</td><td>47.38</td><td>1.28</td><td>18-05-2018</td><td>ILCE-6000</td><td>22</td><td>30</td><td>0.25</td><td>72</td><td>7454</td><td>1505</td></tr><tr><td>5</td><td>Villedieu</td><td>47.88</td><td>1.53</td><td>28-06-2016</td><td>n/a</td><td>n/a</td><td>n/a</td><td>0.27</td><td>26</td><td>2390</td><td>2159</td></tr><tr><td>6</td><td>Blois</td><td>47.56</td><td>1.32</td><td>18-05-2018</td><td>ILCE-6000</td><td>25</td><td>30</td><td>0.33</td><td>20</td><td>1746</td><td>1419</td></tr><tr><td>7</td><td>Tartas</td><td>43.80</td><td>-0.79</td><td>08-06-2019</td><td>FC6540</td><td>20</td><td>24</td><td>0.32</td><td>22</td><td>2151</td><td>1336</td></tr><tr><td>8</td><td>Selommes</td><td>47.76</td><td>1.19</td><td>17-05-2019</td><td>LID-20c</td><td>16.2</td><td>10.26</td><td>0.27</td><td>14</td><td>1105</td><td>891</td></tr><tr><td>9</td><td>Tartas</td><td>43.80</td><td>-0.79</td><td>08-06-2019</td><td>FC6540</td><td>40</td><td>24</td><td>0.63</td><td>24</td><td>2151</td><td>437</td></tr><tr><td>10</td><td>Selommes</td><td>47.76</td><td>1.19</td><td>17-05-2019</td><td>LID-20c</td><td>30</td><td>10.26</td><td>0.66</td><td>14</td><td>1105</td><td>156</td></tr></tbody></table>

T<sup>h</sup> is the training high-resolution dataset, V<sup>h</sup> is the validation high-resolution dataset, and V<sup>l</sup> is the validation low-resolution dataset. *For this site, the microplot extracts were resampled to GSD = 0.25 cm before annotating.
<figure><img src="image_3.png" /><figcaption>FIGURE 2: Examples of maize plants extracted from the in the eight sites used in this study. The image titles indicate the location of the sites. *T*<sup>*h*</sup>, *V*<sup>*h*</sup>, and *V*<sup>*l*</sup> are the training high-resolution dataset, validation high-resolution dataset, and the validation low-resolution dataset, respectively.</figcaption></figure>

RCNN model with a ResNet50 backbone was trained for 12 epochs with a batch size of 2. The weights were optimized using an SGD optimizer (stochastic gradient descent) with a learning rate of 0.02. For the model training, ten patches of 512 × 512 pixels were randomly extracted from each microplot in the training sites. Standard data augmentation strategies such as rotate, flip, scale, brightness/contrast, and jpeg compression were applied.

## 2.5. *Experimental Plan.*
To evaluate the effect of the resolution on the reliability of maize plant detection, we compared Faster-RCNN performances over training and validation datasets made of images of high (GSD ≈ 0.30 cm) and low (GSD ≈ 0.60 cm) resolution. Three training datasets built from *T*<sup>*h*</sup> (Table 1) were considered: (1) the original *T*<sup>*h*</sup> dataset with around 0.32 cm GSD; (2) a dataset, *T*<sub>gm</sub><sup>*h→l*</sup> where the images from *T*<sup>*h*</sup> were downsampled to 0.64 cm GSD using a Gaussian filter and motion blur that mimics the actual low-resolution imagery acquired at higher altitude as described later (Section 2.6.1); and (3) a dataset, where the original *T*<sup>*h*</sup> high-resolution dataset was merged with its low-resolution transform, *T*<sub>gm</sub><sup>*h→l*</sup>. This *T*<sup>*h*</sup> + *T*<sub>gm</sub><sup>*h→l*</sup> is expected to provide robustness of the model towards changes in GSD. Note that we did not investigate the training with the native low-resolution images because the labeling of low-resolution images is often difficult because plants are not easy to identify visually and to draw accurately the corresponding bounding box. Further, only two flights were available at the high altitudes (Table 1) that were reserved for the validation. A specific model was trained over each of the three training datasets considered (Table 2) and then evaluated over independent high- and low-resolution validation datasets.

We considered three validation datasets for the high-resolution images: (1) the native high-resolution validation dataset, *V*<sup>*h*</sup> acquired at low altitude with GSD around

0.30 cm (Table 1); (2) a synthetic high-resolution dataset of GSD around 0.30 cm obtained by upsampling the native low-resolution dataset, acquired at high altitude, using a bicubic interpolation algorithm as described in Section 2.6.2, and it will be called *V*<sub>bc</sub><sup>*l→h*</sup>; and (3) a synthetic high-resolution dataset, *V*<sub>sr</sub><sup>*l→h*</sup>, obtained by applying a super resolution algorithm (see Section 2.6.3) to the native low-resolution dataset *V*<sup>*l*</sup> and resulting in images with a GSD around 0.30 cm. Finally two low-resolution datasets will be also considered: (1) the native low-resolution validation dataset, *V*<sup>*l*</sup> (Table 1), with a GSD around 0.60 cm and (2) a synthetic low-resolution dataset, *V*<sub>gm</sub><sup>*h→l*</sup>, obtained by applying a Gaussian filter to downsample (see Section 2.6.1) the original high-resolution dataset, *V*<sup>*h*</sup>, and get a GSD around 0.60 cm.

## 2.6. *Methods for Image Up- and Downsampling*

### 2.6.1. *Gaussian Filter Downsampling*
To create the synthetic low-resolution datasets *T*<sub>gm</sub><sup>*h→l*</sup> and *V*<sub>gm</sub><sup>*h→l*</sup>, a Gaussian filter with a sigma = 0.63 and a window size = 9 followed by a motion blur with a kernel size = 3 and angle = 45 were applied to downsample the native high-resolution datasets *T*<sup>*h*</sup> and *V*<sup>*h*</sup> by a factor of 2. This solution was preferred to the commonly used bicubic downsampling method because it provides low-resolution images more similar to the native low-resolution UAV images (Figure 3). This was confirmed by comparing the image variance over the Selommes and Tartas sites where both native high- and low-resolution images were available: the variance of the *V*<sub>gm</sub><sup>*h→l*</sup> was closer to that of *V*<sup>*l*</sup> whereas the bicubic downsampled dataset had a larger variance corresponding to sharper images. This is consistent with [38, 51] who used the same method to realistically downsample high-resolution images.
<table><caption>TABLE 2: Description of the training and validation datasets.</caption><thead><tr><th></th><th>Dataset name</th><th>No. microplots</th><th>No. plants</th><th>Comment</th></tr></thead><tbody><tr><td rowspan="3">Training</td><td>T<sup>h</sup></td><td>202</td><td>19,841</td><td>Native high-resolution training dataset</td></tr><tr><td>T<sup>h→l</sup><sub>gm</sub></td><td>202</td><td>19,841</td><td>Downsampling T<sup>h</sup> with Gaussian filter and motion blur</td></tr><tr><td>T<sup>h</sup>+T<sup>h→l</sup><sub>gm</sub></td><td>404</td><td>39,682</td><td>Merging T<sup>h</sup> and T<sup>h→l</sup></td></tr><tr><td rowspan="5">Validation</td><td>V<sup>h</sup></td><td>36</td><td>3256</td><td>Native high-resolution validation dataset</td></tr><tr><td>V<sup>l→h</sup><sub>bc</sub></td><td>36</td><td>3256</td><td>Upsampling V<sup>l</sup> with bicubic algorithm</td></tr><tr><td>V<sup>l→h</sup><sub>sr</sub></td><td>36</td><td>3256</td><td>Upsampling V<sup>l</sup> with Cycle-ESRGAN super resolution</td></tr><tr><td>V<sup>l</sup></td><td>36</td><td>3256</td><td>Native low-resolution validation dataset</td></tr><tr><td>V<sup>h→l</sup><sub>gm</sub></td><td>36</td><td>3256</td><td>Downsampling V<sup>h</sup> with gaussian filter and motion blur</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>FIGURE 3: Visual comparison of the extract of the same plant from the Tartas site between different versions of low resolution. Native low resolution, synthetic low resolution from bicubic downsampling, and synthetic low resolution from Gaussian downsampling (sigma = 0.63, window = 9) followed by a motion blur (kernel size = 3 and angle = 45).</figcaption></figure>

### 2.6.2. Bicubic Upsampling.
The bicubic interpolation algorithm was used to generate $V_{bc}^{l→h}$ by upsampling the native low-resolution UAV images, $V^l$. The bicubic interpolation available within Pillow, the Python Imaging Library [52], was used to resample the images.

### 2.6.3. Super Resolution Images Derived from Cycle-ESRGAN.
The super resolution (SR) is an advanced technique that artificially enhances the textural information while upsampling images. We used a SR model inspired from [53]. It is a two-stage network composed of a CycleGAN network that generates synthetic paired data and a ESRGAN network capable of image upsampling. The CycleGAN [54] performs unsupervised domain mapping between the native low-resolution and bicubic downsampled low-resolution domains. Thus, for any given input image, CycleGAN is trained to add realistic image noise typical of low-resolution images. The ESRGAN-type super resolution network [42] upsamples by a factor of two low-resolution images.

During the training phase, the CycleGAN stage of the network was trained to generate “real-world” noise from an unpaired dataset of native low-resolution and bicubically downsampled images. The second stage of the network consisting of the ESRGAN was then trained using a paired dataset of native high-resolution image and a “realistic” low-resolution image generated by the CycleGAN (Figure 4). The CycleGAN stage of the network was initially

trained for a few epochs following which the two stages (CycleGAN+ESRGAN) were trained together simultaneously. It should be noted that during inference, only the ESRGAN stage of the network would be activated. Hence, for a given input image, the ESRGAN network would upsample the input by a factor of 2. The training parameters and losses reported by Han et al. [53] were used for the model training. The model weights were initialized over the Div2k dataset [55] and finetuned on the UAV dataset detailed below. The Cycle-ESRGAN network was implemented using Keras [56] deep learning library in Python. The codes will be made available on Github at the following link: https://github.com/kaaviyave/Cycle-ESRGAN.

A dedicated training dataset for the super resolution network was prepared using UAV imagery belonging to the following two domains:

(i) Native high-resolution domain: 2234 microplot extractions from four sites with an average GSD of less than 0.33 cm. Some of the sites belonging to the $T^h$ dataset was used as a part of the training

(ii) Native low-resolution domain: 1713 microplot extractions from three sites with an average GSD of 0.46 cm per site

None of the validation sites ($V^l$ and $V^h$ in Table 1) were used in the training of the super resolution model. The
<figure><img src="image_5.png" /><figcaption>FIGURE 4: An illustration of the super resolution pipeline. $I^h$ and $I^l$ are the native high- and low-resolution images, respectively. $I_{bc}^{h→l}$ and $I_{cyclegan}^{h→l}$ are the synthetic low-resolution images prepared from the native high-resolution images by bicubic downsampling and by the CycleGAN network, respectively. $I_{sr}^{l→h}$ is the synthetic high-resolution image generated by the super resolution network.</figcaption></figure>

synthetic downsampled dataset used to train the CycleGAN was prepared by bicubic downsampling the native high-resolution domain by a factor of 2. The images were split into patches of size 256 × 256 pixels for the high-resolution domain and into 128 × 128 pixels for the low-resolution domain.

## 2.7. Evaluation Metrics.
In this study, the average precision (AP), root mean squared error (RMSE), and accuracy will be utilized for the evaluation of the Faster-RCNN models for the purpose of maize plant detection and counting.

AP is a frequently used metric for the evaluation of object detection models and can be considered the area under the precision-recall curve.

$$ 
\begin{aligned}
\text{Precision} &= \frac{\text{TP}}{\text{TP} + \text{FP}}, \\
\text{Recall} &= \frac{\text{TP}}{\text{TP} + \text{FN}},
\end{aligned}
\qquad (1) $$

where TP is the number of true positive, FP is the number of false positive, and the FN is the number of false negative. For the calculation of AP, a predicted bounding box is considered: true positive (TP) if its intersection area over union area (IoU) with the corresponding labeled bounding box is larger than a given threshold. Depending on the objective of the study, different variations exist in the AP metric calculation and the choice of IOU threshold used to qualify a predicted

bounding box as TP. After considering several IoU threshold values, we decided to use an IoU threshold of 0.25 to compute AP. This will be later justified. The Python COCO API was used for the calculation of the AP metric [57].

Accuracy evaluates the model's performance by calculating the ratio of correctly identified plants to all the predictions made by the model. A predicted bounding box is considered true positive if it has a confidence score of more than 0.5 and an IoU threshold of 0.25. Accuracy is then calculated as

$$ \mathrm{Ac} = \frac{\mathrm{TP}}{\mathrm{TP} + \mathrm{FP} + \mathrm{FN}}. \qquad (2) $$

The relative root mean square error (rRMSE) between the number of labeled and detected plants across all images belonging to the same dataset:

$$ \mathrm{rRMSE} = \frac{\sqrt{\sum_{i}^{n} (P_{o,i} - P_{p,i})^2 / n}}{\bar{P}_{o,1}}, \qquad (3) $$

where $P_{o,i}$ is the number of plants labeled on image and $P_{p,i}$ is the number of images predicted by the CNN (confidence score > 0.5 and an IoU > 0.25) and $\bar{P}_{o,1}$ is the average number of labeled plants per image.
<figure><img src="image_6.png" /><figcaption>FIGURE 5: Results of the model trained on the native HR dataset, *T*<sup>*h*</sup>, and applied to the HR validation dataset, *V*<sup>*h*</sup>. (a) Elements of the confusion matrix—true positives (TP), false negatives (FN), and false positives (FP) for the Selommes and Tartas sites. (b) An example of false positive and false negative observed in the two validation sites. The ground truth bounding boxes are shown in green, and the predicted bounding boxes are shown in blue. The green text indicates the IoU of the predicted bounding box with the ground truth, and the blue text indicates the confidence score of the predictions.</figcaption></figure>

**3. Results and Discussion**

3.1. *Faster-RCNN Detects Plants with High Accuracy at High Spatial Resolution.* Very good performances (rRMSE = 0.08; Ac = 0.88, AP = 0.95) are achieved when the model is trained over the high-resolution images (*T*<sup>*h*</sup>) and applied on high-resolution images taken on independent sites (*V*<sup>*h*</sup>). The good performances are explained by the high rate of true positives (Figure 5(a)). However, the detector performs slightly differently on the two sites used for the validation: in Selommes, an overdetection (false positives, FP) is observed for a small number of plants, when the detector splits a plant into two different objects (Figure 5(b)). Conversely, in the Tartas site, some underdetection (false negatives, FN) is observed, with a small number of undetected plants (Figure 5).

A detailed analysis of the precision-recall curves for the configuration [*T*^h, *V*^h] at different IoU (Figure 6) shows a drastic degradation of the detector performances when the IoU is higher than 0.3. This indicates that the model is not accurate when determining the exact dimensions of maize plants. This is partly explained by the difficulty of separating the green from the ground in the shadowed parts of the images. As a consequence, some shaded leaves are excluded from the bounding boxes proposed by the detector and, conversely, some shadowed grounds are wrongly included in the bounding boxes proposed (Figure 5(b)). Further, when a single plant is split into two separate objects by the detector, the resulting bounding boxes are obviously smaller than the corresponding plant (Figure 5(b)). As a consequence, we proposed to use an IoU threshold of 0.25 to evaluate the model

<figure><img src="image_7.png" /><figcaption>FIGURE 6: Precision-recall curve at different IoU thresholds for the plant detection model trained and applied with high-resolution images (*T*^h, *V*^h).</figcaption></figure>

performance to better account for the smaller size of the
detected bounding boxes. This contrasts from most object
detection applications where an IoU threshold of 0.5 or
<table><caption>TABLE 3: Comparison of the performance of the Faster-RCNN models trained and validated over datasets with different resolutions.</caption><thead><tr><th rowspan="2"></th><th rowspan="2"></th><th colspan="4">High resolution (V<sup>h</sup>)</th><th colspan="3">Validation</th></tr><tr><th>rRMSE</th><th>Ac</th><th>AP</th><th>rRMSE</th><th>Ac</th><th>AP</th></tr></thead><tbody><tr><td rowspan="2">Training</td><td>T<sup>h</sup></td><td>0.08</td><td>0.88</td><td>0.95</td><td>0.48</td><td>0.54</td><td>0.64</td></tr><tr><td>T<sup>h→l</sup><sub>gm</sub></td><td>0.52</td><td>0.56</td><td>0.71</td><td>0.29</td><td>0.76</td><td>0.81</td></tr></tbody></table>

<figure><img src="image_8.png" /><figcaption>FIGURE 7: Results of maize plant detection when trained and evaluated across different resolution domains. T<sup>h</sup>: native high-resolution training dataset; T<sup>h→l</sup><sub>gm</sub>: low-resolution training dataset by downsampling T<sup>h</sup> using Gaussian motion blur; V<sup>h</sup>: native high-resolution validation dataset; V<sup>l</sup>: native low-resolution validation dataset; V<sup>l→h</sup><sub>bc</sub>: high-resolution dataset by upsampling V<sup>l</sup> using bicubic interpolation.</figcaption></figure>

0.75 is commonly used to evaluate the performance of the methods [58, 59]. The observed degradation of the model performance for IoU above 0.3 indicates that the method presented provides less accurate localization than in other object detection studies, including both, real world objects and phenotyping applications [49, 60, 61]. An inaccurate estimation of plant dimensions is not critical for those applications assessing germination or emergence rates and uniformity, where plant density is the targeted phenotypic trait. If the focus is to additionally assess the plant size in early developmental stages as well, mask-based RCNN [62, 63] could be used instead. In contrast to algorithms trained on rectangular regions like Faster-RCNN, mask-based algorithms have the potentials to more efficiently manage the shadow projected on the ground by plants, limiting therefore the possible confusion between shaded leaves and ground during the training. However, generating mask annotations is time consuming, increasing the effort needed to generate a diverse training dataset.

These results provide slightly better performances as those reported by David et al. [20] with Ac ≈ 0.8 and rRMSE ≈ 0.1 when using the “out-domain” approach as the one used in this study, i.e. when the training and validation

sites are completely independent. They used images with a spatial resolution around 0.3 cm as in our study. This is also consistent with the results of Karami et al. [26] who obtained an accuracy of 0.82 with a spatial resolution of around 1 cm. They used the anchor-free few shot learning (FSL) method which identifies and localizes the maize plants by estimating the central position. They claim that their method is a little sensitive to object size and thus to the spatial resolution of the images. The accuracy increases up to 0.89 when introducing few images from the validation sites in the training dataset. Kitano et al. [27] proposed a two-step method: they first segment the images using a CNN-based method and then count the segmented objects. They report an average rRMSE of 0.24 over a test dataset where many factors including image resolution vary (ranging from GSD ≈ 0.3 cm to 0.56 cm). They report that their method is sensitive to the size and density of the objects. In the following, we will further investigate the dependency of the performances to image resolution.

3.2. *The Faster-RCNN Model Is Sensitive to Image Resolution and Apparent Plant Size.* The performances of the model were evaluated when it is trained and validated over images
<figure><img src="image_9.png" /><figcaption>FIGURE 8: Effects of the hyperspecialization of Faster-RCNN trained with synthetic low-resolution images ($T_{gm}^{h→l}$) and applied to a high-resolution dataset ($V^l$). (a) Relationship between the size of the ground truth bounding boxes and the average number of predicted bounding boxes intersecting with them. (b) Example of overdetection of maize plants due to different object size. The ground truth bounding boxes are shown in green, and the predicted boxes are shown in blue. The green text indicates the IoU of the predicted boxes with the ground truth, and the blue text indicates the confidence score of the predictions.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>FIGURE 9: An example showing the same plants extracted from the exact same locations in two versions of the validation dataset: native LR ($V^l$) and the synthetic LR obtained from Gaussian downsampling ($V_{gm}^{h→l}$). The first and third columns show the raw images while the second and fourth columns show the detector predictions. The ground truth bounding boxes are shown in green, and the predicted bounding boxes are shown in blue. The green text indicates the IoU of the predicted box with the ground truth, and the blue text indicates the confidence score of the predictions.</figcaption></figure>

<table><caption>TABLE 4: Comparison of the performance of the Faster-RCNN models trained over the augmented data ($T^h + T_{gm}^{h→l}$) and validated over datasets with different resolutions.</caption><thead><tr><th></th><th>rRMSE</th><th>Ac</th><th>AP</th></tr></thead><tbody><tr><td>V<sup>h</sup></td><td>0.06</td><td>0.91</td><td>0.96</td></tr><tr><td>V<sup>h→l</sup><sub>gm</sub></td><td>0.11</td><td>0.92</td><td>0.91</td></tr><tr><td>V<sup>l</sup></td><td>0.48</td><td>0.64</td><td>0.72</td></tr></tbody></table>

<table><caption>TABLE 5: Comparison of the performance of the Faster-RCNN models trained over high-resolution images and applied to the native low-resolution images ($V^l$), the synthetic high-resolution images that is upsampled/transformed using either bicubic ($V_{bc}^{l→h}$) or super resolution ($V_{sr}^{l→h}$) techniques.</caption><thead><tr><th></th><th>rRMSE</th><th>Ac</th><th>AP</th></tr></thead><tbody><tr><td>V<sup>l</sup></td><td>0.58</td><td>0.54</td><td>0.64</td></tr><tr><td>V<sup>l→h</sup><sub>bc</sub></td><td>0.43</td><td>0.63</td><td>0.77</td></tr><tr><td>V<sup>l→h</sup><sub>sr</sub></td><td>0.22</td><td>0.80</td><td>0.85</td></tr></tbody></table>
<figure><img src="image_11.png" /><figcaption>FIGURE 10: A comparison between the performance of the Faster-RCNN model trained on the HR dataset, $T^h$, and applied to the synthetic high-resolution datasets. (a) Model trained on $T^h$ and evaluated on the synthetic HR dataset $V_{bc}^{l→h}$, from the bicubic upsampling. (b) Model trained on $T^h$ applied to the synthetic HR dataset $V_{sr}^{l→h}$, from the super resolution technique.</figcaption></figure>

with different resolution. When Faster-RCNN is trained on the high-resolution domain ($T^h$) and applied to a dataset with low resolution ($V^l$), both AP and Ac decrease almost by 30% (Table 3) compared to the results where the model is trained and applied over high-resolution images. The rate of true positive drops because of the drastic increase of false negatives indicating a high rate of misdetection (Figure 7, $[T^h, V^l]$). This degradation of the detection performances impacts highly the rRMSE that increases up to 0.48. This indicates that the model is sensitive to the resolution of the images. We further investigated if this was linked to the apparent size of the plants and therefore upsampled the validation low-resolution images with a bicubic interpolation method ($V_{bc}^{l→h}$) to get plants with the same size as in the native high-resolution images ($V^h$). Results show that Ac increases from 0.54 to 0.63 and AP from 0.64 to 0.77. However, because of the high imbalance between FN and FP (Figure 7, $[T^h, V_{bc}^{l→h}]$), the counting performances remain poor with rRMSE = 0.49.

When the model is trained over simulated low-resolution images ($T_{gm}^{h→l}$), the detection and counting performances evaluated on high-resolution images ($V^h$) also degrades drastically (Table 3). The rate of true positive is relatively high, but the rate of false positive increases drastically (Figure 7 ($T_{gm}^{h→l}, V^h$)). We observe that the average number of predicted bounding boxes overlapping each labeled box increases linearly with its size (Figure 8). For example, the model identifies on average two plants inside plants larger than 4000 pixels. The imbalance between FN and FP explains the very poor counting performances with rRMSE = 0.52 (Table 3). This result confirms the importance to keep consistent the resolution

and plant size between the training and the application datasets since Faster-RCNN tends to identify objects that have a similar size to the objects used during the training.

We thus evaluated whether data augmentation may improve the performances on the low-resolution images ($V^l$): the Faster-RCNN model trained on the simulated low-resolution images ($T_{gm}^{h→l}$) shows improved detection performances as compared to the training over the native high-resolution images (Table 3) with a decrease of the rRMSE down to 0.29 (Table 3). When this model trained with synthetic low-resolution images ($T_{gm}^{h→l}$) is applied to a dataset downscaled to a similar resolution ($V_{gm}^{h→l}$), the performances improve dramatically with Ac increasing from 0.56 to 0.89 and AP from 0.71 to 0.90 while the rRMSE drops to 0.10. However, when this model trained with synthetic low-resolution images ($T_{gm}^{h→l}$) is applied to the native low-resolution images ($V^l$), moderate detection performances are observed which degrades the counting estimates with rRMSE = 0.29 (Table 3).

The performances of the model trained over the synthetic low-resolution images ($T_{gm}^{h→l}$) are quite different when evaluated over the native images ($V^l$) or the synthetic ones ($V_{gm}^{h→l}$) with the latter yielding results almost comparable to the high-resolution configurations with AP = 0.90 (Table 3). This indicates that the low-resolution synthetic images contain enough information to detect accurately the maize plants. Conversely, the native low-resolution image, $V^l$, has probably lost part of the textural information. In addition, the model trained on the synthetic low-resolution images is not able to extract the remaining pertinent plant
<figure><img src="image_12.png" /><figcaption>FIGURE 11: Illustration of the performance of the Faster-RCNN model on the synthetic high-resolution and native low-resolution datasets. (a, c) Images belonging to three datasets: native low resolution $V^l$, bicubically upsampled $V_{bc}^{l→h}$, and finally the upsampling by super resolution technique $V_{sr}^{l→h}$. (b, d) The results predicted by the model trained on $T_{gm}^{h→l}$ applied to $V^l$ (first column) and the model trained on $T^h$ applied to the synthetic high-resolution datasets $V_{bc}^{l→h}$ and $V_{sr}^{l→h}$ (second and third columns). The ground truth bounding boxes are shown in green, and the predicted bounding boxes are shown in blue. The green text indicates the IoU of the predicted bounding box with the ground truth, and the blue text indicates the confidence score of the predictions.</figcaption></figure>

descriptors from the native low-resolution images. We can observe that the native low-resolution images contain less details as compared to the synthetic ones (Figure 9): some plants are almost not visible in the $V^l$ images, as the textural information vanishes and even the color of maize leaves can not be clearly distinguished from the soil background. This explains why the model was not able to detect the plants, even when it is trained with the synthetic low-resolution images ($T_{gm}^{h→l}$).

Contrary to vectors that operate at an almost constant height like ground vehicles [16, 64–66] or fixed cameras [67–70], camera settings (aperture, focus and integration time) in UAVs need to be adapted to the flight conditions,

especially flight altitude, to maximize image quality. Further, the jpg recording format of the images may also significantly impact image quality. Recording the images in raw format would thus improve the detection capability at the expense of increased data volume and sometimes image acquisition frequency.

### 3.3. Data Augmentation Makes the Model More Resistant to Changes in Image Resolution.
We finally investigated whether mixing high- and low-resolution images in the training dataset would make the model more resistant to changes in the image resolution. Results show that merging native high-resolution with synthetic low-resolution images
$(T^h + T_{gm}^{h \to l})$ provides (Table 4) performances similar to those observed when the model is trained only over high $(T^h)$ or synthetic low $(T_{gm}^{h \to l})$ and validated on the same resolution $(V^h)$ or $(V_{gm}^{h \to l})$ (Table 3). This proves that data augmentation could be a very efficient way to deal with images having different resolutions. Further, this model trained on augmented data $(T^h + T_{gm}^{h \to l})$ (Table 4) surprisingly beats the performances of the model trained only on the high-resolution images $(T^h)$ as displayed in Table 3. This is probably a side effect of the increase of the size of the training dataset (Table 2). Nevertheless, when validating on the native low-resolution images $(V^l)$ (Table 4), the performances are relatively poor as compared to the model trained only on the synthetic low-resolution images $(T_{gm}^{l \to h})$. This is explained by the lower quality of the native low-resolution images as already described in the previous section.

### 3.4. Upsampling with the Super Resolution Method Improves the Performances of Plant Detection on the Native Low-Resolution Images.
If the training is difficult with the native low-resolution images because plants are visually difficult to identify and label, the training should be done over low-resolution images derived from the high-resolution images using a more realistic upsampling method than the standard bicubic interpolation one. Alternatively, the training could be done using the high-resolution images and the low-resolution dataset may be upsampled to a synthetic high-resolution domain using bicubic interpolation or super resolution techniques.

Results show that the super resolution technique improved plant detection very significantly as compared to the native low-resolution $(V^l)$ and bicubic upsampled $(V_{bc}^{l \to h})$ images (Table 5). This impacts positively the counting performances while not reaching the performances obtained with the high-resolution images $(V^h)$. The super resolution reduces drastically the underdetection of maize plants particularly on the Tartas site (Figure 10), where as mentioned in Section 3.2, these native low-resolution images have lower textural information and green fraction per plant.

The super resolution approach enhances the features used to identify maize plants, with colors and edges more pronounced than in the corresponding native LR images (Figure 11). Maize plants are visually easier to recognize in the superresolved images as compared to both the native low-resolution and the bicubically upsampled images.

Nevertheless, although easier to interpret, the images generated by super resolution do not appear natural with some exaggerated textural features of the soil background (Figures 11(c) and 11(d)). In few cases, super resolution images show new features—e.g., coloring some pixels in green—in leaf-shaped shadows or tractor tracks in the background (Figure 12) leading to an increase in the proportion of false positives in certain microplots of the Tartas site (Figure 10(b)). Training the super resolution model with a larger dataset might help the generator network to limit those artifacts. Alternatively, some studies [39, 71, 72] have proposed to integrate the training of the super resolution model

<figure><img src="image_13.png" /><figcaption>**FIGURE 12:** An example where the super resolution approach adds undesired artifacts in the image leading to false positives during detection. The SR model adds a few green pixels to the robot tracks on the soil which have "leaf-like" texture. This is wrongly detected as leaf by the Faster-RCNN model. The ground truth boxes are shown in green, and the predicted boxes are shown in blue. The green text indicates the IoU of the predicted box with the ground truth, and the blue text indicates the confidence score of the predictions.</figcaption></figure>

with the training of the Faster-RCNN. The use of a combined detection loss would provide additional information on the location of the plants, thus forcing the super resolution network to differentiate between plants and background while upsampling the images.

In terms of computation speed, the super resolution network takes approximately 20 s to upsample a low-resolution image of 878 × 250 pixels. The detection of the maize plants using the Faster-RCNN model takes approximately 1.5 s for a low-resolution image of 878 × 250 pixels with *x* objects whereas it takes roughly 4 s for its high-resolution counterpart of size 623 × 2337 pixels. Thus, the super resolution and prediction on a high-resolution image is almost 10 times computationally more expensive than predicting directly on a low-resolution image. All the computations were clocked using a graphical processing unit NVIDIA GEFORCE 1080ti with a memory of 12 GB.

## 4. Conclusion

We evaluated the performances of automatic maize plant detection from UAV images using deep learning methods. Our results show that the Faster-RCNN model achieved very good plant detection and counting (rRMSE = 0.08) performances when high-resolution images (GSD ≈ 0.3 cm) are used both for training and validation. However, when this model is applied to the low-resolution images acquired at higher altitudes, the detection and counting performances degrade drastically with rRMSE = 0.48. We demonstrated that this was mostly due to the hyperspecialization of Faster-RCNN that is expecting plants of similar size as in the training dataset. The sensitivity of the detection method to the object size is a critical issue for plant phenotyping applications, where datasets can be generated from different platforms (UAVs, ground vehicles, portable imaging systems, etc.) each one of them providing images within at a specific ground resolution. Concurrently, it would be optimal to share labeled images to get a wide training dataset. Data augmentation techniques where high- and low-resolution images populate the training
dataset were proved to be efficient and provide performances similar to the ones achieved when the model is trained and validated over the same image resolution. However, the native low-resolution images acquired from the UAV have significant low quality that prevents accurate plant detection. In some cases, the images are difficult to visually interpret which poses a problem both for their labeling and for the detector to localize plants due to the lack of pertinent information. These low-quality images were characterized by a loss of image texture that could come from camera intrinsic performances, inadequate settings, and the jpg recording format. It is thus recommended to pay great attention to the camera choice, settings, and recording format when the UAV is flying at altitudes that provide resolution coarser than 0.3 cm for maize plant counting.

In the future studies, it would be worth evaluating the model performances over datasets acquired over a range of flying altitudes, to identify an optimal flying altitude and data quality for plant detection. In our study, it was demonstrated that the quality of the synthetic low-resolution dataset was highly dependent on the low-altitude images used for resampling. Thus, evaluating the model performances at a range of GSDs using additional synthetic low-resolution datasets from Gaussian motion blur resampling would not be representative of the real-world acquisition conditions. It would hence be more pertinent to develop specific metrics allowing to evaluate the richness of textural information contained in the images, rather than using GSD as the main criteria to evaluate image quality.

Finally, we evaluated a super resolution Cycle-ESRGAN-based method to partially overcome the problem of suboptimal image quality. The super resolution method significantly improved the results on the native low-resolution dataset compared to the classic bicubic upsampling strategies. However, the performances when applied to the native low-resolution images were moderate and far poorer than those obtained with the native high-resolution images with simulated superresolved images showing sometimes artifacts. A future direction to reduce the artifacts of such super resolution algorithms can be to integrate the GAN training along with the training of the plant detection network. Another direction would be to introduce some labeled low-resolution images in the training dataset to possibly integrate their features in model. It would also be worth evaluating the performance of more recent object detection networks for plant counting tasks. For instance, one-stage object detection networks such as Yolo-v5 [73]/Yolo-v4 [74] or RetinaNet [75] that outperform Faster-RCNN, would increase the data processing throughput. However, their ability to handle small-scale objects and sensitivity to data quality needs to be studied.

## Data Availability

The UAV data (microplot extractions with bounding boxes) used to support the findings of this study are available from the corresponding authors upon request.

## Conflicts of Interest

The authors declare that there is no conflict of interest regarding the publication of this article.

## Authors' Contributions

KV, RL, SM, and FB designed the study. KV implemented the super resolution pipeline and conducted the analysis. KV and RL contributed extensively to the writing of the article. RL and FB supervised the study. AC and WG participated in discussions and in writing the manuscript. All authors read, revised, and approved the final manuscript.

## Acknowledgments

We would like to thank the CAPTE team for contributing to the construction of the labeled dataset used in this study. We received support from ANRT for the CIFRE grant of KV, cofunded by Hiphen.

## References

[1] R. F. Holt and D. R. Timmons, "Influence of precipitation, soil water, and plant population interactions on corn grain yields¹," *Agronomy Journal*, vol. 60, no. 4, pp. 379-381, 1968.

[2] W. Xu, C. Liu, K. Wang et al., "Adjusting maize plant density to different climatic conditions across a large longitudinal distance in China," *Field Crops Research*, vol. 212, pp. 126-134, 2017.

[3] Y. Zhao, S. Xing, Q. Zhang, F. Zhang, and W. Ma, "Causes of maize density loss in farmers' fields in Northeast China," *Journal of Integrative Agriculture*, vol. 18, no. 8, pp. 1680-1689, 2019.

[4] Y. Gan, E. H. Stobbe, and J. Moes, "Relative date of wheat seedling emergence and its impact on grain yield," *Crop Science*, vol. 32, no. 5, pp. 1275-1281, 1992.

[5] J. S. Graybill, W. J. Cox, and D. J. Otis, "Yield and quality of forage maize as influenced by hybrid, planting date, and plant density," *Agronomy Journal*, vol. 83, no. 3, pp. 559-564, 1991.

[6] R. A. Fischer, O. H. Moreno Ramos, I. Ortiz Monasterio, and K. D. Sayre, "Yield response to plant density, row spacing and raised beds in low latitude spring wheat with ample soil resources: an update," *Field Crops Research*, vol. 232, pp. 95-105, 2019.

[7] G. Maddonni, M. Otegui, and A. Cirilo, "Plant population density, row spacing and hybrid effects on maize canopy architecture and light attenuation," *Field Crops Research*, vol. 71, no. 3, pp. 183-193, 2001.

[8] L. Li, J. Sun, F. Zhang, X. Li, S. Yang, and Z. Rengel, "Wheat/maize or wheat/soybean strip intercropping: I. Yield advantage and interspecific interactions on nutrients," *Field Crops Research*, vol. 71, no. 2, pp. 123-137, 2001.

[9] J. Olsen, L. Kristensen, and J. Weiner, "Influence of sowing density and spatial pattern of spring wheat (Triticum aestivum) on the suppression of different weed species," *Weed Biology and Management*, vol. 6, no. 3, pp. 165-173, 2006.

[10] D. B. Egli and M. Rucker, "Seed vigor and the uniformity of emergence of corn seedlings," *Crop Science*, vol. 52, no. 6, pp. 2774-2782, 2012.
[11] N. W. Hopper, J. R. Overholt, and J. R. Martin, "Effect of culti-
var, temperature and seed size on the germination and emer-
gence of soya beans (Glycine max (L.) Merr.)," *Annals of
Botany*, vol. 44, no. 3, pp. 301–308, 1979.

[12] P. V. V. Prasad, K. J. Boote, J. M. G. Thomas, L. H. Allen, and D. W. Gorbet, "Influence of soil temperature on seedling emergence and early growth of peanut cultivars in field conditions," *Journal of Agronomy and Crop Science*, vol. 192, no. 3, pp. 168–177, 2006.

[13] Z. Berzsenyi and I. S. Tokatlidis, "Density dependence rather than maturity determines hybrid selection in dryland maize production," *Agronomy Journal*, vol. 104, no. 2, pp. 331–336, 2012.

[14] D. Deery, J. Jimenez-Berni, H. Jones, X. Sirault, and R. Furbank, "Proximal remote sensing buggies and potential applications for field-based phenotyping," *Agronomy*, vol. 4, no. 3, pp. 349–379, 2014.

[15] S. Madec, F. Baret, B. de Solan et al., "High-throughput phenotyping of plant height: comparing unmanned aerial vehicles and ground LiDAR estimates," *Frontiers in Plant Science*, vol. 8, p. 2002, 2017.

[16] A. Ruckelshausen, P. Biber, M. Dorna et al., "BoniRob-an autonomous field robot platform for individual plant phenotyping," *Precision Agriculture*, vol. 9, no. 841, p. 1, 2009.

[17] Y. Shi, J. A. Thomasson, S. C. Murray et al., "Unmanned aerial vehicles for high-throughput phenotyping and agronomic research," *PLoS One*, vol. 11, no. 7, article e0159781, 2016.

[18] G. Yang, J. Liu, C. Zhao et al., "Unmanned aerial vehicle remote sensing for field-based crop phenotyping: current status and perspectives," *Frontiers in Plant Science*, vol. 8, p. 1111, 2017.

[19] X. Jin, S. Liu, F. Baret, M. Hemerlé, and A. Comar, "Estimates of plant density of wheat crops at emergence from very low altitude UAV imagery," *Remote Sensing of Environment*, vol. 198, pp. 105–114, 2017.

[20] E. David, G. Daubige, F. Joudelat et al., "Plant detection and counting from high-resolution RGB images acquired from UAVs: comparison between deep-learning and handcrafted methods with application to maize, sugar beet, and sunflower crops," *bioRxiv*, 2021.

[21] D. Stroppiana, M. Pepe, M. Boschetti et al., "Estimating crop density from multi-spectral uav imagery in maize crop," *ISPRS - Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci.*, vol. XLII-2/W13, pp. 619–624, 2019.

[22] P. Randerlović, V. Đorđević, S. Milić et al., "Prediction of soybean plant density using a machine learning model and vegetation indices extracted from RGB images taken with a UAV," *Agronomy*, vol. 10, no. 8, p. 1108, 2020.

[23] B. Li, X. Xu, J. Han et al., "The estimation of crop emergence in potatoes by UAV RGB imagery," *Plant Methods*, vol. 15, no. 1, p. 15, 2019.

[24] D. S. Shrestha and B. L. Steward, "Shape and size analysis of corn plant canopies for plant population and spacing sensing," *Applied Engineering in Agriculture*, vol. 21, no. 2, pp. 295–303, 2005.

[25] S. Liu, F. Baret, B. Andrieu, P. Burger, and M. Hemmerlé, "Estimation of wheat plant density at early stages using high resolution imagery," *Frontiers in Plant Science*, vol. 8, p. 739, 2017.

[26] A. Karami, M. Crawford, and E. J. Delp, "Automatic plant counting and location based on a few-shot learning technique,"

*IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing*, vol. 13, pp. 5872–5886, 2020.

[27] B. T. Kitano, C. C. T. Mendes, A. R. Geus, H. C. Oliveira, and J. R. Souza, "Corn plant counting using deep learning and UAV images," *IEEE Geoscience and Remote Sensing Letters*, pp. 1–5, 2019.

[28] J. Ribera, Y. Chen, C. Boomsma, and E. J. Delp, "Counting plants using deep learning," in *2017 IEEE Global Conference on Signal and Information Processing (GlobalSIP)*, pp. 1344–1348, Montreal, QC, Canada, November 2017.

[29] M. Z. Alom, T. M. Taha, C. Yakopcic et al., "A state-of-the-art survey on deep learning theory and architectures," *Electronics*, vol. 8, no. 3, p. 292, 2019.

[30] E. David, S. Madec, P. Sadeghi-Tehran et al., "Global wheat head detection (GWHD) dataset: a large and diverse dataset of high-resolution RGB-labelled images to develop and benchmark wheat head detection methods," *Plant Phenomics*, vol. 2020, article 3521852, 12 pages, 2020.

[31] A. Dobrescu, M. V. Giuffrida, and S. A. Tsaftaris, "Leveraging multiple datasets for deep leaf counting," in *2017 IEEE International Conference on Computer Vision Workshops (ICCVW)*, pp. 2072–2079, Venice, Italy, October 2017.

[32] S. Ghosal, B. Zheng, S. C. Chapman et al., "A weakly supervised deep learning framework for sorghum head detection and counting," *Plant Phenomics*, vol. 2019, article 1525874, 14 pages, 2019.

[33] D. P. Hughes and M. Salathe, "An open access repository of images on plant health to enable the development of mobile disease diagnostics," 2015, https://arxiv.org/abs/1511.08060.

[34] M. Minervini, A. Fischbach, H. Scharr, and S. A. Tsaftaris, "Finely-grained annotated datasets for image-based plant phenotyping," *Pattern Recognition Letters*, vol. 81, pp. 80–89, 2016.

[35] M. Kisantal, Z. Wojna, J. Murawski, J. Naruniec, and K. Cho, "Augmentation for small object detection," 2019, https://arxiv.org/abs/1902.07296.

[36] K. Tong, Y. Wu, and F. Zhou, "Recent advances in small object detection based on deep learning: A review," *Image and Vision Computing*, vol. 97, p. 103910, 2020.

[37] D. Dai, Y. Wang, Y. Chen, and L. Van Gool, "Is Image Super-Resolution Helpful for Other Vision Tasks?," in *IEEE Winter Conference on Applications of Computer Vision (WACV)*, pp. 1–9, 2016.

[38] M. Fromm, M. Berrendorf, E. Faerman, Y. Chen, and B. Sch, *XD-STOD: cross-domain super resolution for tiny object detection*, pp. 142–148, 2019.

[39] V. Magoulianitis, D. Ataloglou, A. Dimou, D. Zarpalas, and P. Daras, "Does deep super-resolution enhance UAV detection?," in *2019 16th IEEE International Conference on Advanced Video and Signal Based Surveillance (AVSS)*, pp. 1–6, Taipei, Taiwan, September 2019.

[40] C. Dong, C. C. Loy, K. He, and X. Tang, "Image super-resolution using deep convolutional networks," *IEEE Transactions on Pattern Analysis and Machine Intelligence*, vol. 38, no. 2, pp. 295–307, 2016.

[41] A. Lugmayr, M. Danelljan, and R. Timofte, "Unsupervised learning for real-World world super-resolution," in *IEEE/CVF International Conference on Computer Vision Workshop (ICCVW)*, pp. 3408–3416, 2019.

[42] X. Wang, K. Yu, S. Wu et al., "ESRGAN: enhanced super-resolution generative adversarial networks," in *Lecture Notes*
in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics), vol. 11133, pp. 63–79, LNCS, 2019.

[43] K. Zhang, W. Zuo, and L. Zhang, “Learning a single convolutional super-resolution network for multiple degradations,” in *2018 IEEE/CVF Conference on Computer Vision and Pattern Recognition*, pp. 3262–3271, Salt Lake City, UT, USA, June 2018.

[44] M. Fritsche, S. Gu, and R. Timofte, “Frequency separation for real-world super-resolution,” in *2019 IEEE/CVF International Conference on Computer Vision Workshop (ICCVW)*, pp. 3599–3608, Seoul, Korea (South), October 2019.

[45] Tzutalin, “LabelImg,” Git code, 2015, https://github.com/tzutalin/labelImg.

[46] S. Ren, K. He, R. Girshick, and J. Sun, “Faster R-CNN: towards real-time object detection with region proposal networks,” *IEEE Transactions on Pattern Analysis and Machine Intelligence*, vol. 39, no. 6, pp. 1137–1149, 2017.

[47] X. Jin, S. Madec, D. Dutartre, B. de Solan, A. Comar, and F. Baret, “High-throughput measurements of stem characteristics to estimate ear density and above-ground biomass,” *Plant Phenomics*, vol. 2019, article 4820305, 10 pages, 2019.

[48] Y. Liu, C. Cen, Y. Che, R. Ke, Y. Ma, and Y. Ma, “Detection of maize tassels from UAV RGB imagery with faster R-CNN,” *Remote Sens.*, vol. 12, no. 2, p. 338, 2020.

[49] S. Madec, X. Jin, H. Lu et al., “Ear density estimation from high resolution RGB imagery using deep learning technique,” *Agric. For. Meteorol.*, vol. 264, pp. 225–234, 2019.

[50] K. Chen, J. Wang, J. Pang et al., “MMDetection: open MMLab detection toolbox and benchmark,” 2019, arXiv preprint arXiv:1906.07155.

[51] J. Shermeyer and A. Van Etten, “The effects of super-resolution on object detection performance in satellite imagery,” in *Proceedings of the IEEE/CVF Conference on Computer Vision and Pattern Recognition (CVPR) Workshops*, 2019.

[52] A. Clark, “Pillow (PIL Fork) documentation,” 2015, https://buildmedia.readthedocs.org/media/pdf/pillow/latest/pillow.pdf.

[53] H. Zhen, E. Dai, X. Jia et al., *Unsupervised image super-resolution with an indirect supervised path*, 2019, arXiv preprint arXiv:1910.02593.

[54] J.-Y. Zhu, T. Park, P. Isola, and A. A. Efros, “Unpaired image-to-image translation using cycle-consistent adversarial networks,” in *2017 IEEE International Conference on Computer Vision (ICCV)*, pp. 2242–2251, Venice, Italy, October 2017.

[55] E. Agustsson and R. Timofte, “NTIRE 2017 challenge on single image super-resolution: methods and results,” in *2017 IEEE Conference on Computer Vision and Pattern Recognition Workshops (CVPRW)*, pp. 126–135, Honolulu, HI, USA, July 2017.

[56] F. Chollet, “Keras,” GitHub, 2015, https://keras.io.

[57] T. Y. Lin, M. Maire, S. Belongie et al., “Microsoft COCO: common objects in context,” in *Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics)*, vol. 8693, PART 5, pp. 740–755, LNCS, 2014.

[58] R. Padilla, S. L. Netto, and E. A. B. da Silva, “A survey on performance metrics for object-detection algorithms,” in *2020 International Conference on Systems, Signals and Image Processing (IWSSIP)*, pp. 237–242, Niterói, Brazil, July 2020.

[59] T. Kong, A. Yao, Y. Chen, and F. Sun, “HyperNet: towards accurate region proposal generation and joint object detection,” in *2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR)*, pp. 845–853, Las Vegas, NV, USA, June 2016.

[60] H. Jiang and E. Learned-Miller, “Face detection with the faster R-CNN,” in *2017 12th IEEE International Conference on Automatic Face & Gesture Recognition (FG 2017)*, pp. 650–657, Washington, DC, USA, May 2017.

[61] S. Zhang, R. Wu, K. Xu, J. Wang, and W. Sun, “R-CNN-based ship detection from high resolution remote sensing imagery,” *Remote Sensing*, vol. 11, no. 6, p. 631, 2019.

[62] K. He, G. Gkioxari, P. Dollár, and R. Girshick, “Mask R-CNN,” *IEEE Transactions on Pattern Analysis and Machine Intelligence*, vol. 42, no. 2, pp. 386–397, 2020.

[63] M. Machefer, F. Lemarchand, V. Bonnefond, A. Hitchins, and P. Sidiropoulos, “Mask R-CNN refitting strategy for plant counting and sizing in uav imagery,” *Remote Sensing*, vol. 12, no. 18, p. 3015, 2020.

[64] A. Comar, P. Burger, B. de Solan, F. Baret, F. Daumard, and J. F. Hanocq, “A semi-automatic system for high throughput phenotyping wheat cultivars in-field conditions: description and first results,” *Functional Plant Biology*, vol. 39, no. 11, pp. 914–924, 2012.

[65] J. W. White and M. M. Conley, “A flexible, low-cost cart for proximal sensing,” *Crop Science*, vol. 53, no. 4, pp. 1646–1649, 2013.

[66] G. Quaglia, C. Visconte, L. S. Scimmi, M. Melchiorre, P. Cavallone, and S. Pastorelli, “Design of a UGV powered by solar energy for precision agriculture,” *Robotics*, vol. 9, no. 1, p. 13, 2020.

[67] R. Khanna, J. Rehder, M. Moeller, E. Galceran, and R. Siegwart, *Studying phenotypic variability in crops using a hand-held sensor platform*, IROS Work. Agri-Food Robot, 2015.

[68] J. L. Crain, Y. Wei, J. Barker III et al., “Development and deployment of a portable field phenotyping platform,” *Crop Science*, vol. 56, no. 3, pp. 965–975, 2016.

[69] S. Wu, W. Wen, Y. Wang et al., “MVS-Pheno: a portable and low-cost phenotyping platform for maize shoots using multi-view stereo 3D reconstruction,” *Plant Phenomics*, vol. 2020, article 1848437, 17 pages, 2020.

[70] N. Virlet, K. Sabermanesh, P. Sadeghi-Tehran, and M. J. Hawkesford, “Field Scanalyzer: an automated robotic field phenotyping platform for detailed crop monitoring,” *Functional Plant Biology*, vol. 44, no. 1, p. 143, 2017.

[71] M. Haris, G. Shakhnarovich, and N. Ukita, *Task-driven super resolution: object detection in low-resolution images*, 2018, arXiv preprint arXiv:1803.11316.

[72] H. Ji, Z. Gao, T. Mei, and B. Ramesh, “Improved faster R-CNN with multiscale feature fusion and homography augmentation for vehicle detection in remote sensing images,” *IEEE Geoscience and Remote Sensing Letters*, pp. 1–5, 2019.

[73] G. Jocher, A. Stoken, J. Borovec, L. Changyu, and A. Hogan, *ultralytics/yolov5: v3.1 - bug fixes and performance improvements*, Zenodo, 2020.

[74] A. Bochkovskiy, C.-Y. Wang, and H.-Y. M. Liao, “YOLOv4: optimal speed and accuracy of object detection,” 2020, https://arxiv.org/abs/2004.10934.

[75] T.-Y. Lin, P. Goyal, R. Girshick, K. He, and P. Dollár, “Focal loss for dense object detection,” in *Proceedings of the IEEE international conference on computer vision*, pp. 2980–2988, 2017.