# A method to estimate plant density and plant spacing heterogeneity: application to wheat crops

- Shouyang Liu
^{1}Email author, - Fred Baret
^{1}, - Denis Allard
^{2}, - Xiuliang Jin
^{1}, - Bruno Andrieu
^{3}, - Philippe Burger
^{4}, - Matthieu Hemmerlé
^{5}and - Alexis Comar
^{5}

**13**:38

**DOI: **10.1186/s13007-017-0187-1

© The Author(s) 2017

**Received: **25 September 2016

**Accepted: **2 May 2017

**Published: **17 May 2017

## Abstract

### Background

Plant density and its non-uniformity drive the competition among plants as well as with weeds. They need thus to be estimated with small uncertainties accuracy. An optimal sampling method is proposed to estimate the plant density in wheat crops from plant counting and reach a given precision.

### Results

Three experiments were conducted in 2014 resulting in 14 plots across varied sowing density, cultivars and environmental conditions. The coordinates of the plants along the row were measured over RGB high resolution images taken from the ground level. Results show that the spacing between consecutive plants along the row direction are independent and follow a gamma distribution under the varied conditions experienced. A gamma count model was then derived to define the optimal sample size required to estimate plant density for a given precision. Results suggest that measuring the length of segments containing 90 plants will achieve a precision better than 10%, independently from the plant density. This approach appears more efficient than the usual method based on fixed length segments where the number of plants are counted: the optimal length for a given precision on the density estimation will depend on the actual plant density. The gamma count model parameters may also be used to quantify the heterogeneity of plant spacing along the row by exploiting the variability between replicated samples. Results show that to achieve a 10% precision on the estimates of the 2 parameters of the gamma model, 200 elementary samples corresponding to the spacing between 2 consecutive plants should be measured.

### Conclusions

This method provides an optimal sampling strategy to estimate the plant density and quantify the plant spacing heterogeneity along the row.

### Keywords

Wheat Gamma-count model Density RGB imagery Sampling strategy Plant spacing heterogeneity## Background

Plant density at emergence is governed by the sowing density and the emergence rate. For a given plant density, the uniformity of plant distribution at emergence may significantly impact the competition among plants as well as with weeds [1, 2]. Plant density and uniformity is therefore a key factor explaining production, although a number of species are able to compensate for low plant densities by a comparatively significant development of individual plants during the growth cycle. For wheat crops which are largely cultivated over the globe, tillering is one of the main mechanisms used by the plant to adapt its development to the available resources that are partly controlled by the number of tillers per unit area. The tillering coefficient therefore appears as an important trait to be measured. It is usually computed as the ratio of the number of tillers per unit area divided by the plant density [3]. Plant density is therefore one of the first variables measured commonly in most agronomical trials.

Crops are generally sown in rows approximately evenly spaced by seedling devices. Precision seedling systems mostly used for crops with plants spaced on the row by more than few centimeters (e.g. maize, sunflower or soybean) distribute seeds relatively evenly along the row. Conversely, for most crops with short distances among plants on the row, e.g. wheat, barley or canola, seeds are distributed non-evenly along the row. This can be attributed both to the mechanisms that free, at a variable frequency, the seed from the seed tank, and the trajectory of the seed that may also vary in the pipe that drives it from the seed tank to the soil. Further, once reaching the soil, the seed may also move with the soil displaced by the sowing elements penetrating the soil surface. Finally, some seeds may abort or some young plants may die because of pests or too extreme local environmental conditions (excess or deficit of moisture, low temperature etc.). The population density and its non-uniformity are therefore recognized as key traits of interests to characterize the canopy at the emergence stage. However, very little work documents the plant distribution pattern along the row, which is partly explained by the lack of dedicated device for accurate plant position measurement [4]. Electromagnetic digitizers are very low throughput and not well adapted to such field measurements [5]. Alternatively, algorithms have been developed to measure the inter-plant spacing along the row for maize crops from top-view RGB (Red Green Blue) images [6, 7]. Improvements were then proposed by using three dimensional sensors [8–10]. However, these algorithms were only validated on maize crops that show relatively simple plant architecture with generally fixed inter-plant spacing along the row.

