- Open Access
Unmanned aerial platform-based multi-spectral imaging for field phenotyping of maize
© Zaman-Allah et al. 2015
Received: 30 April 2015
Accepted: 9 June 2015
Published: 24 June 2015
Recent developments in unmanned aerial platforms (UAP) have provided research opportunities in assessing land allocation and crop physiological traits, including response to abiotic and biotic stresses. UAP-based remote sensing can be used to rapidly and cost-effectively phenotype large numbers of plots and field trials in a dynamic way using time series. This is anticipated to have tremendous implications for progress in crop genetic improvement.
We present the use of a UAP equipped with sensors for multispectral imaging in spatial field variability assessment and phenotyping for low-nitrogen (low-N) stress tolerance in maize. Multispectral aerial images were used to (1) characterize experimental fields for spatial soil-nitrogen variability and (2) derive indices for crop performance under low-N stress. Overall, results showed that the aerial platform enables to effectively characterize spatial field variation and assess crop performance under low-N stress. The Normalized Difference Vegetation Index (NDVI) data derived from spectral imaging presented a strong correlation with ground-measured NDVI, crop senescence index and grain yield.
This work suggests that the aerial sensing platform designed for phenotyping studies has the potential to effectively assist in crop genetic improvement against abiotic stresses like low-N provided that sensors have enough resolution for plot level data collection. Limitations and future potential uses are also discussed.
To ensure improved agricultural productivity, the development and deployment of phenotyping technologies that enable monitoring of phenotypic changes of crop plants in the field is a critical component [1, 2]. Satellite imaging technologies have become an extremely useful tool for collecting data useful for various agricultural applications. However, the major challenges that limit their application in the area of crop improvement are the high cost and the lack of resolution for plot level crop data collection as well as the large revisit periods. The use of manned airborne remote sensing has demonstrated capabilities for large scale crop condition monitoring or for example yield and quality forecasting due to the high spatial and spectral resolutions of the sensors mounted. However in the case of breeding and except for big seed companies its high operating costs and the operational complexity involved have usually limited its use so far to research activities . Unmanned aerial vehicle platforms (UAPs) equipped with sensors are emerging as an important, albeit affordable, component of precision agriculture and crop improvement . The use of these platforms is becoming critical in crop phenotyping because of their ability to rapidly phenotype large numbers of plots and field trials in a dynamic way that can assist the identification and definition of the genetics behind crop yield variability. In addition, the traditional methods currently used, like visual senescence and plant vigor scorings deliver ranking that are variable, depending on the training and subjective appreciations of the staff devoted to that task. With optimum spatial and spectral resolutions, remote sensing from satellite and conventional aerial platforms can provide spatially and spectrally derived parameters for various purposes including crop condition [5–7], crop forecasting and yield predictions [8–10], disease detection and nutrient deficiency [11–13], and photosynthetic pigment content [14–16]. This becomes extremely important with regard to the increasing demand to support and accelerate progress in breeding for novel traits which at the same time requires to accurately measure increasingly large numbers of plants. With improvements in spatial, spectral and temporal resolution of aerial remote sensing, UAPs will enable near real-time visual assessment for crop monitoring in the field yield predictions, crop status mapping, weed detection, and disease and nutrient deficiency detection. Moreover the development of these miniaturized, affordable light-weight unmanned aerial platforms, with better flight control, have enabled the acquisition of high resolution images for various remote sensing applications. Preliminary reports on identification of damaged leaves using the normalized difference vegetation index (NDVI) showed a good similarity between the NDVI values as predicted by remote sensing using UAPs with that of ground truth . Similarly, good potential from multispectral imaging sensors mounted on UAV platforms for physiological condition assessment [18, 19] and stress detection in different crops was reported  including hyperspectral imaging. In addition, these studies reported a good relationship between the predicted and validated values of leaf area index (LAI) in maize (r2 = 0.5) plants as well as chlorophyll concentration at the crop level.
Spectral measurements enable to derive a number of reflectance vegetation indices which have been introduced in both field research and breeding programs for large-scale phenotyping and dynamic estimation of biomass, greenness, nitrogen content, pigment composition, photosynthetic status, and water content [20, 21]. However so far their use in plant phenotyping under field conditions remains far more novel than their implementation under controlled (e.g. greenhouse or growth chamber) conditions .
Another critical area where aerial remote sensing can be useful is the characterization of spatial field variability which results usually from crop management history, spatial changes in soil characteristics and elevation gradients affecting water and nutrients movement. Spatial variability is a serious limitation to breeding efficiency because it creates variation of the stress level imposed within trials, which decreases the heritability of the phenotypic traits evaluated  and obstruct the detection of the genetic signal . Spatial variability in crop productivity is even more evident when differences in resources such as soil N become limiting . In low-input management systems Verhulst et al.  showed that standard deviation and coefficient of variation of NDVI values were high. Recently, Cairns et al.  have reported a very large residual and genotype × trial variances in the drought and combined drought and heat treatments relative to the well-watered treatment which resulted in reduced heritability for means estimates from the stressed trials. In addition, a combined analysis of the southern Africa regional trials of CIMMYT and partners also pointed out plot residual variance to be much higher under managed stress relative to nonstress trials . These results highlight the need for measures to reduce the effects of field variability so as to increase the genetic signal to noise ratio. One way of addressing this spatial field variability problem is to collect soil information; but this has proven to be a laborious process. Therefore, aerial spectral imaging could be a quick and low-cost method for experimental field characterization.
To be relevant for breeding, plant phenotyping, should allow to objectively select key trait(s) under the least spatial field variability conditions. This underlines the critical need to use the “right” tools for data collection and for minimizing the spatial variability. To date, only few studies have reported attempts to use UAP’s remote sensing for spatial field variation assessment and crop phenotyping in the field. This work reports a proof of concept exercise on how a UAV-based remote sensing platform equipped with sensor for multispectral imaging could be used for experimental field spatial variability and maize phenotyping in the field under low-nitrogen stress conditions.
Results and discussion
Experimental field characterization
Descriptive statistics for the NDVI values at CIMMYT’s Harare low-N fields, Harare, Zimbabwe
Plant stress detection
Nitrogen status of field crops can be assessed using leaf or canopy spectral reflectance data [6, 30]. Several studies have found that nondestructive measurements of leaf or canopy reflectance can be used for detecting N-deficient stress in maize [31, 32], rice , and wheat . Using aerial multispectral imaging with a UAP, we evaluated the capabilities for remotely sensing low-N stress in maize hybrids. NDVI data generated from multispectral imaging were used to compute a low-N stress index.
Our data have shown that remote sensing using a UAP can be used to detect subtle differences of N stress within a field to a resolution of a single row and yet evaluate an entire field (Figure 2). This level of resolution together with the fast data collection that the platform allow, open an avenue for in-season descriptive analysis of plant development variables. With more exploration, this will permit determination of the critical growth stage to consider, the appropriate spectral bands for crop performance analysis, and the integration of co-variate in statistical designs or even the integration of data in crop models.
Comparison between ground-based measurements and UAV-based remote sensing
Leaf senescence affects the plant’s ability to fill the grains by reducing the grain filling duration. Crop senescence can be used to indirectly assess the ability of a genotype to maintain a higher plant photosynthetic capacity under N deficiency conditions.
In addition the variation of the leaf senescence index in relation to N rates show that at low N range, the slope is ten times higher than at high N rates and the relationship is stronger.
Plant leaf senescence decreased with increase in N application rate (Figure 6a). As a consequence of the delay in senescence plant photosynthetic capacity prolongs the effective leaf area duration . On the opposite side, N deficiency accelerates leaf senescence rates throughout the life cycle . This has been shown to be important for ear and kernel initiation, contributing to define maize sink capacity  and maintaining functional kernels throughout grain filling with positive impact on the number of developed kernels and kernel final size . Under low soil N conditions, the photosynthetic capacity is reduced due to early senescence as compared to high soil N. The crop senescence index was 20% smaller at 160 kg ha−1 AN relative to 0 kg ha−1 AN. It was reported that leaf senescence could explain 47% of genotypic variation in nitrogen use efficiency in the field experiments .
Relationship with grain yield
Descriptive statistics for the Grain yield values (t ha−1) in 10 maize hybrids grown under 6 N application rates
N-application rate (Kg ha−1 AN)
Coefficients of correlation between grain yield (GY) and (1) NDVI extracted from multispectral images taken with the UAV-platform (NDVI-UAP) and (2) leaf senescence index
Crop senescence index
All N application rates
0 and 10 kg ha−1 AN
80 and 160 kg ha−1 AN
Regarding crop senescence index, the correlation with grain yield followed the same trend as NDVI (Table 3). The correlation was stronger under low N conditions (0–10 kg.ha−1) as compared to sufficient N conditions. N deficiency is known to accelerate leaf senescence rate  which was shown to be important for ear and kernel initiation, contributing to define maize sink capacity  and maintaining functional kernels throughout grain filling .
The use of UAP for field spatial variability and field-based crop phenotyping is novel, but is expected to become an important tool for improving efficiency in crop breeding. Currently, most of the limitations in the deployment of these platforms in breeding are related to the cost of sensors, spatial resolution of imagery, data processing, management and complexity of operation. The results of the current study suggest that remote sensing from UAP has a great potential for field and crop trait characterization under field conditions. To date, only few studies have been carried out attempting to use UAP’s remote sensing for spatial field variation assessment and crop phenotyping in the field. Our results showed that, this type of platform can be used in low N stress detection/senescence as well as for estimating final yield in maize. However, to effectively deploy this type of platform in a breeding program, there is need to measure its breeding value through selection indices and to have a well-designed data processing and management plan.
Experimental set up
To address the relevance of UAV-based remote sensing platforms for field crop phenotyping, we set up two experiments at the CIMMYT-Harare research station (−17.725787 S, 31.016457 E) on 5-years nitrogen depleted fields referred here as managed low-N fields. The average Nitrate–N concentration of the managed low-N fields (using the spectrophotometric method) was 8 mg kg−1 soil. The first experiment focused on the characterization of the field site for spatial variability. The field was sown with a single wheat variety (Triticum aestivum L. cv. SC-Stallion). The planting was done during the winter season 2012 under controlled irrigation to ensure that germination and growth were homogenous. The second experiment focused on assessing N-response of maize (Zea mays L.) hybrids. A total of 10 hybrids were used, 5 hybrids were classified as low-N tolerant and 5 hybrids as low-N susceptible based on previous multi-location experiments. Six N fertilization levels were used: 0, 10, 20, 40, 80 and 160 kg ha−1 of Ammonium Nitrate (AN). Except in the case of 80 and 160 kg ha−1 of AN, where a split application was used (50% at knee-high and 50% prior to anthesis), N fertilization was performed at knee-high stage. The experiment was a split-plot design with hybrids as main plots and N levels as subplots, replicated three times in a randomized complete block design. The experiment was planted on the 23rd December 2013 as 2 rows per plot. The rows were 75 cm apart and 4 m long with 17 planting stations per row. To reduce border effects, 2-row plots of a commercial variety were planted on all sides of each treatment. Single super phosphate fertilizer (14% P2O5: 7% K2O) was applied to all the plots at a rate of 400 kg ha−1, to supply the crop with phosphorous and potassium. Hand weeding was used to control weeds.
Unmanned aerial platform
This platform has the capability of flying at a wide range of altitudes and is capable of carrying a payload of up to 1.5 kg and flying for 30 min. It can cover up to 40 ha for a 30 min flight at an altitude of 150 m above ground level.
Imagery acquisition and data collection
Images were collected from the wheat experiment to assess spatial variability 1 week prior to booting and for the whole N-response trial on two dates (anthesis and 2 weeks after anthesis). They were acquired using the ADC-Lite multispectral camera mounted on the UAV platform near midday under cloud-free conditions. The ADC-Lite is a specialized light weight agricultural multispectral camera (ADC-Lite, Tetracam, Inc., Chatsworth, CA, USA) with a 3.2 megapixel Complementary Metal Oxide Semiconductor (CMOS) sensor. The camera’s sensor has a ground resolution of 60 mm per pixel and 123 × 92 m Field of View at an altitude of 150 m above ground level. The camera is optimized to simultaneously capture Green, Red and NIR channels with bands approximately equal to Landsat Thematic Mapper TM2 (520–600 nm), TM3 (630–690 nm) and TM4 (760–900 nm) and has a lens focal length of 8.0 mm. The UAP followed a flight plan made of 6 waypoints covering the entire field at an altitude of 150 m above terrain and a velocity of 45 km h−1. The images (2048 × 1536 pixels) were taken at the rate of one every 5 s.
Different ground-based sensors were used to collect data. A leaf chlorophyll meter (Minolta SPAD-502, Spectrum Technologies Inc, Plainfield, IL, USA) a spectroradiometer provided with an active sensor (GreenSeeker handheld crop sensor, Trimble, USA) as well as a conventional digital camera (Cyber-shot DSC-WX80, Sony, Japan) were used to measure spectroradiometrical leaf canopy (NDVI) and RGB (red/green/blue) image derived vegetation indices, respectively, as detailed in the data procession section.
Image data processing
Color infra-red multispectral images
A radiometric calibration was performed after applying a channel decomposition, which consists of demosaicing the infrared color filter array (CFA) to reconstruct each G-R-NIR sample from the undersampled ones. The imagery was synchronized through the GPS position and triggering time recorded for each image, without the usage of any additional inertial units. For this study, only absolute positions were used to generate the ortho-rectified mosaics, following by the steps of image registration, calibration and blending.
The multispectral images acquired by the UAP enabled identification of each individual row (Figure 2) and the extraction of maize plant reflectance values. For each row, a region of interest (ROI) was established manually in the center of the row based on visual detection of the edges/soil to extract reflectance values (Figure 2a) using only vegetation pixels. The image reflectance data extracted from each treatment field, coinciding with each flight time, were subsequently used to compute the vegetation indices used in the analysis.
The NDVI and yield data were subjected to classical statistical analysis to obtain (1) descriptive statistics mean, minimum, maximum, standard deviation and coefficient of variation and (2) correlation coefficients using GraphPad Prism (GraphPad Software Inc., San Diego, CA, USA, 1996).
ZM and CJ planned, carried out the trials and prepared the manuscript. TA and MC participated in the development and selection of the hybrids used. VO was involved in data collection. AJL, Z-TPJ and HA participated in the design of the paltform, sensors calibration and data analysis process. DB, CP, OM and PMB helped to draft the manuscript. All authors read and approved the final manuscript.
This work was supported by CRP-Maize Grant from CIMMYT, funding (AGL2012-40053-C03-01 and AGL2013-44147-R.) from the Spanish Ministerio de Economia y Competitividad and the Improved Maize for African Soils (IMAS) project.
Compliance with ethical guidelines
Competing interests The authors declare that they have no competing interests.
Open AccessThis 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.
- White JW, Andrade-Sanchez P, Gore MA, Bronson KF, Coffelt TA, Conley MM et al (2012) Field-based phenomics for plant genetics research. Field Crops Res 133:101–112View ArticleGoogle Scholar
- Walter A, Liebisch F, Hund A (2015) Plant phenotyping: from bean weighing to image analysis. Plant Methods 11:14PubMed CentralPubMedView ArticleGoogle Scholar
- Berni J, Zarco-Tejada PJ, Suarez L, Fereres E (2009) Thermal and narrowband multispectral remote sensing for vegetation monitoring from an unmanned aerial vehicle. IEEE Trans Geosci Remote Sens 47(3):722–738View ArticleGoogle Scholar
- Frank L, Norbert K, David S, Achim W, Andreas H (2015) Remote, aerial phenotyping of maize traits with a mobile multi-sensor approach Andreas. Plant Methods 11:9View ArticleGoogle Scholar
- Boissard P, Pointel JG, Huet P (1993) Reflectance, green leaf area index and ear hydric status of wheat from anthesis until maturity. Int J Remote Sens 14:2713–2729View ArticleGoogle Scholar
- Potdar MB (1993) Sorghum yield modeling based on crop growth-parameters determined from visible and near-IR channel NOAA AVHRR data. Int J Remote Sens 14:895–905View ArticleGoogle Scholar
- Blackmer TM, Schepers JS, Varvel GE (1994) Light reflectance compared with other nitrogen stress measurements in corn leaves. Agron J 86(6):934–938View ArticleGoogle Scholar
- Tucker CJ, Holben BN, Elgin JH, Mcmurtrey J (1980) Relationship of spectral data to grain-yield variation. Photogramm Eng Remote Sens 46:657–666Google Scholar
- Clevers JGPW (1997) A simplified approach for yield prediction of sugar beet based on optical remote sensing data. Remote Sens Environ 61(2):221–228View ArticleGoogle Scholar
- Moran M, Maas S, Pinter P Jr (1995) Combining remote sensing and modeling for estimating surface evaporation and biomass production. Remote Sens Rev 12:335–353View ArticleGoogle Scholar
- Malthus TJ, Madeira AC (1993) High resolution spectroradiometry: spectral reflectance of field bean leaves infected by Botrytis Fabae. Remote Sens Environ 45:107–116View ArticleGoogle Scholar
- Adams ML, Philpot WD, Norvell WA (1999) Yellowness index: an application of spectral second derivatives to estimate chlorosis of leaves in stressed vegetation. Int J Remote Sens 20(18):3663–3675View ArticleGoogle Scholar
- Adams ML, Norvell WA, Philpot WD, Peverly JH (2000) Toward the discrimination of manganese, zinc, copper, and iron deficiency in “Bragg” soybean using spectral detection methods. Agron J 92:268–274View ArticleGoogle Scholar
- Sampson PH, Zarco-Tejada PJ, Mohammed GH, Miller JR, Noland TL (2003) Hyperspectral remote sensing of forest condition: estimating chlorophyll content in tolerant hardwoods. For Sci 49:381–391Google Scholar
- Zarco-Tejada P, Berjón A, López-Lozano R, Miller J, Martín P, Cachorro V et al (2005) Assessing vineyard condition with hyperspectral González indices: leaf and canopy reflectance simulation in a row-structured discontinuous canopy. Remote Sens Environ 99:271–287View ArticleGoogle Scholar
- Martín P, Zarco-Tejada P, González M, Berjón A (2007) Using hyperspectral remote sensing to map grape quality in “Tempranillo” vineyards affected by iron deficiency chlorosis. Vitis 46:7–14Google Scholar
- Nebikar S, Annen A, Schurrer M, Oesch D (2008) A light-weight multispectral sensor for micro UAV-opportunities for very high resolution airborne remote sensing. The International Archives of the Photogrammetry, Remote sensing and Spatial Information Science, Vol. XXXVII, Part B, Beijing, ChinaGoogle Scholar
- Zarco-Tejada PJ, González-Dugo V, Berni JAJ (2012) Fluorescence, temperature and narrow-band indices acquired from a UAV platform for water stress detection using a micro-hyperspectral imager and a thermal camera. Remote Sens Environ 117:322–337View ArticleGoogle Scholar
- Zarco-Tejada PJ, Catalina A, González MR, Martín P (2013) Relationships between net photosynthesis and steady-state chlorophyll fluorescence retrieved from airborne hyperspectral imagery. Remote Sens Environ 136:247–258View ArticleGoogle Scholar
- Cairns JE, Sanchez C, Vargas M, Ordoñez RA, Araus JL (2012) Dissecting maize productivity: ideotypes associated with grain yield under drought stress and well-watered conditions. J Integr Plant Biol 54:1007–1020PubMedView ArticleGoogle Scholar
- Fiorani F, Rascher U, Jahnke S, Schurr U (2012) Imaging plants dynamics in heterogenic environments. Curr Opin Biotechnol 23:227–235PubMedView ArticleGoogle Scholar
- Araus JL, Cairns J (2014) Field high-throughput phenotyping—the new crop breeding frontier. Trends Plant Sci 19:52–61PubMedView ArticleGoogle Scholar
- Masuka B, Araus JL, Sonder K, Das B, Cairns JE (2012) Phenotyping for abiotic stress tolerance in maize. J Integr Plant Biol 54:238–249PubMedView ArticleGoogle Scholar
- Campos H, Heard JE, Ibañez M, Luethy MH, Peters TJ, Warner DC (2011) Effective and efficient platforms for crop phenotype characterisation under drought. In: Monneveux P, Ribaut JM (eds) Drought phenotyping in crops: from theory to practice. CGIAR Generation Challenge Programme, Texcoco, pp 39–47Google Scholar
- Bänziger M, Edmeades GO, Beck D, Bellon M (2000) Breeding for drought and nitrogen stress tolerance in Maize: from theory to practice 2000. CIMMYT, MexicoGoogle Scholar
- Verhulst N, Govaerts B, Sayre KD, Deckers J, Francois IM, Dendooven L (2008) Using NDVI and soil quality analysis to assess influence of agronomic management on within plot spatial variability and factors limiting production. Plant Soil J. 317:41–59View ArticleGoogle Scholar
- Cairns JE, Crossa J, Zaidi P, Grudloyma P, Sanchez C, Araus JL et al (2013) Identification of drought, heat, and combined drought and heat tolerant donors in maize (Zea mays L.). Crop Sci 53(4):1335–1346View ArticleGoogle Scholar
- Weber VS, Melchinger AE, Magorokosho C, Makumbi D, Bänziger M, Atlin GN (2012) Efficiency of managed stress screening of elite maize hybrids under drought and low nitrogen for yield under rainfed conditions in southern Africa. Crop Sci 2012(52):1011–1020View ArticleGoogle Scholar
- Lopez-Lozano R, Casterad MA, Herrero J (2010) Site-specific management units in a commercial maize plot delineated using very high resolution remote sensing and soil properties mapping. Comput Electron Agric 73:219–229View ArticleGoogle Scholar
- Osborne SL, Schepers JS, Francis DD, Schlemmer MR (2002) Use of spectral radiance to estimate in-season biomass and grain yield in nitrogen and water-stressed corn. Crop Sci 42:165–171PubMedView ArticleGoogle Scholar
- Blackmer TM, Shepers JS (1996) Aerial photography to detect nitrogen stress in corn. J Crop Phys 148:440–444Google Scholar
- Graeff S, Steffens D, Schubert S (2001) Use of reflectance measurements for the early detection of N, P, Mg, and Fe deficiencies in Zea mays L. J Plant Nutr Soil Sci 164:445–450View ArticleGoogle Scholar
- Wang K, Shen ZQ, Wang RC (1998) Effects of N nutrition on the spectral reflectance characteristics of rice leaf and canopy. J Zhejiang Agric Univ 24:93–97Google Scholar
- Voullot MO, Huet P, Boissard P (1998) Early detection of N deficiency in a wheat crop using physiological and radiometric methods. Agronomie 18:117–130View ArticleGoogle Scholar
- Graeff S, Claupein W (2003) Quantifying nitrogen status of corn (Zea mays L.) in the field by reflectance measurements. Europ J Agron 19:611–618View ArticleGoogle Scholar
- Wiesler F, Behrens T, Horst WJ (2001) Nitrogen efficiency of contrasting rape ideotypes. In: Horst WJ, Schenk MK, Bürkert A, Claassen N, Flessa H, Frommer WB, Goldbach H, Olfs H-W, Römheld V, Sattelmacher B, Schmidhalter U, Schubert S, Wirén Nv, Wittenmayer L (eds) Plant nutrition. Developments in plant and soil sciences, vol 92. Springer, Netherlands, pp 60–61Google Scholar
- Earl HJ, Tollenaar M (1997) Maize leaf absorptance of photosynthetically active radiation and its estimation using a chlorophyll meter. Crop Sci 37(2):436–440View ArticleGoogle Scholar
- Wolfe DH, Henderson DW, Hsiao TC, Alvino A (1988) Integrative water and nitrogen effects on the senescence of maize. I. Leaf area duration, nitrogen distribution and yield. Agron J 80(4):859–864View ArticleGoogle Scholar
- Tollenaar M, McCullough DE, Dwyer LM (1994) Physiological basis of the genetic improvement of corn. In: Slafer GA (ed) Genetic improvement of field crops. M. Dekker, New York, pp 183–236Google Scholar
- Jones RJ, Schreibe BM, Roessler J (1996) Kernel sink capacity in maize: genotypic and maternal regulation. Crop Sci 36:301–306View ArticleGoogle Scholar
- Schulte auf’m Erley G, Wijaya KA, Ulas A, Becker H, Wiesler F, Horst WJ (2007) Leaf senescence and N uptake parameters as selection traits for nitrogen efficiency of oilseed rape cultivars. Physiol Plant 130:519–531View ArticleGoogle Scholar
- Rawson HM, Hindmarsh JH, Fisher RA, Stockman YM (1983) Changes in leaf photosynthesis with plant ontogeny and relationships with yield per ear in wheat cultivars and 120 progeny. Aust J Plant Physiol 10:503–514View ArticleGoogle Scholar
- Ellen J (1987) Effects of plant density and nitrogen fertilization in winter wheat (Triticum aestivum L.): I. Production pattern and grain yield. Neth J Agric Sci 35:137–153Google Scholar
- Tucker CJ (1979) Red and photographic infrared linear combinations for monitoring vegetation. Remote Sensing Environ 8:127–150View ArticleGoogle Scholar
- Stone ML, Solie JB, Whitney RW, Raun WR, Lees HL (1996) Sensors for the detection of nitrogen in winter wheat; SAE Paper No. 961757; 1996 SAE, Warrendale PAGoogle Scholar
- Aparicio N, Villegas D, Casadesus J, Araus JL, Royo C (2000) Spectral vegetation indices as nondestructive tools for determining durum wheat yield. Agron J 92:83–91View ArticleGoogle Scholar
- Casadesús J, Villegas D (2014) Conventional digital cameras as a tool for assessing leaf area index and biomass for cereal breeding. J Integr Plant Biol 56(1):7–14PubMedView ArticleGoogle Scholar
- Casadesús J, Kaya Y, Bort J, Nachit MM, Araus JL, Amor S et al (2007) Using vegetation indices derived from conventional digital cameras as selection criteria for wheat breeding in water-limited environments. Ann Appl Biol 150:227–236View ArticleGoogle Scholar