Search

Accurate estimation of crop biophysical and biochemical parameters during crop growing seasons is essential for improving site-specific management and yield estimation. The potential ability of fully polarimetric synthetic aperture radar (SAR) data in estimating above-ground biomass of oilseed rape was investigated in this study. The temporal profile of different scattering intensity and polarimetric features during the entire growing season was identified with ground measurements. A polarimetric feature, relying on the polarimetric decomposition method, was put forward to estimate the biomass of oilseed rape. Validation results revealed great potential with a determination coefficient (R2) of 0.85, root mean squared error (RMSE) of 41.6 g/m2, and relative error (RE) of 28.5% for dry biomass, and an R2 of 0.76, RMSE of 527.4 g/m2 and RE of 28.6% for fresh biomass. Moreover, the use of full polarization SAR data was compared with single and dual polarization SAR data. The results suggest that when full polarization SAR data is available, a simpler model, higher saturation point and better accuracy can be achieved in biomass estimation of oilseed rape, which highlights the importance and value of polarimetry information in quantitative crop monitoring. This study provides guidelines for in-season monitoring of crop growth parameters with SAR data, which further improves crop monitoring capability in adverse weather conditions. [ABSTRACT FROM AUTHOR]

Abstract When dealing with multiple growth stages, estimates of above-ground biomass (AGB) based on optical vegetation indices (VIs) are difficult for two reasons: (i) optical VIs saturate at medium-to-high canopy cover, and (ii) organs that grow vertically (e.g., biomass of reproductive organs and stems) are difficult to detect by canopy spectral VIs. Although several significant improvements have been made for estimating AGB by using narrow-band hyperspectral VIs, synthetic aperture radar, laser intensity direction and ranging, the crop surface model technique, and combinations thereof, applications of these new techniques have been limited by cost, availability, data-processing difficulties, and high dimensionality. The present study thus evaluates the use of ultrahigh-ground-resolution image textures, VIs, and combinations thereof to make multiple temporal estimates and maps of AGB covering three winter-wheat growth stages. The selected gray-tone spatial-dependence matrix-based image textures (e.g., variance, entropy, data range, homogeneity, second moment, dissimilarity, contrast, correlation) are calculated from 1-, 2-, 5-, 10-, 15-, 20-, 25-, and 30-cm-ground-resolution images acquired by using an inexpensive RGB sensor mounted on an unmanned aerial vehicle (UAV). Optical-VI data were obtained by using a ground spectrometer to analyze UAV-acquired RGB images. The accuracy of AGB estimates based on optical VIs varies, with validation R 2: 0.59–0.78, root mean square error (RMSE): 1.22–1.59 t/ha, and mean absolute error (MAE): 1.03–1.27 t/ha. The most accurate AGB estimate was obtained by combining image textures and VIs, which gave R 2: 0.89, MAE: 0.67 t/ha, and RMSE: 0.82 t/ha. The results show that (i) the eight selected textures from ultrahigh-ground-resolution images were significantly related to AGB, (ii) the combined use of image textures from 1- to 30-cm-ground-resolution images and VIs can improve the accuracy of AGB estimates as compared with using only optical VIs or image textures alone; and (iii) high AGB values from winter-wheat reproductive growth stages can be accurately estimated by using this method; (iv) high estimates of winter-wheat AGB (8–14 t/ha) using the proposed combined method (DIS1, SE30, B460, B560, B670, EVI2 using MSR) show a 22.63% (nRMSE) improvement compared with using only spectral VIs (LCI, NDVI using MSR), and a 21.24% (nRMSE) improvement compared with using only image textures (COR1, DIS1, SE30, EN30 using MSR). Thus, the combined use of image textures and VIs can help improve estimates of AGB under conditions of high canopy coverage. [ABSTRACT FROM AUTHOR]

Accurate and high-throughput identification of the initial anthesis of soybean varieties is important for the breeding and screening of high-quality soybean cultivars in field trials. The objectives of this study were to identify the initial day of anthesis (IADAS) of soybean varieties based on remote sensing multispectral time-series images acquired by unmanned aerial vehicles (UAVs), and analyze the differences in the initial anthesis of the same soybean varieties between two different climatic regions, Shijiazhuang (SJZ) and Xuzhou (XZ). First, the temporal dynamics of several key crop growth indicators and spectral indices were analyzed to find an effective indicator that favors the identification of IADAS, including leaf area index (LAI), above-ground biomass (AGB), canopy height (CH), normalized-difference vegetation index (NDVI), red edge chlorophyll index (CIred edge), green normalized-difference vegetation index (GNDVI), enhanced vegetation index (EVI), two-band enhanced vegetation index (EVI2) and normalized-difference red-edge index (NDRE). Next, this study compared several functions, like the symmetric gauss function (SGF), asymmetric gauss function (AGF), double logistic function (DLF), and fourier function (FF), for time-series curve fitting, and then estimated the IADAS of soybean varieties with the first-order derivative maximal feature (FDmax) of the CIred edge phenology curves. The relative thresholds of the CIred edge curves were also used to estimate IADAS, in two ways: a single threshold for all of the soybean varieties, and three different relative thresholds for early, middle, and late anthesis varieties, respectively. Finally, this study presented the variations in the IADAS of the same soybean varieties between two different climatic regions and discussed the probable causal factors. The results showed that CIred edge was more suitable for soybean IADAS identification compared with the other investigated indicators because it had no saturation during the whole crop lifespan. Compared with DLF, AGF and FF, SGF provided a better fitting of the CIred edge time-series curves without overfitting problems, although the coefficient of determination (R2) and root mean square error (RMSE) were not the best. The FDmax of the SGF-fitted CIred edge curve (SGF_CIred edge) provided good estimates of the IADAS, with an RMSE and mean average error (MAE) of 3.79 days and 3.00 days, respectively. The SGF-fitted_CIred edge curve can be used to group the soybean varieties into early, middle and late groups. Additionally, the accuracy of the IADAS was improved (RMSE = 3.69 days and MAE = 3.09 days) by using three different relative thresholds (i.e., RT50, RT55, RT60) for the three flowering groups compared to when using a single threshold (RT50). In addition, it was found that the IADAS of the same soybean varieties varied greatly when planted in two different climatic regions due to the genotype–environment interactions. Overall, this study demonstrated that the IADAS of soybean varieties can be identified efficiently and accurately based on UAV remote sensing multispectral time-series data. [ABSTRACT FROM AUTHOR]