Manual field counting in wheat crops is still extensively employed as the reference method. Measurements of plant population density should be completed when the majority of plants have just emerged and before the beginning of tillering when individual plants start to be difficult to be identified. Plants are counted over elementary samples corresponding either to a quadrat or to a segment [11]. The elementary samples need to be replicated in the plot to provide a more representative value [12]. For wheat crops, [3] suggested that at least a total of 3 m of rows (0.5 m segment length repeated 6 times) should be counted, while [13] proposed to sample a total of 6 m (segments made of 2 consecutive rows by one meter repeated 3 times in the plot). [14] proposed to repeat at least 4 times the counting in 0.25 m^{2} quadrats corresponding roughly to a total of 6.7 m length of rows (assuming the rows are spaced by 0.15 m). In this case, quadrats may be considered as a set of consecutive row segments with the same length when the quadrat is oriented parallel to the row direction or with variable lengths when the quadrat is oriented differently. Although these recommendations are simple and easy to apply, they may not correspond to an optimal sampling designed to target a given precision level. They may either provide low precision if under sampled or correspond to a waste of human resources in the opposite case.

The sample size required to reach a given precision of the plant density will depend on the population density and the heterogeneity of plant positions along the row that may be described by the distribution of the distances between consecutive plants. This distribution is more likely to be skewed, which could be described by an exponential distribution or a more general one such as the Weibull or the gamma distributions. Fitting such random distribution functions provides not only access to the plant density at the canopy level, but also to its local variation that may impact the development of neighboring plants as discussed earlier.

The objective of this study is to propose an optimal sampling method for plant density estimation and to quantify the heterogeneity of plant spacing along the row. For this purpose, a model is first developed to describe the distribution of the plants along the row. The model is then calibrated over a number of ground experiments. Further, the model is used to compare several plant counting strategies and to evaluate the optimal sampling size to reach a given precision. Finally, the model was also exploited to design a method for quantifying the non-uniformity of plant distribution.

## Methods

### Field experiment

The experimental design in 2014 over the three sites

