Predicting cover crop biomass by lightweight UAS-based RGB and NIR photography: an applied photogrammetric approach

Easy-to-capture and robust plant status indicators are important factors when implementing precision agriculture techniques on fields. In this study, aerial red, green and blue color space (RGB) photography and near-infrared (NIR) photography was performed on an experimental field site with nine different cover crops. A lightweight unmanned aerial system (UAS) served as platform, consumer cameras as sensors. Photos were photogrammetrically processed to orthophotos and digital surface models (DSMs). In a first validation step, the spatial precision of RGB orthophotos (x and y, ± 0.1 m) and DSMs (z, ± 0.1 m) was determined. Then, canopy cover (CC), plant height (PH), normalized differenced vegetation index (NDVI), red edge inflection point (REIP), and green red vegetation index (GRVI) were extracted. In a second validation step, the PHs derived from the DSMs were compared with ground truth ruler measurements. A strong linear relationship was observed (R ² = 0.80?0.84). Finally, destructive biomass samples were taken and compared with the remotely-sensed characteristics. Biomass correlated best with plant height (PH), and good approximations with linear regressions were found (R ² = 0.74 for four selected species, R ² = 0.58 for all nine species). CC and the vegetation indices (VIs) showed less significant and less strong overall correlations, but performed well for certain species. It is therefore evident that the use of DSM-based PHs provides a feasible approach to a species-independent non-destructive biomass determination, where the performance of VIs is more species-dependent.     

Retrieval of biophysical vegetation parameters using simultaneous inversion of high resolution remote sensing imagery constrained by a vegetation index

This study proposes a new method for inverting radiative transfer models to retrieve canopy biophysical parameters using remote sensing imagery. The inversion procedure is improved with respect to standard inversion, and achieves simultaneous inversion of leaf area index (LAI), soil reflectance (ρ_soil), chlorophyll content (C_a+b) and average leaf angle (ALA). In this approach, LAI is used to constrain modelling conditions during the inversion process, providing information about the phenological state of each plot under study. Due to the small area of the vegetation plots used for the inversion procedure and in order to avoid redundant information and improve computation efficiency, existing plot segmentation was used. All retrieved biophysical parameters, except LAI, were assumed to be invariant within each plot. The proposed methodology, based on the combination of PROSPECT and SAILH models, was tested over 16 cereal fields and 51 plots, on two dates, which were chosen to ensure crop assessment at different phenological stages. Plots were selected to provide a wide range of LAI between 0 and 6. Field measurements of LAI, ALA and C_a+b were conducted and used as ground truth for validation of the proposed model-inversion methodology. The approach was applied to very high spatial resolution remote sensing data from the QuickBird 2 satellite. The inversion procedure was successfully applied to the imagery and retrieved LAI with R ² = 0.83 and RMSE = 0.63 when compared to LAI2000 ground measurements. Separate inversions for barley and wheat yielded R ² = 0.89 (RMSE = 0.64) and R ² = 0.56 (RMSE = 0.61), respectively.     

Improving in-season estimation of rice yield potential and responsiveness to topdressing nitrogen application with Crop Circle active crop canopy sensor