• Vegetation indices show a strong correlation with nitrogen status in tea plants. • The combination of multispectral data and hyperspectral data is effective in monitoring nitrogen level. • Fusing data of image information and spectral information has remarkable ability to predict nitrogen content in tea plants. Nitrogen (N) plays a pivotal role in management of tea plantation, with significant impacts on the growth, productivity, and nutrition status of tea plants. The existing methods for N content monitoring of tea leaves are complicated and can not realize in suite and in real time way. This study proposed a method for estimating the N content of tea plants in field conditions based on a combination of a multispectral imaging system and hyperspectral data. A total of 32 parameters were extracted from five tea gardens using calibrated multispectral images of the tea plant canopy, and 27 indices were selected by Pearson correlation analysis. A total of 28 wavelengths selected by competitive adaptive reweighted sampling from hyperspectral data were combined with 27 multispectral indices as the original data. Subsequently, five variables of fused data (H, VOG, BGI, 1664 nm and 1665 nm) were selected by variable combination population analysis based on the 55 combination parameters. Partial least squares regression, random forest regression, and support vector machine regression (SVR) models all showed excellent performance for both the calibration and prediction sets. The overall results indicated that the infused data of multispectral and hyperspectral data combined with SVR are effective in monitoring the N level under field conditions, and the R2 (coefficient of determination) and root mean square error values of the prediction were 0.9186 and 0.0560, respectively. The findings of this study are important in retaining the nutritional and quality attributes of agricultural commodities. [ABSTRACT FROM AUTHOR]

• Accurately forecast the first flowering period of two typical apple producing areas in China. • It fills in the research of forecasting the first flowering period of apple from point scale to regional scale. • The phenological survey method of replacing time with space is very effective for the calibration of sequential model. • The prediction of the apple first flowering on a regional scale has important application value for optimizing orchard management. China is one of the largest apple-producing countries in the world, with large orchards and diverse climates. Accurately forecasting the first-flowering time of apple trees can assist orchard managers in their deciding when to apply anti-freeze. The temperature-driven sequential model from previous studies can be used to forecast the flowering phenology of deciduous fruit trees. However, this model requires many years of observational data for calibration, so flowering forecasts based on traditional phenological models cannot be implemented in areas that lack such historical data. To overcome this problem, the present work combines a spatial rather than a temporal phenological survey method with 1-km-gridded temperature products to calibrate the chill and heat requirement parameters of the sequential model. We then use the model to forecast the first-flowering on a regional scale for Luochuan and Linyi, which are two main apple-producing areas of China. The results show that the proposed method accurately forecasts regional flowering. The root mean squared errors (RMSE) for Luochuan and Linyi were 4.7 and 4.4 days, respectively, and the normalized RMSEs were all less than 5.19%. We expect the proposed regional first-flowering forecast method to be an important aid to optimize orchard management. [ABSTRACT FROM AUTHOR]

Accurate determination of phenological information of crops is essential for field management and decision-making. Remote sensing time-series data are widely used for extracting phenological phases. Existing methods mainly extract phenological phases directly from individual remote sensing time-series, which are easily affected by clouds, noise, and mixed pixels. This paper proposes a novel method of phenological phase extraction based on the time-weighted dynamic time warping (TWDTW) algorithm using MODIS Normalized Difference Vegetation Index (NDVI) 5-day time-series data with a spatial resolution of 500 m. Firstly, based on the phenological differences between winter wheat and other land cover types, winter wheat distribution is extracted using the TWDTW classification method, and the results show that the overall classification accuracy and Kappa coefficient reach 94.74% and 0.90, respectively. Then, we extract the pure winter-wheat pixels using a method based on the coefficient of variation, and use these pixels to generate the average phenological curve. Next, the difference between each winter-wheat phenological curve and the average winter-wheat phenological curve is quantitatively calculated using the TWDTW algorithm. Finally, the key phenological phases of winter wheat in the study area, namely, the green-up date (GUD), heading date (HD), and maturity date (MD), are determined. The results show that the phenological phase extraction using the TWDTW algorithm has high accuracy. By verification using phenological station data from the Meteorological Data Sharing Service System of China, the root mean square errors (RMSEs) of the GUD, HD, and MD are found to be 9.76, 5.72, and 6.98 days, respectively. Additionally, the method proposed in this article is shown to have a better extraction performance compared with several other methods. Furthermore, it is shown that, in Hebei Province, the GUD, HD, and MD are mainly affected by latitude and accumulated temperature. As the latitude increases from south to north, the GUD, HD, and MD are delayed, and for each 1° increment in latitude, the GUD, HD, and MD are delayed by 4.84, 5.79, and 6.61 days, respectively. The higher the accumulated temperature, the earlier the phenological phases occur. However, latitude and accumulated temperature have little effect on the length of the phenological phases. Additionally, the lengths of time between GUD and HD, HD and MD, and GUD and MD are stable at 46, 41, and 87 days, respectively. Overall, the proposed TWDTW method can accurately determine the key phenological phases of winter wheat at a regional scale using remote sensing time-series data. [ABSTRACT FROM AUTHOR]