Sites | Latitude | Longitude | Cultivar | Density (seeds m |
---|---|---|---|---|

Toulouse | 43.5°N | 1.5°E | Apache | 100, 200, 300, 400, 600 |

Grignon | 48.8°N | 1.9°E | Premio; Attlass; Flamenko; Midas; Koréli | 150 |

Avignon | 43.9°N | 4.8°E | Apache | 100, 200, 300, 400 |

### Image processing

The coordinates \( {\text{x}}_{\text{n}} \) of plant *n* (noted \( {\text{Plant}}_{\text{n}} \)) along the row axis allow to compute the spacing \( \Delta {\text{x}}_{\text{n}} = ({\text{x}}_{\text{n}} - {\text{x}}_{{{\text{n}} - 1}} ) \) between \( {\text{Plant}}_{\text{n}} \) and \( {\text{Plant}}_{{{\text{n}} - 1}} . \) The actual plant density expressed in plants per square meter horizontal ground (plants m^{−2}) was computed simply as the number of plants counted on the segments, divided by the product of the length of the segments and the row spacing.

### Development and calibration of the plant distribution model

#### Distribution of plant spacing

*p*value of Chi square test (Table 2) [17]. Besides, the tail of the Weibull distribution tends toward zero less rapidly than that of the gamma distribution: Weibull may show few samples with very large values [18], increasing the risk of overestimation for the larger plant spacing. The gamma distribution was therefore used in the following and writes [19]:

Parameters of the fitted distributions

Sites | Sowing density (seeds m | Cultivar | Gamma | Weibull | ||||
---|---|---|---|---|---|---|---|---|

a | b | p value of Chi square test | a | b | p value of Chi square test | |||

Avignon | 100 | Apache | 1.14 | 6.38 | 0.27 | 7.44 | 1.07 | 0.29 |

200 | Apache | 1.25 | 4.04 | 0.62 | 5.29 | 1.13 | 0.05 | |

300 | Apache | 0.99 | 2.53 | 0.38 | 2.51 | 1.00 | 0.56 | |

400 | Apache | 0.96 | 1.50 | 0.22 | 1.39 | 0.94 | 0.57 | |

Toulouse | 100 | Apache | 1.07 | 5.01 | 0.12 | 5.32 | 0.99 | 0.10 |

200 | Apache | 1.39 | 1.95 | 0.17 | 2.86 | 1.15 | 0.12 | |

300 | Apache | 1.21 | 2.28 | 0.94 | 2.89 | 1.12 | 0.94 | |

400 | Apache | 1.24 | 1.37 | 0.51 | 1.76 | 1.10 | 0.40 | |

600 | Apache | 1.16 | 0.96 | 0.37 | 1.14 | 1.09 | 0.21 | |

Grignon | 150 | Premio | 1.12 | 3.37 | 0.70 | 3.85 | 1.06 | 0.68 |

150 | Attlass | 1.13 | 2.48 | 0.69 | 2.87 | 1.05 | 0.67 | |

150 | Flamenko | 1.11 | 3.3 | 0.92 | 3.75 | 1.05 | 0.92 | |

150 | Midas | 1.24 | 3.03 | 0.21 | 3.92 | 1.12 | 0.24 | |

150 | Koréli | 1.15 | 2.89 | 0.24 | 3.48 | 1.15 | 0.18 |

#### Modeling the distribution of the number of plants per row segment

_{1}is the number of plants over the segment of length l, and \( {\text{IG}}\left( {{\text{a}} \cdot {\text{n}},\frac{l}{b}} \right) \) is the incomplete gamma function:

_{1}, estimated over a segment of length l can be expressed by introducing the row spacing distance, r, assumed to be known:

The expectation, \( {\text{E}}\left( {{\text{D}}_{\text{l}} } \right) \), converges toward the actual density of the population when 1 → ∞.

*a*and

*b*may lead to the same plant density, but with variations in their distribution along the row (Fig. 5). The fitting of parameters

*a*and

*b*over the 14 plots using the transformed gamma-count model (Eq. 9) shows that the shape parameter, a, varies from 0.96 to 1.39 and is quite stable. Conversely, the scale parameter, b, appears to vary widely from 0.96 to 6.38, mainly controlling the plant density (Fig. 5). Since the CV depends only on the shape parameter a (Eq. 4), it should not vary much across the 14 plots considered. This was confirmed by applying a one-way analysis of variance on the CV values of the 14 plots available (F = 1.09, P = 0.3685): no significant differences are observed. This result may be partly explained by the fact that the same type of seed drill was used for all the three sites.

## Results

### Optimal sample size to reach a given precision for plant density estimation

The transformed gamma-count model provides a convenient way to investigate the effect of the sampling size on the precision of the density estimates. The precision will be quantified here using the coefficient of variation (CV). The sample size can be expressed either as a given length of the segments where the (variable) number of plants should be counted, or as a (variable) length of the segment to be measured corresponding to a given number of consecutive plants. The two alternative sampling approaches will be termed FLS (Fixed Length of Segments) for the first one, and FNP (Fixed Number of Plants) for the second one.

*L*, where plants need to be counted. The optimal

*L*value for a given target precision quantified by the CV will mainly depend on the current density as demonstrated in Fig. 6a: longer segments are required for the low densities. Conversely, shorter segments are needed for high values of the plant density to reach the same precision. The scale parameter,

*b*, that controls the plant density drives therefore the optimal segment length

*L*(Fig. 6a). Counting plants over

*L*= 5 m (500 cm) provides a precision better than 10% for densities larger than 150 plants·m

^{−2}for the most common conditions characterized by a shape coefficient

*a*> 0.9. These figures agree well with the usual practice for plant counting as reviewed in the introduction [3, 13, 14]. Increasing the precision quantified by the CV will require longer segments

*L*to be sampled (Fig. 7a).

When considering the FNP approach, the sample size is driven by the number, *N*, of consecutive plants that defines to a row segment whose length need to be measured. The simulations of the model (Fig. 6b) show that *N* is mainly independent from the plant density. For the 14 plots considered in this study, segments with 70 < *N* < 110 plants should be measured to reach a precision of CV = 10%. The shape parameter *a* influences dominantly the sample size: more heterogeneous distribution of plants characterized by small values of the shape parameter will require more plants to be counted (Fig. 6b). To increase the precision (lower CV), more plants will also need to be counted (Fig. 7b).

Optimal sampling size for FSL and FNP over different densities (100, 150, 200, 300, 400 and 600 seeds m^{−2}) and precisions (5, 10 and 15%)

Sowing density (seeds m | Parameters | CV = 5% | CV = 10% | CV = 15% | ||||
---|---|---|---|---|---|---|---|---|

a | b | FSL (cm) | FNP (Nb. Plt) | FSL (cm) | FNP (Nb. Plt) | FSL (cm) | FNP (Nb. Plt) | |

100 | 1.11 | 5.70 | 2478 | 363 | 620 | 90 | 250 | 39 |

150 | 1.15 | 3.01 | 1406 | 348 | 351 | 88 | 130 | 37 |

200 | 1.32 | 3.00 | 1398 | 308 | 350 | 78 | 130 | 33 |

300 | 1.10 | 2.41 | 1162 | 363 | 291 | 90 | 110 | 39 |

400 | 1.10 | 1.44 | 774 | 363 | 194 | 90 | 60 | 39 |

600 | 1.16 | 0.96 | 584 | 348 | 146 | 85 | 50 | 37 |

### Sampling strategy to quantify plant spacing variability on the row

The previous sections demonstrated that the scale and shape parameters could be estimated from the observed distribution of the plant spacing. However, the measurement of individual plant spacing is tedious and prone to errors as outlined earlier. The estimation of these parameters from the variability observed between small row segments containing a fixed number of plants will therefore be investigated here. This FNP approach is preferred here to the FSL one because there will be no additional uncertainties introduced by the position of the first and last plants of the segment with the corresponding start and end of the segment. These uncertainties may be significant in case of small segments in the FLS approach.

_{n}corresponding to

*n*plant spacing between (

*n*+ 1) consecutive plants with \( {\text{L}}_{\text{n}} = \sum \nolimits_{{{\text{i}} = 1}}^{\text{n}} \Delta {\text{x}}_{\text{i}} , \) as a gamma distribution with

*n*·

*a*as shape parameter and the same scale parameter

*b*as the one describing the distribution of ∆X.

The parameters *a* and *b* will therefore be estimated by adjusting the gamma model described in Eq. 12 for the given value of *n* + 1 consecutive plants.

*a*and

*b*parameters estimation was further investigated. A numerical experiment based on a Monte-Carlo approach was conducted considering a standard case corresponding to the average of the 14 plots sampled in 2014 with

*a*= 1.10 and

*b*= 2.27. The sampling size is defined by the number of consecutive plants for the FNP approach considered here and by the number of replicates. For each sampling size 300 samples were generated by randomly drawing in the gamma distribution (Eq. 12) and parameters

*a*and

*b*were estimated. The standard deviation between the 300 estimates of

*a*and

*b*parameters was finally used to compute the corresponding CV. This process was applied to a number of replicates varying between 20 to 300 by steps of 10 and a number of plants per segment varying between 2 (i.e. spacing between two consecutive plants) to 250 within 12 steps. This allows describing the variation of the coefficient of estimated values of parameters

*a*and

*b*as a function of the number of replicates and the number of plants (Fig. 8).

Results show that the sensitivity of the CV of estimates of parameters *a* and *b* are very similar (Fig. 8). The sensitivity of parameters *a* and *b* is dominated by the number of replicates: very little variation of CV is observed when the number of plants per segment varies (Fig. 8). Parameters *a* and *b* require about 200 replicates independently from the number of plants per segment. It seems therefore more interesting to make very small segments to decrease the total number of plants to count.

Additional investigations not shown here for the sake of brevity, confirmed the independency of the number of replicates to the number of plants per segment when parameters *a* and *b* are varying. Further, the number of replicates need to be increased as expected when the shape parameter *a* decreases (i.e. when the plant spacing is more variable) to keep the same precision on estimates of *a* and *b* parameters.

## Discussion and conclusions

A method was proposed to estimate plant density and sowing pattern from high resolution RGB images taken from the ground. The method appears to be much more comfortable as compared with the standard outdoor methods based on plant counting in the field. Images should ideally be taken around Haun stage 1.5 for wheat crops when most plants have already emerged and tillering has not yet started. Great attention should be paid to the geometric correction in order to get accurate ortho-images where distances can be measured accurately. The processing of images here was automatic except the last step corresponding to the interactive visual extraction of the plants’ coordinates in the image. However, recent work [24, 25] suggests that it will be possible to automatize this last step to get a fully high-throughput method.

The method proposed is based on the modeling of the plant distribution along the row. It was first demonstrated that the plant spacing between consecutive plants are independent which corresponds to a very useful simplifying assumption. The distribution of plant spacing was then proved to follow a gamma distribution. Although the Weibull distribution showed similar good performance, it was not selected because of the comparatively heavier tails of the distribution that may create artefacts. Further the Weibull model does not allow to simply derive the distribution law of the length of segments containing several consecutive plants [26]. The gamma model needs a scale parameter that drives mostly the intensity of the process, i.e. the plant density, and a shape parameter that governs the heterogeneity of plant spacing. This model was transformed into a count data model to investigate the optimal sampling required to get an estimate of plant density for a given precision level.

^{−2}over the 14 plots available. The model performs better for the low density with a RMSE of 21 plants m

^{−2}for density lower than 400 plants m

^{−2}. These discrepancies may be mainly explained by the accuracy in the measurement of the position of individual plants (around 1–2 mm). Uncertainties on individual plant spacing will be high in relative values as compared to that associated with the measurement of the length of the segment used in the simple method to get the ‘reference’ plant density. Hence it is obviously even more difficult to get a good accuracy in plant spacing measurements for high density, i.e. with a small distance among plants. In addition, small deviations from the gamma-count model are still possible, although the previous results were showing very good performance.

The model proposed here concerns mainly relatively nominal sowing, i.e. when the sowing was successful on average on the row segments considered: portions of rows with no plants due to sowing problems or local damaging conditions (pests, temperature and moisture). The sowing was considered as nominal on most of the plots investigated in this study, with no obvious ‘accidents’. However, it is possible to automatically identify from the images the unusual row segments with missing plants or excessive concentration of plants [25]. Rather than describing blindly the bulk plant density, it would be then preferred to get a nested sampling strategy: the unusual segments could be mapped extensively, and the plant density of nominal and unusual segments could be described separately using the optimal sampling proposed here.

This study investigated the sampling strategy to estimate the plant density with emphasis on the variability of plant spacing along the row, corresponding to the sampling error. However additional sources of error should be accounted for including measurement biases, uncertainties in row spacing or non-randomness in the sample selection [27– 29]. Unlike sampling error, it could not be minimized by increasing sampling size. The non-sampling error may be reduced by combining a random sampling selection procedure with a measurement method ensuring high accuracy including accounting for the actual values of the row spacing measured over each segment [30].

Optimal sampling requires a tradeoff between minimum sampling error obtained with maximum sampling size and minimum cost obtained with minimum sampling size [31]. The optimal sampling strategy should first be designed according to the precision targeted here quantified by the coefficient of variation (CV) characterizing the relative variability of the estimated plant density between several replicates of the sampling procedure. The term ‘optimal’ should therefore be understood as the minimum sampling effort to be spent to achieve the targeted precision. Two approaches were proposed: the first one considers a fixed segment length (FSL) over which the plants have to be counted; the second one considers a fixed number of successive plants (FNP) defining a row segment, the length of which needs to be measured. The first method (FLS) is the one generally applied within most field experiments. However, we demonstrated that it is generally sub-optimal: since the segment length required to achieve a given CV depends mainly on the actual plant density: the sampling will be either too large for the targeted precision, or conversely too small, leading to possible degradation of the precision of plant density estimates. Nevertheless, for the plant density (>100 plants m^{−2}) and shape parameter (*a* > 0.9) usually experienced, a segment length of 6 m will ensure a precision better than 10%. The second approach (FNP) appears generally more optimal: it aims at measuring the length of the segment corresponding to a number of consecutive plants that will depend mainly on the targeted precision. Results demonstrate that in our conditions, the density should be evaluated over segments containing 90 plants to achieve a 10% precision. The sampling size will always be close to optimal as compared to the first approach where optimality requires the knowledge of the plant density that is to be estimated. Further, the FNP approach is probably more easy to implement with higher reliability: as a matter of facts, measuring the length of a segment defined by plants at its two extremities is easier than counting the number of plants in a fixed length segment, where the extremities could be in the vicinity of a plant and its inclusion or not in the counting could be prone to interpretation biases by the operator. The total number of plants required in a segment could be split into subsamples containing smaller number of plants that will be replicated to get the total number of plants targeted. This will improve the spatial representativeness. Overall, the method proposed meets the requirements defined by [32, 33] for the next genearation of phenotyping tools: increase the accuracy, the precision and the throughput while reducing the labor and budgetary costs.

The gamma-count model proved to be well suited to describe the plant spacing distribution along the row over our contrasted experimental situations. It can thus be used to describe the heterogeneity of plant spacing as suggested by [20]. This may be applied for detailed canopy architecture studies or to quantify the impact of the sowing pattern heterogeneity on inter-plant competition [1, 2]. The heterogeneity of plant spacing may be described by the scale and shape parameters of the gamma model. Quantification of the heterogeneity of plant spacing requires repeated measurements over segments defined by a fixed number of plants. Our results clearly show that the precision on estimates of the gamma count parameters depends only marginally on the number of plants in each segment. Conversely, it depends mainly on the number of segments (replicates) to be measured. For the standard conditions experienced in this study, the optimal sampling strategy to get a CV lower than 10% on the two parameters of the gamma distribution would be to repeat 200 times the measurement of plant spacing between 2 consecutive plants.

## Declarations

### Authors’ contributions

SL and FB designed the experiment and XJ, BA, PB, MH and AC contributed to the field measurement in different experimental sites. DA significantly contributed to the method development. The manuscript was written by SL and significantly improved by FB. All authors read and approved the final manuscript.

### Acknowledgements

We thank the people from Grignon, Toulouse and Avignon who participated to the experiments. Great thanks to Paul Bataillon and Jean-Michel Berceron from UE802 Toulouse INRA for their help in field experiment. The work was completed within the UMT-CAPTE funded by the French Ministry of Agriculture.

### Competing interests

The authors declare that they have no competing interests.

### Availability of data and materials

All data analyzed during this study are presented in this published article.

### Funding

This study was supported by “Programme d’investissement d’Avenir” PHENOME (ANR-11-INBS-012) and Breedwheat (ANR-10-BTR-03) with participation of France Agrimer and “Fonds de Soutien à l’Obtention Végétale”. The grant of the principal author was funded by the Chinese Scholarship Council.

### Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Olsen J, Kristensen L, Weiner J. Influence of sowing density and spatial pattern of spring wheat (
*Triticum aestivum*) on the suppression of different weed species. Weed Biol Manag. 2006;6:165–73.View ArticleGoogle Scholar - Olsen J, Weiner J. The influence of Triticum aestivum density, sowing pattern and nitrogen fertilization on leaf area index and its spatial variation. Basic Appl Ecol. 2007;8:252–7.View ArticleGoogle Scholar
- Reynolds MP, Ortiz-Monasterio JI, McNab A. CIMMYT: application of physiology in wheat breeding. Mexico: CIMMYT; 2001.Google Scholar
- Van der Heijden G, De Visser PHB, Heuvelink E. Measurements for functional-structural plant models. Frontis. 2007;22:13–25.Google Scholar
- Ma Y, Wen M, Guo Y, Li B, Cournede P-H, de Reffye P. Parameter optimization and field validation of the functional-structural model GREENLAB for maize at different population densities. Ann Bot. 2007;101:1185–94.View ArticlePubMedPubMed CentralGoogle Scholar
- Tang L, Tian L. Plant identification in mosaicked crop row images for automatic emerged corn plant spacing measurement. Trans ASABE. 2008;6:2181–91.View ArticleGoogle Scholar
- Tang L, Tian LF. Real-time crop row image reconstruction for automatic emerged corn plant spacing measurement. Trans ASABE. 1079;2008:51.Google Scholar
- Jin J, Tang L. Corn plant sensing using real-time stereo vision. J Field Robot. 2009;26:591–608.View ArticleGoogle Scholar
- Nakarmi AD, Tang L. Automatic inter-plant spacing sensing at early growth stages using a 3D vision sensor. Comput Electron Agric. 2012;82:23–31.View ArticleGoogle Scholar
- Nakarmi AD, Tang L. Within-row spacing sensing of maize plants using 3D computer vision. Biosys Eng. 2014;125:54–64.View ArticleGoogle Scholar
- Norman DW:
*The farming systems approach to development and appropriate technology generation.*Food & Agriculture Org.; 1995. - Marshall MN. Sampling for qualitative research. Fam Pract. 1996;13:522–6.View ArticlePubMedGoogle Scholar
- Gate P. Ecophysiologie du blé. Paris: Tec & Doc Lavoisier; 1995.Google Scholar
- Whaley JM, Sparkes DL, Foulkes MJ, Spink JH, Semere T, Scott RK. The physiological response of winter wheat to reductions in plant density. Ann Appl Biol. 2000;137:165–77.View ArticleGoogle Scholar
- Haun J. Visual quantification of wheat development. Agron J. 1973;65:116–9.View ArticleGoogle Scholar
- Brown M, Lowe DG. Automatic panoramic image stitching using invariant features. Int J Comput Vision. 2006;74:59–73.View ArticleGoogle Scholar
- Plackett RL. Karl pearson and the Chi squared test. Int Stat Rev/Revue Internationale de Statistique. 1983;51:59–72.Google Scholar
- Papalexiou SM, Koutsoyiannis D, Makropoulos C. How extreme is extreme? An assessment of daily rainfall distribution tails. Hydrol Earth Syst Sci. 2013;17:851–62.View ArticleGoogle Scholar
- Winkelmann R. A count data model for gamma waiting times. Stat Pap. 1996;37:177–87.View ArticleGoogle Scholar
- Liu S, Baret F, Andrieu B, Abichou M, Allard D, Solan Bd, Burger P. Modeling the distribution of plants on the row for wheat crops: consequences on the green fraction at the canopy level. Comput Electron Agric. 2016;136:147–56.View ArticleGoogle Scholar
- Winkelmann R. Duration dependence and dispersion in count-data models. J Bus Econ Stat. 1995;13:467–74.Google Scholar
- Zeviani WM, Ribeiro PJ Jr, Bonat WH, Shimakura SE, Muniz JA. The Gamma-count distribution in the analysis of experimental underdispersed data. J Appl Stat. 2014;41:2616–26.View ArticleGoogle Scholar
- Steutel FW, Van Harn K. Infinite divisibility of probability distributions on the real line. Boca Raton: CRC Press; 2003.View ArticleGoogle Scholar
- Jin X, Liu S, Baret F, Hemerlé M, Comar A: Estimates of plant density from images acquired from UAV over wheat crops at emergence. Remote Sens Environ. 2016 (
**accepted**). - Liu S, Baret F, Andrieu B, Burger P, Hemerlé M: Estimation of plant density from high resolution RGB imagery over wheat crops. Front Plant Sci. 2016 (
**accepted**). - Rinne H. The Weibull distribution: a handbook. Boca Raton: CRC Press; 2008.View ArticleGoogle Scholar
- Lohr S. Sampling: design and analysis. Boston: Cengage Learning; 2009.Google Scholar
- Israel GD:
*Determining sample size.*University of Florida Cooperative Extension Service, Institute of Food and Agriculture Sciences, EDIS; 1992. - Särndal C-E, Swensson B, Wretman J. Model assisted survey sampling. New York: Springer Science & Business Media; 2003.Google Scholar
- Scheuren F, Association AS: What is a survey? In. American Statistical Association; 2004.
- Biemer PP. Total survey error: design, implementation, and evaluation. Public Opin Q. 2010;74:817–48.View ArticleGoogle Scholar
- Cobb JN, DeClerck G, Greenberg A, Clark R, McCouch S. Next-generation phenotyping: requirements and strategies for enhancing our understanding of genotype–phenotype relationships and its relevance to crop improvement. Theor Appl Genet. 2013;126:867–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Fiorani F, Schurr U. Future scenarios for plant phenotyping. Annu Rev Plant Biol. 2013;64:267–91.View ArticlePubMedGoogle Scholar