In-season site-specific nitrogen (N) management is a promising strategy to improve crop N use efficiency and reduce risks of environmental contamination. To successfully implement such precision management strategies, it is important to accurately estimate yield potential without additional topdressing N application (YP₀) as well as precisely assess the responsiveness to additional N application (RI) during the growing season. Previous research has mainly used normalized difference vegetation index (NDVI) or ratio vegetation index (RVI) obtained from GreenSeeker active crop canopy sensor with two fixed bands in red and near-infrared (NIR) spectrums to estimate these two parameters. The development of three-band Crop Circle active sensor provides a potential to improve in-season estimation of YP₀ and RI. The objectives of this study were twofold: (1) identify important vegetation indices obtained from Crop Circle ACS-470 sensor for estimating rice YP₀ and RI; and (2) evaluate their potential improvements over GreenSeeker NDVI and RVI. Four site-years of field N rate experiments were conducted in 2012 and 2013 at the Jiansanjiang Experiment Station of China Agricultural University located in Northeast China. The GreenSeeker and Crop Circle ACS-470 active canopy sensor with green, red edge, and NIR bands were used to collect rice canopy reflectance data at different key growth stages. The results indicated that both the GreenSeeker (best R² = 0.66 and 0.70, respectively) and Crop Circle (best R² = 0.71 and 0.77, respectively) sensors worked well for estimating YP₀ and RI at the stem elongation stage. At the booting stage, Crop Circle red edge optimized soil adjusted vegetation index (REOSAVI, R² = 0.82) and green ratio vegetation index (R² = 0.73) explained 26 and 22 % more variability in YP₀ and RI, respectively, than GreenSeeker NDVI or RVI. At the heading stage, the GreenSeeker sensor indices became saturated and consequently could not be used for YP₀ or RI estimation, while Crop Circle REOSAVI and normalized green index could still explain more than 70 % of YP₀ and RI variability. It is concluded that both sensors performed similarly at the stem elongation stage, but significantly better results were obtained by the Crop Circle sensor at the booting and heading stages. Furthermore, the results revealed that Crop Circle green band-based vegetation indices performed well for RI estimation while the red edge-based vegetation indices were the best for estimating YP₀ at later growth stages.     

A double-layer model for improving the estimation of wheat canopy nitrogen content from unmanned aerial vehicle multispectral imagery

The accurate and rapid estimation of canopy nitrogen content (CNC) in crops is the key to optimizing in-season nitrogen fertilizer application in precision agriculture. However, the determination of CNC from field sampling data for leaf area index (LAI), canopy photosynthetic pigments (CPP; including chlorophyll a, chlorophyll b and carotenoids) and leaf nitrogen concentration (LNC) can be time-consuming and costly. Here we evaluated the use of high-precision unmanned aerial vehicle (UAV) multispectral imagery for estimating the LAI, CPP and CNC of winter wheat over the whole growth period. A total of 23 spectral features (SFs; five original spectrum bands, 17 vegetation indices and the gray scale of the RGB image) and eight texture features (TFs; contrast, entropy, variance, mean, homogeneity, dissimilarity, second moment, and correlation) were selected as inputs for the models. Six machine learning methods, i.e., multiple stepwise regression (MSR), support vector regression (SVR), gradient boosting decision tree (GBDT), Gaussian process regression (GPR), back propagation neural network (BPNN) and radial basis function neural network (RBFNN), were compared for the retrieval of winter wheat LAI, CPP and CNC values, and a double-layer model was proposed for estimating CNC based on LAI and CPP. The results showed that the inversion of winter wheat LAI, CPP and CNC by the combination of SFs+TFs greatly improved the estimation accuracy compared with that by using only the SFs. The RBFNN and BPNN models outperformed the other machine learning models in estimating winter wheat LAI, CPP and CNC. The proposed double-layer models (R²=0.67–0.89, RMSE=13.63–23.71 mg g^–1, MAE=10.75–17.59 mg g^–1) performed better than the direct inversion models (R²=0.61–0.80, RMSE=18.01–25.12 mg g^–1, MAE=12.96–18.88 mg g^–1) in estimating winter wheat CNC. The best winter wheat CNC accuracy was obtained by the double-layer RBFNN model with SFs+TFs as inputs (R²=0.89, RMSE=13.63 mg g^–1, MAE=10.75 mg g^–1). The results of this study can provide guidance for the accurate and rapid determination of winter wheat canopy nitrogen content in the field.     

基于无人机多光谱遥感的春玉米叶面积指数和地上部生物量估算模型比较研究