With the increase in the frequency of extreme weather events in recent years, apple growing areas in the Loess Plateau frequently encounter frost during flowering. Accurately assessing the frost loss in orchards during the flowering period is of great significance for optimizing disaster prevention measures, market apple price regulation, agricultural insurance, and government subsidy programs. The previous research on orchard frost disasters is mainly focused on early risk warning. Therefore, to effectively quantify orchard frost loss, this paper proposes a frost loss assessment model constructed using meteorological and remote sensing information and applies this model to the regional-scale assessment of orchard fruit loss after frost. As an example, this article examines a frost event that occurred during the apple flowering period in Luochuan County, Northwestern China, on 17 April 2020. A multivariable linear regression (MLR) model was constructed based on the orchard planting years, the number of flowering days, and the chill accumulation before frost, as well as the minimum temperature and daily temperature difference on the day of frost. Then, the model simulation accuracy was verified using the leave-one-out cross-validation (LOOCV) method, and the coefficient of determination (R2), the root mean square error (RMSE), and the normalized root mean square error (NRMSE) were 0.69, 18.76%, and 18.76%, respectively. Additionally, the extended Fourier amplitude sensitivity test (EFAST) method was used for the sensitivity analysis of the model parameters. The results show that the simulated apple orchard fruit number reduction ratio is highly sensitive to the minimum temperature on the day of frost, and the chill accumulation and planting years before the frost, with sensitivity values of ≥0.74, ≥0.25, and ≥0.15, respectively. This research can not only assist governments in optimizing traditional orchard frost prevention measures and market price regulation but can also provide a reference for agricultural insurance companies to formulate plans for compensation after frost. [ABSTRACT FROM AUTHOR]

Rapid and accurate crop aboveground biomass estimation is beneficial for high-throughput phenotyping and site-specific field management. This study explored the utility of high-definition digital images acquired by a low-flying unmanned aerial vehicle (UAV) and ground-based hyperspectral data for improved estimates of winter wheat biomass. To extract fine textures for characterizing the variations in winter wheat canopy structure during growing seasons, we proposed a multiscale texture extraction method (Multiscale_Gabor_GLCM) that took advantages of multiscale Gabor transformation and gray-level co-occurrency matrix (GLCM) analysis. Narrowband normalized difference vegetation indices (NDVIs) involving all possible two-band combinations and continuum removal of red-edge spectra (SpeCR) were also extracted for biomass estimation. Subsequently, non-parametric linear (i.e., partial least squares regression, PLSR) and nonlinear regression (i.e., least squares support vector machine, LSSVM) analyses were conducted using the extracted spectral features, multiscale textural features and combinations thereof. The visualization technique of LSSVM was utilized to select the multiscale textures that contributed most to the biomass estimation for the first time. Compared with the best-performing NDVI (1193, 1222 nm), the SpeCR yielded higher coefficient of determination (R2), lower root mean square error (RMSE), and lower mean absolute error (MAE) for winter wheat biomass estimation and significantly alleviated the saturation problem after biomass exceeded 800 g/m2. The predictive performance of the PLSR and LSSVM regression models based on SpeCR decreased with increasing bandwidths, especially at bandwidths larger than 11 nm. Both the PLSR and LSSVM regression models based on the multiscale textures produced higher accuracies than those based on the single-scale GLCM-based textures. According to the evaluation of variable importance, the texture metrics "Mean" from different scales were determined as the most influential to winter wheat biomass. Using just 10 multiscale textures largely improved predictive performance over using all textures and achieved an accuracy comparable with using SpeCR. The LSSVM regression model based on the combination of the selected multiscale textures, and SpeCR with a bandwidth of 9 nm produced the highest estimation accuracy with R2val = 0.87, RMSEval = 119.76 g/m2, and MAEval = 91.61 g/m2. However, the combination did not significantly improve the estimation accuracy, compared to the use of SpeCR or multiscale textures only. The accuracy of the biomass predicted by the LSSVM regression models was higher than the results of the PLSR models, which demonstrated LSSVM was a potential candidate to characterize winter wheat biomass during multiple growth stages. The study suggests that multiscale textures derived from high-definition UAV-based digital images are competitive with hyperspectral features in predicting winter wheat biomass. [ABSTRACT FROM AUTHOR]

The branches of fruit trees provide support for the growth of leaves, buds, flowers, fruits, and other organs. The number and length of branches guarantee the normal growth, flowering, and fruiting of fruit trees and are thus important indicators of tree growth and yield. However, due to their low height and the high number of branches, the precise management of fruit trees lacks a theoretical basis and data support. In this paper, we introduce a method for extracting topological and structural information on fruit tree branches based on LiDAR (Light Detection and Ranging) point clouds and proved its feasibility for the study of fruit tree branches. The results show that based on Terrestrial Laser Scanning (TLS), the relative errors of branch length and number are 7.43% and 12% for first-order branches, and 16.75% and 9.67% for second-order branches. The accuracy of total branch information can reach 15.34% and 2.89%. We also evaluated the potential of backpack-LiDAR by comparing field measurements and quantitative structural models (QSMs) evaluations of 10 sample trees. This comparison shows that in addition to the first-order branch information, the information about other orders of branches is underestimated to varying degrees. The root means square error (RMSE) of the length and number of the first-order branches were 3.91 and 1.30 m, and the relative root means square error (NRMSE) was 14.62% and 11.96%, respectively. Our work represents the first automated classification of fruit tree branches, which can be used in support of precise fruit tree pruning, quantitative forecast of yield, evaluation of fruit tree growth, and the modern management of orchards. [ABSTRACT FROM AUTHOR]