明确基于无人机多光谱遥感的玉米叶面积指数（LAI）和地上部生物量的最优估算模型对获取即时、无损、可靠的长势关键参量具有重要意义。2018—2019年,以郑单958（ZD958）和先玉335（XY335）为研究对象,设置4个施氮处理,通过无人机搭载多光谱相机获取多光谱影像,分析两品种LAI和地上部生物量与植被指数相关性,分别构建了基于植被指数的LAI和地上部生物量预测模型。结果表明：同一植被指数在两品种中对施氮量的变化响应规律不同;在吐丝期,幂函数对ZD958的LAI和地上部生物量估算效果最好,指数函数对XY335的LAI估算效果好,幂函数对地上部生物量估算效果好;在灌浆期,幂函数对两品种的LAI估算效果最佳,而指数函数对两品种的地上部生物量估算效果最好。研究结果为进一步提高春玉米长势监测的精度提供了重要依据。     

新疆喀纳斯国家自然保护区植被叶面积指数观测与遥感估算

叶面积指数(Leaf Area Index,LAI)是重要的植被结构参数,调控着植被与大气之间的物质与能量交换,在生态环境脆弱的我国西北部开展植被LAI的研究对阐明该地区植被对气候变化和人类活动的响应特征具有重要的科学意义。利用LAI-2200和TRAC仪器观测了新疆喀纳斯国家级自然保护区森林和草地的有效叶面积指数(LAI_e)和真实LAI,构建了其遥感估算模型,生成了研究区LAI_e和LAI的空间分布图。在此基础上,分析了LAI随地形因子(海拔、坡度、坡向)的变化特征,探讨了将其应用于估算研究区森林生物量密度的可行性,并评估了研究区MODIS LAI产品的精度。结果表明:研究区阔叶林、针阔混交林、针叶林、草地LAI_e的平均值分别为4.40、3.18、2.57、1.76,LAI的平均值分别为4.76、3.93、3.27、2.30。LAI_e和LAI的高值主要集中分布在湖泊和河流附近;植被LAI随海拔、坡度和坡向的变化表现出明显的垂直地带性的特点。LAI随海拔和坡度的增加呈现先增加后减小的变化趋势,坡向对针叶林和草地LAI的影响明显,但对阔叶林和针阔混交林LAI的影响较弱;森林生物量密度(BD)随LAI增加而线性增加(BD=44.396LAI-25.946,R²=0.83),研究区森林生物量密度平均值为120.3 t/hm²,估算的总生物量为5.0×10⁶t;MODIS LAI产品与利用TM数据生成的LAI之间具有一定的相似性(森林R²=0.42,草地R²=0.53),但森林和草地的MODIS LAI产品分别比利用TM数据生成的LAI偏低16.5%和24.4%。     

Leaf area index estimation in vineyards using a ground-based LiDAR scanner

Estimation of grapevine vigour using mobile proximal sensors can provide an indirect method for determining grape yield and quality. Of the various indexes related to the characteristics of grapevine foliage, the leaf area index (LAI) is probably the most widely used in viticulture. To assess the feasibility of using light detection and ranging (LiDAR) sensors for predicting the LAI, several field trials were performed using a tractor-mounted LiDAR system. This system measured the crop in a transverse direction along the rows of vines and geometric and structural parameters were computed. The parameters evaluated were the height of the vines (H), the cross-sectional area (A), the canopy volume (V) and the tree area index (TAI). This last parameter was formulated as the ratio of the crop estimated area per unit ground area, using a local Poisson distribution to approximate the laser beam transmission probability within vines. In order to compare the calculated indexes with the actual values of LAI, the scanned vines were defoliated to obtain LAI values for different row sections. Linear regression analysis showed a good correlation (R ² = 0.81) between canopy volume and the measured values of LAI for 1 m long sections. Nevertheless, the best estimation of the LAI was given by the TAI (R ² = 0.92) for the same length, confirming LiDAR sensors as an interesting option for foliage characterization of grapevines. However, current limitations exist related to the complexity of data process and to the need to accumulate a sufficient number of scans to adequately estimate the LAI.     

Estimation of bioenergy crop yield and N status by hyperspectral canopy reflectance and partial least square regression

The objective of this study was to compare performance of partial least square regression (PLSR) and best narrowband normalize nitrogen vegetation index (NNVI) linear regression models for predicting N concentration and best narrowband normalize different vegetation index (NDVI) for end of season biomass yield in bioenergy crop production systems. Canopy hyperspectral data was collected using an ASD FieldSpec FR spectroradiometer (350–2500 nm) at monthly intervals in 2012 and 2013. The cropping systems evaluated in the study were perennial grass {mixed grass [50 % switchgrass (Panicum virgatum L.), 25 % Indian grass “Cheyenne” (Sorghastrum nutans (L.) Nash) and 25 % big bluestem “Kaw” (Andropogon gerardii Vitman)] and switchgrass “Alamo”} and high biomass sorghum “Blade 5200” (Sorghum bicolor (L.) Moench) grown under variable N applications rates to estimate biomass yield and quality. The NNVI was computed with the wavebands pair of 400 and 510 nm for the high biomass sorghum and 1500 and 2260 nm for the perennial grass that were strongly correlated to N concentration for both years. Wavebands used in computing best narrowband NDVI were highly variable, but the wavebands from the red edge region (710–740 nm) provided the best correlation. Narrowband NDVI was weakly correlated with final biomass yield of perennial grass (r² = 0.30 and RMSE = 1.6 Mg ha^?1 in 2012 and r² = 0.37 and RMSE = 4.0 Mg ha^?1, but was strongly correlated for the high biomass sorghum in 2013 (r² = 0.72 and RMSE = 4.6 Mg ha^?1). Compared to the best narrowband VI, the RMSE of the PLSR model was 19–41 % lower for estimating N concentration and 4.2–100 % lower for final biomass. These results indicates that PLSR might be best for predicting the final biomass yield using spectral sample obtained in June to July, but narrowband NNVI was more robust and useful in predicting N concentration.     

融合无人机光谱信息与纹理信息的冬小麦生物量估测

【目的】生物量是表征植被生命活动的重要参数,对植被长势监测、产量预测有重要意义。以无人机为平台的高光谱遥感技术,具有机动灵活、成本低、空间覆盖广的优势,能够及时准确地估测植被生物量,已成为遥感估算研究的热点之一。由于光谱特征反演生物量存在饱和问题,因此,本研究尝试结合纹理特征与植被指数构建一种"图-谱"融合指标,探究"图-谱"融合指标的抗饱和能力及生物量估测能力。【方法】首先,利用无人机高光谱影像,提取其光谱信息和纹理信息,分别基于植被指数和纹理特征构建生物量模型。其次,针对光谱特征存在的饱和问题,将植被指数与对生物量敏感的纹理指标相乘或相除两种形式构建"图-谱"融合指标,分析"图-谱"融合指标的饱和性,并基于"图-谱"融合指标构建生物量估算模型。最后,对比不同指标构建的生物量模型的估测效果,来分析"图-谱"融合指标估测生物量的能力。【结果】(1)植被指数多在LAI=5时出现饱和现象,而"图-谱"融合指标VI×sm658,VI/ent658,VI/dis658,VI/con658,VI/dis514,VI/con514,VI/var514,VI×con802,VI×dis802均在LAI5时才出现饱和现象,相比之下,这些"图-谱"融合指标一定程度上改善了饱和问题;(2)与植被指数相比(除了GNDVI、NDVI之外),抗饱和能力提高的"图-谱"融合指标VI×sm658、VI/ent658、VI/dis658、VI/con658、VI/dis514、VI/con514、VI/var514、VI×con802、VI×dis802,其与生物量的相关性也相对提高,所构建的生物量模型精度较高(R2=0.81,RMSE=826.02 kg·hm-2)。(3)对比单一植被指数、纹理特征,将纹理特征与光谱特征相结合的"图-谱"融合指标估算小麦生物量的能力相对最强,模型精度明显高于单一植被指数(R2=0.69)和单一纹理特征(R2=0.71)构建的生物量模型。【结论】"图-谱"融合指标的抗饱和能力明显提高,其构建的生物量模型精度也有效提高,实现了结合光谱信息和纹理信息的冬小麦生物量遥感估测,为生物量定量反演提供一种新思路。     

Methodology for measuring fAPAR in crops using a combination of active optical and linear irradiance sensors: a case study in Triticale (X Triticosecale Wittmack)

The amount of photosynthetically active radiation (PAR, 0.4–0.7 μm) absorbed by plants for photosynthesis relative to incident radiation is defined as the fraction of absorbed photosynthetically active radiation (fAPAR). This is an important variable in both plant biomass production and plant growth modeling. This study investigates the application of a newly developed, linear irradiance sensor (LightScout Quantum Bar Sensor, LightScout, Spectrum Technologies, Inc. USA), to quantify fAPAR for a demonstrator crop, Triticale (X Triticosecale Wittmack). A protocol was devised for sensor placement to determine reflected PAR components of fAPAR and to determine the optimal time of day and sensor orientation for data collection. Coincident, top of canopy, normalized difference vegetation index (NDVI) measurements were also acquired with a CropCircle? ACS-210 sensor and measurements correlated with derived fAPAR values. The optimum height of the linear irradiance sensor above soil or plant canopy was found to be 0.4 m while measuring reflected PAR. Measurement of fAPAR was found to be stable when conducted within 1 h of local solar noon in order to avoid significant bidirectional effects resulting from diurnal changes of leaf orientation relative to the vertically-placed sensor. In the row crop studied, averaging fAPAR readings derived from the linear irradiance sensor orientated across and along the plant row provided an R² = 0.81 correlation with above-canopy NDVI. Across row sensor orientation also gave a similar correlation of R² = 0.76 allowing the user to reduce sampling time.     

Crop height variability detection in a single field by multi-temporal terrestrial laser scanning

Information on crop height, crop growth and biomass distribution is important for crop management and environmental modelling. For the determination of these parameters, terrestrial laser scanning in combination with real-time kinematic GPS (RTK–GPS) measurements was conducted in a multi-temporal approach in two consecutive years within a single field. Therefore, a time-of-flight laser scanner was mounted on a tripod. For georeferencing of the point clouds, all eight to nine positions of the laser scanner and several reflective targets were measured by RTK–GPS. The surveys were carried out three to four times during the growing periods of 2008 (sugar-beet) and 2009 (mainly winter barley). Crop surface models were established for every survey date with a horizontal resolution of 1 m, which can be used to derive maps of plant height and plant growth. The detected crop heights were consistent with observations from panoramic images and manual measurements (R² = 0.53, RMSE = 0.1 m). Topographic and soil parameters were used for statistical analysis of the detected variability of crop height and significant correlations were found. Regression analysis (R² < 0.31) emphasized the uncertainty of basic relations between the selected parameters and crop height variability within one field. Likewise, these patterns compared with the normalized difference vegetation index (NDVI) derived from satellite imagery show only minor significant correlations (r < 0.44).     

Phenotyping a diversity panel of quinoa using UAV-retrieved leaf area index,SPAD-based chlorophyll and a random forest approach

Given its high nutritional value and capacity to grow in harsh environments, quinoa has significant potential to address a range of food security concerns. Monitoring the development of phenotypic traits during field trials can provide insights into the varieties best suited to specific environmental conditions and management strategies. Unmanned aerial vehicles (UAVs) provide a promising means for phenotyping and offer the potential for new insights into relative plant performance. During a field trial exploring 141 quinoa accessions, a UAV-based multispectral camera was deployed to retrieve leaf area index (LAI) and SPAD-based chlorophyll across 378 control and 378 saline-irrigated plots using a random forest regression approach based on both individual spectral bands and 25 different vegetation indices (VIs) derived from the multispectral imagery. Results show that most VIs had stronger correlation with the LAI and SPAD-based chlorophyll measurements than individual bands. VIs including the red-edge band had high importance in SPAD-based chlorophyll predictions, while VIs including the near infrared band (but not the red-edge band) improved LAI prediction models. When applied to individual treatments (i.e. control or saline), the models trained using all data (i.e. both control and saline data) achieved high mapping accuracies for LAI (R²?=?0.977–0.980, RMSE?=?0.119–0.167) and SPAD-based chlorophyll (R²?=?0.983–0.986, RMSE?=?2.535–2.861). Overall, the study demonstrated that UAV-based remote sensing is not only useful for retrieving important phenotypic traits of quinoa, but that machine learning models trained on all available measurements can provide robust predictions for abiotic stress experiments.

     

Leaf area index based nitrogen diagnosis in irrigated lowland rice

Leaf area index(LAI)is used for crop growth monitoring in agronomic research,and is promising to diagnose the nitrogen(N)status of crops.This study was conducted to develop appropriate LAI-based N diagnostic models in irrigated lowland rice.Four field experiments were carried out in Jiangsu Province of East China from 2009 to 2014.Different N application rates and plant densities were used to generate contrasting conditions of N availability or population densities in rice.LAI was determined by LI-3000,and estimated indirectly by LAI-2000 during vegetative growth period.Group and individual plant characters(e.g.,tiller number(TN)and plant height(H))were investigated simultaneously.Two N indicators of plant N accumulation(NA)and N nutrition index(NNI)were measured as well.A calibration equation(LAI=1.7787LAI_(2000)–0.8816,R~2=0.870~(**))was developed for LAI-2000.The linear regression analysis showed a significant relationship between NA and actual LAI(R~2=0.863~(**)).For the NNI,the relative LAI(R~2=0.808~(**))was a relatively unbiased variable in the regression than the LAI(R~2=0.33~(**)).The results were used to formulate two LAI-based N diagnostic models for irrigated lowland rice(NA=29.778LAI–5.9397;NNI=0.7705RLAI+0.2764).Finally,a simple LAI deterministic model was developed to estimate the actual LAI using the characters of TN and H(LAI=–0.3375(TH×H×0.01)~2+3.665(TH×H×0.01)–1.8249,R~2=0.875~(**)).With these models,the N status of rice can be diagnosed conveniently in the field.     

A light-weight multi-spectral aerial imaging system for nitrogen crop monitoring

Image-based remote sensing is one promising technique for precision crop management. In this study, the use of an ultra light aircraft (ULA) equipped with broadband imaging sensors based on commercial digital cameras was investigated to characterize crop nitrogen status in cases of combined nitrogen and water stress. The acquisition system was composed of two Canon^? EOS 400D digital cameras: an original RGB camera measuring luminance in the Red, Green and Blue spectral bands, and a modified camera equipped with an external band-pass filter measuring luminance in the near-infrared. A 5?month experiment was conducted on a sugarcane (Saccharum officinarum) trial consisting of three replicates. In each replicate, two sugarcane cultivars were grown with two levels of water input (rainfed/irrigated) and three levels of nitrogen (0, 65 and 130?kg/ha). Six ULA flights, coupled with ground crop measurements, took place during the experiment. For nitrogen status characterisation, three indices were tested from the closed canopy: the normalised difference vegetation index (NDVI), the green normalised difference vegetation index (GNDVI), and a broadband version of the simple ratio pigment index (hereafter referred to as the SRPI_b), calculated from the ratio between blue and red bands of the digital camera. The indices were compared with two nitrogen crop variables: leaf nitrogen content (N_L) and canopy nitrogen content (N_C). SRPI_b showed the best correlation (R ²?=?0.7) with N_L, independently of the water and the N treatment. NDVI and GNDVI were best correlated with N_C values with correlation coefficients of 0.7 and 0.64 respectively, but the regression coefficients were dependent on the water and N treatment. These results showed that SRPI_b could characterise the nitrogen status of sugarcane crop, even in the case of combined stress, and that such acquisition systems are promising for crop nitrogen monitoring.     

Prediction of leaf area index in almonds by vegetation indexes

Three levels of scale for determining leaf area index (LAI) were explored within an almond orchard of alternating rows of Nonpareil and Monterey varieties using hemispherical photography and mule lightbar (MLB) at ground level up to airborne and satellite imagery. We compared LAI estimates of 56 fisheye photos strategically placed in the orchard to validate 500,000 MLB point scans of a small portion of the aisles between tree rows to water and vegetation indexes of MASTER (MODIS/ASTER simulator) and Landsat 5 imagery. The high correlation of fisheye photo LAI to MLB LAI estimates establishes this new method against the measurement standard within the plant community while significantly increasing sample size. MLB LAI and MASTER vegetation indexes, such as NDWI (normalized difference water index), GMI (Gitelson–Merzlyak index) and NDVI (normalized difference vegetation index), were highly correlated (r² = 0.90). In addition, a high correlation (r² = 0.80) between the MLB measured LAI and selected Landsat derived vegetation indexes (VI) was found. This scaling and validation of LAI estimate expands the spatial area and frequency of determination for time series analysis of crop phenology studies.     

Spectral indices from aerial images and their relationship with properties of a corn crop

Identification of areas with similar restrictions to crop productivity could improve the efficiency to manage agricultural systems, guarantee stable yields, and reduce the effect of droughts in rainfed systems. The ability of any vegetation index to discriminate N and moisture-related changes in leaf reflectance would present an important advantage over the present diagnostic system which involves soil-testing for moisture and available N. The purpose of the study was to calibrate different vegetation indices regarding their capacity to identify water and nitrogen availability for rainfed corn crops in the semiarid Pampas of Argentina. A field experiment with corn with a control without fertilization (N0), and fertilized with 120 kg ha^?1 of nitrogen (N120) was used. Two sites, Low (L) and High (H), were identified within the field, according to their altimetry, a multi-spectral aerial photography was taken from a manned airplane during flowering stage of the corn crop, and four spectral indices were calculated (NDVI, green NDVI, NGRDI, (NIR/GREEN)-1). At six georeferenced points at each site soil texture, organic matter, available phosphorus, nitrogen and moisture contents as well as corn aerial biomass and grain yield were determined. The two sites differed in most of the evaluated soil properties, crop biomass and grain yield. The spectral information obtained at crop flowering showed clear differences between sites H and L for all four indices, indicating that any of these would be able to detect the differences in soil moisture and fertility among these environments. Both (NIR/GREEN)-1 and green NDVI had the best correlation with crop yield determined in the field, and therefore could be considered most appropriate for estimating corn yields from images taken at flowering. For estimation of N requirements, green NDVI differentiated best between fertilized and non-fertilized crop in the moisture limited environment (H), while (NIR/GREEN)-1 performed better in the site where soil moisture was non-limiting (L).     

Evaluation of spectral indices and continuous wavelet analysis to quantify aphid infestation in wheat

Wheat aphid, Sitobion avenae F. is one of the most destructive insects infesting winter wheat and appears almost annually in northwest China. Past studies have demonstrated the potential of remote sensing for detecting crop diseases and insect damage. This study aimed to investigate the spectroscopic estimation of leaf aphid density by applying continuous wavelet analysis to the reflectance spectra (350–2 500 nm) of 60 winter wheat leaf samples. Continuous wavelet transform (CWT) was performed on each of the reflectance spectra to generate a wavelet power scalogram compiled as a function of wavelength location and scale of decomposition. Linear regression between the wavelet power and aphid density was to identify wavelet features (coefficients) that might be the most sensitive to aphid density. The results identified five wavelet features between 350 and 2 500 nm that provided strong correlations with leaf aphid density. Spectral indices commonly used to monitor crop stresses were also employed to estimate aphid density. Multivariate linear regression models based on six sensitivity spectral indices or five wavelet features were established to estimate aphid density. The results showed that the model with five wavelet features (R² = 0.72, RMSE = 16.87) performed better than the model with six sensitivity spectral indices (R² = 0.56, RMSE = 21.19), suggesting that the spectral features extracted through CWT might potentially reflect aphid density. The results also provided a new method for estimating aphid density using remote sensing.     

Development of an in-season estimate of yield potential utilizing optical crop sensors and soil moisture data for winter wheat

When utilizing optical sensors to make in-season agronomic recommendations in winter wheat, one parameter often required is the in-season grain yield potential at the time of sensing. Current estimates use an estimate of biomass, such as normalized difference vegetation index (NDVI), and growing degree days (GDDs) from planting to NDVI data collection. The objective of this study was to incorporate soil moisture data to improve the ability to predict final grain yield in-season. Crop NDVI, GDDs that were adjusted based upon if there was adequate water for crop growth, and the amount of soil profile (0–0.80 m) water were incorporated into a multiple linear regression model to predict final grain yield. Twenty-two site-years of N fertility trials with in-season grain yield predictions for growth stages ranging from Feekes 3 to 10 were utilized to calibrate the model. Three models were developed: one for all soil types, one for loamy soil textured sites, and one for coarse soil textured sites. The models were validated with 11 independent site-years of NDVI and weather data. The results indicated there was no added benefit to having separate models based upon soil types. Typically, the models that included soil moisture, more accurately predicted final grain yield. Across all site years and growth stages, yield prediction estimates that included soil moisture had an R² = 0.49, while the current model without a soil moisture adjustment had an R² = 0.40.     

Shadow effect on multi-spectral image for detection of nitrogen deficiency in corn

The efficiency of side-dressing, a more efficient of nitrogen application method than uniform application in either late Fall or early Spring, relies heavily on the capability of nitrogen deficiency detection on a sprayer. To determine the site-specific yield potential for corn, multi-spectral image analysis including aerial- and ground-based images has been used. Some acceptable calibration relationships between the multi-spectral reflectance and SPAD readings have been found from previous study. In sunny weather conditions there was a shadow in the image made by corn leaf itself. This research investigated the shadow effect on the image for detecting corn nitrogen deficiency based on corn canopy reflectance information. The results indicated that the reflectance of red channel in shadow area showed strong inverse correlation, so the vegetation index NDVI using red and NIR channels also showed strong correlation (R² = 77) compared to the whole leaf and bright area. And the reflectance (green and red) and vegetation index(G_NDVI, NDVI, and ratio) in shadow area showed more consistent correlations than others using these image analysis methods.     

Developing a process-based and remote sensing driven crop yield model for maize(PRYM–Maize) and its validation over the Northeast China Plain

Spatial dynamics of crop yield provide useful information for improving the production. High sensitivity of crop growth models to uncertainties in input factors and parameters and relatively coarse parameterizations in conventional remote sensing(RS) approaches limited their applications over broad regions. In this study, a process-based and remote sensing driven crop yield model for maize(PRYM–Maize) was developed to estimate regional maize yield, and it was implemented using eight data-model coupling strategies(DMCSs) over the Northeast China Plain(NECP). Simulations under eight DMCSs were validated against the prefecture-level statistics(2010–2012) reported by National Bureau of Statistics of China, and inter-compared. The 3-year averaged result could give more robust estimate than the yearly simulation for maize yield over space. A 3-year averaged validation showed that prefecture-level estimates by PRYM–Maize under DMCS8, which coupled with the development stage(DVS)-based grain-filling algorithm and RS phenology information and leaf area index(LAI), had higher correlation(R, 0.61) and smaller root mean standard error(RMSE, 1.33 t ha~(–1)) with the statistics than did PRYM–Maize under other DMCSs. The result also demonstrated that DVS-based grain-filling algorithm worked better for maize yield than did the harvest index(HI)-based method, and both RS phenology information and LAI worked for improving regional maize yield estimate. These results demonstrate that the developed PRYM–Maize under DMCS8 gives reasonable estimates for maize yield and provides scientific basis facilitating the understanding the spatial variations of maize yield over the NECP.