Fast and accurate counting of wheat ears in field conditions is a key element for determining wheat yield. To obtain the number of wheat ears in a field, we propose a new counting algorithm based on computer vision. This algorithm counts wheat ears in remote images through semantic segmentation regression network (SSRNet). SSRNet is a multistage convolutional neural network that we propose to achieve counting problems through regression. In SSRNet, first, the original image is cropped to increase the amount of data. This method effectively solves the small sample dataset. Next, based on the cropping results, we build a fully convolutional neural network (FCNN) to segment wheat ears in field conditions. FCNN increases the accuracy of wheat ears counting by accurately segmenting wheat ears in a complex background. Then, we build a regression convolutional neural network (RCNN) to count wheat ears based on the segmentation results of FCNN. In RCNN, we propose a new activation function positive rectification linear unit (PrLU) to process the last layer of the fully connected layer, so that RCNN can effectively count the number of wheat ears in the image. Finally, a counting strategy is proposed to count the number of wheat ears in the original image. To verify the counting performance of SSRNet, we compare the counting result of SSRNet with the real value of manual statistics. The results show that the average accuracy (Acc), $R^{2}$ , and root mean squared error (RMSE) of the SSRNet count results on the test set in this article are 0.980, 0.996, and 9.437, respectively. It can be seen from the results that our proposed method can accurately count wheat ears in field conditions. At the same time, the counting time (0.11 s) shows that SSRNet can quickly estimate the number of wheat ears in field conditions. We concluded that this study can provide important technical support for the high-throughput field wheat ears counting task in large-scale phenotyping work. [ABSTRACT FROM AUTHOR]

Generating spatial crop yield information is of great significance for academic research and guiding agricultural policy. Existing public yield datasets have a coarse spatial resolution, spanning from 1 to 43 km. Although these datasets are useful for analyzing large-scale temporal and spatial change in yield, they cannot deal with small-scale spatial heterogeneity, which happens to be the most significant characteristic of the Chinese farmers' economy. Hence, we generated a 30 m Chinese winter wheat yield dataset (ChinaWheatYield30m) for major winter-wheat-producing provinces in China for the period 2016–2021 with a semi-mechanistic model (hierarchical linear model, HLM). The yield prediction model was built by considering the wheat growth status and climatic factors. It can estimate wheat yield with excellent accuracy and low cost using a combination of satellite observations and regional meteorological information (i.e., Landsat 8, Sentinel 2 and ERA5 data from the Google Earth Engine (GEE) platform). The results were validated using in situ measurements and census statistics and indicated a stable performance of the HLM based on calibration datasets across China, with a correlation coefficient (r) of 0.81 and a relative root mean square error (rRMSE) of 12.59 %. With regards to validation, the ChinaWheatYield30m dataset was highly consistent with in situ measurement data and statistical data (p<0.01), indicated by an r (rRMSE) of 0.72** (15.34 %) and 0.69** (19.16 %). The ChinaWheatYield30m is a sophisticated dataset with both high spatial resolution and excellent accuracy; such a dataset will provide basic knowledge of detailed wheat yield distribution, which can be applied for many purposes including crop production modeling and regional climate evaluation. The ChinaWheatYield30m dataset generated from this study can be downloaded from 10.5281/zenodo.7360753 (Zhao et al., 2022b). [ABSTRACT FROM AUTHOR]

The estimation of physicochemical crop parameters based on spectral indices depend strongly on planting year, cultivar, and growing period. Therefore, the efficient monitoring of crop growth and nitrogen (N) fertilizer treatment requires that we develop a generic spectral index that allows the rapid assessment of the plant nitrogen content (PNC) of crops and that is independent of year, cultivar, and growing period. Thus, to obtain the best indicator for estimating potato PNC, herein, we provide an in-depth comparative analysis of the use of hyperspectral single-band reflectance and two- and three-band spectral indices of arbitrary bands for estimating potato PNC over several years and for different cultivars and growth periods. Potato field trials under different N treatments were conducted over the years 2018 and 2019. An unmanned aerial vehicle hyperspectral remote sensing platform was used to acquire canopy reflectance data at several key potato growth periods, and six spectral transformation techniques and 12 arbitrary band combinations were constructed. From these, optimal single-, two-, and three-dimensional spectral indices were selected. Finally, each optimal spectral index was used to estimate potato PNC under different scenarios and the results were systematically evaluated based on a correlation analysis and univariate linear modeling. The results show that, although the spectral transformation technique strengthens the correlation between spectral information and potato PNC, the PNC estimation model constructed based on single-band reflectance is of limited accuracy and stability. In contrast, the optimal three-band spectral index TBI 5 (530,734,514) performs optimally, with coefficients of determination of 0.67 and 0.65, root mean square errors of 0.39 and 0.39, and normalized root mean square errors of 12.64% and 12.17% for the calibration and validation datasets, respectively. The results thus provide a reference for the rapid and efficient monitoring of PNC in large potato fields. [ABSTRACT FROM AUTHOR]

Aboveground biomass (AGB) is an important indicator for crop-growth monitoring and yield prediction, and accurate monitoring of AGB is beneficial to agricultural fertilization management and optimization of planting patterns. Imaging spectrometer sensors mounted on unmanned aerial vehicle (UAV) remote-sensing platforms have become an important technical method for monitoring AGB because the method is convenient, rapidly collects data and provides image data with high spatial and spectral resolution. To confirm the feasibility of UAV hyperspectral remote-sensing technology to estimate AGB, this study acquired hyperspectral images and measured AGB data over the potato bud, tuber formation, tuber growth, and starch-storage periods. The canopy spectrum obtained in each growth period was smoothed by using the Savitzky–Golay filtering method, and the spectral-reflection feature parameters, spectral-location feature parameters, and vegetation indexes were extracted. First, a Pearson correlation analysis was performed between the three types of characteristic spectral parameters and AGB, and the spectral parameters that reached a significant level of 0.01 in each growth period were selected. Next, the spectral parameters reaching a significance of 0.01 were optimized and screened by moving window partial least squares (MWPLS), Monte Carlo uninformative variable elimination (MC-UVE), and random frog (RF) methods, and the final model parameters were determined according to the thresholds of the root mean square error of cross-validation (RMSEcv), the reliability index, and the selected probability. Finally, the three optimal characteristic spectral parameters and their combinations were used to estimate the potato AGB in each growth period by combining the partial least squares regression (PLSR) and Gaussian process regression (GPR) methods. The results show that, (i) ranked from high to low, vegetation indexes, spectral-location feature parameters, and spectral-reflection feature parameters in each growth period are correlated with the AGB, and these correlations all first improve and then degrade in going from the budding period to the starch-storage period. (ii) The AGB estimation model based on the characteristic variables screened by the three methods in each growth period is most accurate with RF, less so with MC-UVE, and least accurate with MWPLS. (iii) Estimating the AGB with the same variables combined with the PLSR method in each growth period is more accurate than the corresponding GPR method, but the estimations produced by the two methods both show a trend of first improving and then worsening from the budding period to the starch-accumulation period. The accuracy of the estimation models constructed by PLSR and GPR from high to low is based on comprehensive variables, vegetation indexes, spectral-location feature parameters and spectral-reflection feature parameters. (iv) When combined with the RF-PLSR method to estimate AGB in each growth period, the best R2 values are 0.65, 0.68, 0.72, and 0.67, the corresponding RMSE values are 167.76, 162.98, 160.77, and 169.24 kg/hm2, and the corresponding NRMSE values are 19.76%, 16.01%, 15.04%, and 16.84%. The results of this study show that a variety of characteristic spectral parameters may be extracted from UAV hyperspectral images, that the RF method may be used for optimizing and screening, and that PLSR regression provides accurate estimates of the potato AGB. The proposed approach thus provides a rapid, accurate, and nondestructive way to monitor the growth status of potatoes. [ABSTRACT FROM AUTHOR]

• We developed an RRCNet that fuses the deep and shallow spectral features to map RRC. • RRCNet was trained based on an RTM-based "soil-crop-crop residue" spectra dataset. • Use transfer learning (TL) to recalibrate the fully connected layer of RRCNet. • RRCNet + TL can provide more accurate RRC estimates than SI + statistical regression. Accurate determination of rice residue cover (RRC) can improve the monitoring of tillage information. Currently, the accurate determination of RRC using optical remote sensing is hindered by variations in cropland moisture and cover of following crops. The fractional cover (FC) of the soil (f S), crop (f C), and crop residue (f CR) changes (f S + f C + f CR = 1) after the following crop is planted, which increases the difficulty of remote-sensing RRC estimation. Cropland soil moisture and crop residue moisture affect the values of cropland and crop residue spectral indices (CRSIs), thereby reducing the accuracy of remote-sensing RRC estimation. Deep learning techniques (e.g., convolutional neural networks [CNN] and transfer learning [TL]) have been proven to extract the deep features of input images with distortion invariance, such as displacement and scaling, which are similar to moisture and the following crop effects on remote-sensing CRSIs. This study aimed to evaluate the combined use of deep features of cropland spectra extracted by deep learning techniques to estimate the cropland RRC under the effects of variations in cropland moisture and cover of the following crops. This study proposes an RRCNet CNN that fuses deep and shallow features to improve RRC estimation. A PROSAIL radiative transfer model was employed to simulate a cropland "soil–crop–crop residue" mixed spectra dataset (n = 103,068), considering the variations in cropland moisture and the cover of the following crop. The RRCNet was first pre-trained using the simulated dataset, and then the knowledge from the pre-trained RRCNet was updated based on field experimental FCs, RRCs, and Sentinel-2 image spectra using the TL technique. Our study indicates that RRCNet can incorporate shallow and deep spectral features of cropland "soil–crop–crop residue" mixed spectra, providing high-performance FCs and RRC estimation. The FCs and RRC estimates from RRCNet + TL (FCs: R 2 = 0.939, root mean squared error (RMSE) = 0.071; RRC: R 2 = 0.891, RMSE = 0.083) were more accurate than those from CRSI + multiple linear regression, CRSI + random forest, and CRSI + support vector machine (FCs: R 2 = 0.877–907, RMSE = 0.086–0.101; RRC: R 2 = 0.378–0.714, RMSE = 0.133–0.229). We mapped the multistage RRC based on Sentinel-2 multispectral instrument (MSI) images and RRCNet. Tillage information can be inferred from RRC and RRC difference maps changes. [ABSTRACT FROM AUTHOR]

Vertical heterogeneity of the biochemical characteristics of crop canopy is important in diagnosing and monitoring nutrition, disease, and crop yield via remote sensing. However, the research on vertical isomerism was not comprehensive. Experiments were carried out from the two levels of simulation and verification to analyze the applicability of this recently development model. Effects of winter wheat on spectrum were studied when input different structure parameters (e.g., leaf area index (LAI)) and physicochemical parameters (e.g., chlorophyll content ( C h l a + b ) and water content ( C w )) to the mSCOPE (Soil Canopy Observation, Photochemistry, and Energy fluxes) model. The maximum operating efficiency was 127.43, when the winter wheat was stratified into three layers. Meanwhile, the simulation results also proved that: the vertical profile of LAI had an influence on canopy reflectance in almost all bands; the vertical profile of C h l a + b mainly affected the reflectivity of visible region; the vertical profile of C w only affected the near-infrared reflectance. The verification results showed that the vegetation indexes (VIs) selected of different bands were strongly correlated with the parameters of the canopy. LAI, C h l a + b and C w affected VIs estimation related to LAI, C h l a + b and C w respectively. The Root Mean Square Error (RMSE) of the new-proposed N D V I g r e e n was the smallest, which was 0.05. Sensitivity analysis showed that the spectrum was more sensitive to changes in upper layer parameters, which verified the rationality of mSCOPE model in explaining the law that light penetration in vertical nonuniform canopy gradually decreases with the increase of layers. [ABSTRACT FROM AUTHOR]

Yield evaluation of breeding lines is the key to successful release of cultivars, which is becoming a serious issue due to soil heterogeneity in enlarged field tests. This study aimed at establishing plot-yield prediction models using unmanned aerial vehicle (UAV)-based hyperspectral remote sensing for yield-selection in large-scale soybean breeding programs. Three sets of soybean breeding lines (1103 in total) were tested in blocks-in-replication experiments for plot yield and canopy spectral reflectance on 454~950 nm bands at different growth stages using a UAV-based hyperspectral spectrometer (Cubert UHD185 Firefly). The four elements for plot-yield prediction model construction were studied respectively and concluded as: the suitable reflectance-sampling unit-size in a plot was its 20%–80% central part; normalized difference vegetation index (NDVI) and ration vegetation index (RVI) were the best combination of vegetation indices; the initial seed-filling stage (R5) was the best for a single stage prediction, while another was the best combination for a two growth-stage prediction; and multi-variate linear regression was suitable for plot-yield prediction. In model establishment for each material-set, a random half was used for modelling and another half for verification. Twenty-one two growth-stage two vegetation-index prediction models were established and compared for their modelling coefficient of determination (RM2) and root mean square error of the model (RMSEM), verification RV2 and RMSEV, and their sum RS2 and RMSES. Integrated with the coincidence rate between the model predicted and the practical yield-selection results, the models, MA1-2, MA4-2 and MA6-2, with coincidence rates of 56.8%, 58.5% and 52.4%, respectively, were chosen for yield-prediction in yield-test nurseries. The established model construction elements and methods can be used as local models for pre-harvest yield-selection and post-harvest integrated yield-selection in advanced breeding nurseries as well as yield potential prediction in plant-derived-line nurseries. Furthermore, multiple models can be used jointly for plot-yield prediction in soybean breeding programs. [ABSTRACT FROM AUTHOR]

Timely and effective prediction of nitrogen content in summer maize could provide support data for precise fertilization. In this study, the feasibility and expansibility of predicting the nitrogen mechanism model of summer maize leaves through its entire growth period were investigated on the basis of the theory of leaf radiation transmission mechanism. A complete random test of data from two maize varieties and two nitrogen fertilizer applications in 2017 was conducted. Three versions of the leaf optical PROperties SPECTra (PROSPECT) model, namely, PROSPECT-4, PROSPECT-5, and PROSPECT-D were used to link the established leaf nitrogen density (LND) and chlorophyll-a + b (chl-a + b) models, that is, chl-a + b-LND model. A nitrogen response transfer model (N-RTM) was established by linking the optimal PROSPECT and chl-a + b-LND models. Results were as follows. (1) chl-a + b estimation using the PROSPECT-D model yielded the highest accuracy (the coefficient of determination (R2) = 0.774, the normalized root mean squared error (nRMSE) = 13.19%) among the three PROSPECT models, it shows that the model considering more factors can better reflect the internal law of blade, and could be used as the basic model of N-RTM; (2) Established chl-a + b-LND models based on the dataset from each growth stage showed differences using the confidence interval method, and the R2 values of the optimal regression model at V12, VT, and R3 were 0.794, 0.781, and 0.821, respectively. Based on the changes of chl-a + b and LND during the growth period, a piecewise model was constructed; (3) The R2 and nRMSE values between the measured and estimated LNDs were 0.656% and 22.86%, respectively. The validation results are better than the traditional empirical model. The results showed that the segmented model, which considered the interaction of various factors within the leaves and the change of chl-a + b-LND during the growth period, had better performance in nitrogen monitoring. The constructed nitrogen model in this study preliminarily realized the remote sensing prediction of the nitrogen mechanism model and had a certain mechanism. [ABSTRACT FROM AUTHOR]

• Our method focuses on exploring the effect of texture indices to estimate potato PNC. • An texture index was developed for PNC estimation throughout the potato growth cycle. • O-TTI3 can further improve the estimation accuracy of potato PNC. Plant nitrogen content (PNC) is vital for evaluating crop nitrogen nutrient status and for net primary productivity. Therefore, rapid and accurate acquisition of crop PNC information can help reduce fertilizer and increase efficiency in modern agriculture. This study investigated whether an optimized texture index constructed from hyperspectral characteristic wavelength images acquired from a unmanned aerial vehicle (UAV) may be used to estimate PNC in potatoes over multiple growth stages. A potato field trial conducted in 2019 in Beijing, China, included different planting densities, nitrogen fertilizer gradients, and potato varieties with three replicates. A UAV remote sensing platform served to acquire hyperspectral images during three critical N demand periods for potatoes. In addition, we simultaneously conducted field sampling to obtain ground-truth PNC measurements. Following the classical form of vegetation indices, 12 texture indices were constructed based on hyperspectral texture features, and an arbitrary variable combination optimization algorithm was used to optimize these indices. Finally, the texture index, which had the highest correlation with potato PNC, was used as the best indicator for estimating the PNC status of potato at multiple growth stages and whether this indicator, in combination with spectral information, could further improve the accuracy of potato PNC estimation was subsequently explored. The results showed that (i) the optimal texture index TTI3 (R 494 -Cor, R 578 -Hom, R 514 -Sem) maintained a good linear relationship with PNC at the late stage of potato growth, and the accuracy and stability of the PNC estimation models constructed based on it was significantly better than that of a single texture metric. (ii) Compared with spectral information alone, the texture index combined with spectral features improved the accuracy of potato PNC estimation. More specifically, TTI3 (R 494 -Cor, R 578 -Hom, R 514 -Sem) combined with the three-band spectral index TBI 5 (530, 734, 514) achieved the best estimation accuracy with calibrated R2 , root mean square error (RMSE), and normalized RMSE of 0.77%, 0.28%, and 9.88%, respectively. The results of this study showed that the texture index constructed by combining multiple texture metrics enhanced the association between texture features and potato PNC over multiple growth stages, thus improving the monitoring accuracy of potato nitrogen nutrition status. This study can provide a reference for accurately managing crop nitrogen nutrition. [ABSTRACT FROM AUTHOR]

Grain protein content (GPC) is an important indicator of wheat quality. Earlier estimation of wheat GPC based on remote sensing provided effective decision to adapt optimized strategies for grain harvest, which is of great significance for agricultural production. The objectives of this field study are: (i) To assess the ability of spectral vegetation indices (VIs) of Sentinel 2 data to detect the wheat nitrogen (N) attributes related to the grain quality of winter wheat production, and (ii) to examine the accuracy of wheat N status and GPC estimation models based on different VIs and wheat nitrogen parameters across Analytical Spectra Devices (ASD) and Unmanned Aerial Vehicle (UAV) hyper-spectral data-simulated sentinel data and the real Sentinel-2 data. In this study, four nitrogen parameters at the wheat anthesis stage, including plant nitrogen accumulation (PNA), plant nitrogen content (PNC), leaf nitrogen accumulation (LNA), and leaf nitrogen content (LNC), were evaluated for their relationship between spectral parameters and GPC. Then, a multivariate linear regression method was used to establish the wheat nitrogen and GPC estimation model through simulated Sentinel-2A VIs. The coefficients of determination ( R 2 ) of four nitrogen parameter models were all greater than 0.7. The minimum R 2 of the prediction model of wheat GPC constructed by four nitrogen parameters combined with VIs was 0.428 and the highest R 2 was 0.467. The normalized root mean square error (nRMSE) of the four nitrogen estimation models ranged from 26.333% to 29.530% when verified by the ground-measured data collected from the Beijing suburbs, and the corresponding nRMSE for the GPC-predicted models ranged from 17.457% to 52.518%. The accuracy of the estimated model was verified by UAV hyper-spectral data which had resized to different spatial resolution collected from the National Experimental Station for Precision Agriculture. The normalized root mean square error (nRMSE) of the four nitrogen estimation models ranged from 16.9% to 37.8%, and the corresponding nRMSE for the GPC-predicted models ranged from 12.3% to 13.2%. The relevant models were also verified by Sentinel-2A data collected in 2018 while the minimum nRMSE for GPC invert model based on PNA was 7.89% and the maximum nRMSE of the GPC model based on LNC was 12.46% in Renqiu district, Hebei province. The nRMSE for the wheat nitrogen estimation model ranged from 23.200% to 42.790% for LNC and PNC. These data demonstrate that freely available Sentinel-2 imagery can be used as an important data source for wheat nutrition and grain quality monitoring. [ABSTRACT FROM AUTHOR]

• Potential of land surface temperature to predict flowering phenology was explored. • A spatio-temporal reconstruction approach is developed to compenstate the missing LST. • The predicted first flowering date are validated with field measurements at multiple sites. • Feasibility for generating spatially continuous information on flowering phenology. The first flowering date (FFD) is a critical phenological parameter closely related to the apple yield, so the accurate prediction of the FFD is important for precise orchard production management. Existing methods to predict the FFD are mostly based on air temperature (Ta) measured at meteorological stations, but to great differences in meteorological variations and the ecological conditions, these methods cannot present the differences of FFD under complex meteorological conditions and provide spatially continuous FFD information at the level of a region. Therefore, we propose a method to predict spatially continuous apple FFD from remote sensing land surface temperature (LST) based on flowering prediction model. Firstly, the missing LST data were reconstructed by using spatio-temporal reconstruction (STR) approach developed. Next, new air temperature (NAT) data were generated by using the daily Ta estimation (DTE) model and the reconstructed LST. Finally, apple FFD was predicted by the NAT data and the apple flowering prediction model established based on random forest (RF) algorithm and the phenology sequential model, and the prediction accuracy was verified by comparison with the independently measured apple FFD. The LST reconstructed by using the STR approach has mean absolute error (MAE) ranging from 0.51 to 0.68 °C, and root mean square error (RMSE) ranging from 1.07 to 1.21 °C. The MAE between the NAT data and the High-Resolution Land Surface Data Assimilation System (HR-CLDAS) meteorological data ranges from 2.15 to 3.23 °C, and the RMSE ranges from 2.81 to 4.27 °C. In addition, the determination coefficient (R2) and RMSE between the predicted and measured FFD is 0.72 and 2.96 days, respectively. These results demonstrate that the developed method maximizes the potential of MODIS LST in predicting spatially continuous apple FFD, which is valuable for flower and fruit thinning, to defend against frost disasters, and in general for refined orchard production management. [ABSTRACT FROM AUTHOR]

• Develop a within-growing season yield forecast system with random forest model. • Random forest model performs well in predicting grain yield in China. • We identified the most important stage-specific predictors determining crop yield. • The most important variable influencing yields varied with crop types. Accurate and timely crop yield forecasts can provide essential information to make conclusive agricultural policies and to conduct investments. Recent studies have used different machine learning techniques to develop such yield forecast systems for single crops at regional scales. However, no study has used multiple sources of environmental predictors (climate, soil, and vegetation) to forecast yields for three major crops in China. In this study, we adopted 7-year observed crop yield data (2013–2019) for three major grain crops (wheat, maize, and rice) across China, and three major data sets including climate, vegetation indices, and soil properties were used to develop a dynamic yield forecasting system based on the random forest (RF) model. The RF model showed good performance for estimating yields of all three crops with correlation coefficient (r) higher than 0.75 and normalized root means square errors (nRMSE) lower than 18.0%. Our results also showed that crop yields can be satisfactorily forecasted at one to three months prior to harvest. The optimum lead time for yield forecasting depended on crop types. In addition, we found the major predictors influencing crop yield varied between crops. In general, solar radiation and vegetation indices (especially during jointing to milk development stages) were identified as the main predictor for winter wheat; vegetation indices (throughout the growing season) and drought (especially during emergence to tasseling stages) were the most important predictors for spring maize; soil moisture (throughout the growing season) was the dominant predictor for summer maize, late rice, and mid rice; precipitation (especially during booting to heading stages) was the main predictor for early rice. Our study provides insights into practical crop yield forecasting and the understanding of yield response to environmental conditions at a large scale across China. The methods undertaken in this research can be easily implemented in other countries with available information on climate, soil, and vegetation conditions. [ABSTRACT FROM AUTHOR]

Industrialization production with high quality and effect on winter is an important measure for accelerating the shift from increasing agricultural production to improving quality in terms of grain protein content (GPC). Remote sensing technology achieved the GPC prediction. However, large deviations in interannual expansion and regional transfer still exist. The present experiment was carried out in wheat producing areas of Beijing (BJ), Renqiu (RQ), Quzhou, and Jinzhou in Hebei Province. First, the spectral consistency of Landsat 8 Operational Land Imager (LS8) and RapidEye (RE) was compared with Sentinel-2 (S2) satellites at the same ground point in the same period. The GPC prediction model was constructed by coupling the vegetation index with the meteorological data obtained by the European Center for Medium-range Weather Forecasts using hierarchical linear model (HLM) method. The prediction and spatial expansion of regional GPC were validated. Results were as follows: (1) Spectral information calculated from S2 imagery were highly consistent with LS8 (R2 = 1.00) and RE (R2 = 0.99) imagery, which could be jointly used for GPC modeling. (2) The predicted GPC by using the HLM method (R2 = 0.524) demonstrated higher accuracy than the empirical linear model (R2 = 0.286) and showed higher improvements across inter-annual and regional scales. (3) The GPC prediction results of the verification samples in RQ, BJ, Xiaotangshan (XTS) in 2018, and XTS in 2019 were ideal with root mean square errors of 0.61%, 1.13%, 0.91%, and 0.38%, and relative root mean square error of 4.11%, 6.83%, 6.41%, and 2.58%, respectively. This study has great application potential for regional and inter-annual quality prediction. [ABSTRACT FROM AUTHOR]

Accurate and dynamic monitoring of crop nitrogen status is the basis of scientific decisions regarding fertilization. In this study, we compared and analyzed three types of spectral variables: Sensitive spectral bands, the position of spectral features, and typical hyperspectral vegetation indices. First, the Savitzky-Golay technique was used to smooth the original spectrum, following which three types of spectral parameters describing crop spectral characteristics were extracted. Next, the successive projections algorithm (SPA) was adopted to screen out the sensitive variable set from each type of parameters. Finally, partial least squares (PLS) regression and random forest (RF) algorithms were used to comprehensively compare and analyze the performance of different types of spectral variables for estimating corn leaf nitrogen content (LNC). The results show that the integrated variable set composed of the optimal ones screened by SPA from three types of variables had the best performance for LNC estimation by the validation data set, with the values of R2, root means square error (RMSE), and normalized root mean square error (NRMSE) of 0.77, 0.31, and 17.1%, and 0.55, 0.43, and 23.9% from PLS and RF, respectively. It indicates that the PLS model with optimally multitype spectral variables can provide better fits and be a more effective tool for evaluating corn LNC. [ABSTRACT FROM AUTHOR]