Next Article in Journal
Effectiveness of the Reconstructed MODIS Typical-Angle Reflectances on Forest Biomass Estimation
Next Article in Special Issue
Using Apparent Electrical Conductivity to Delineate Field Variation in an Agroforestry System in the Ozark Highlands
Previous Article in Journal
Preparatory Railway Track Geometry Estimation Based on GNSS and IMU Systems
Previous Article in Special Issue
Potential of ALOS2 Polarimetric Imagery to Support Management of Poplar Plantations in Northern Italy
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transfer-Learning-Based Approach for Yield Prediction of Winter Wheat from Planet Data and SAFY Model

1
National Engineering and Technology Center for Information Agriculture, Nanjing Agricultural University, Nanjing 210095, China
2
Key Laboratory of Quantitative Remote Sensing in Agriculture of Ministry of Agriculture and Rural Affairs, Information Technology Research Center, Beijing Academy of Agriculture and Forestry Sciences, Beijing 100097, China
3
College of Geomatics, Shandong University of Science and Technology, Qingdao 266590, China
4
College of Agriculture, Qingdao Hengxing University of Science and Technology, Qingdao 266100, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2022, 14(21), 5474; https://doi.org/10.3390/rs14215474
Submission received: 14 September 2022 / Revised: 13 October 2022 / Accepted: 24 October 2022 / Published: 31 October 2022
(This article belongs to the Special Issue Application of Remote Sensing in Agroforestry II)

Abstract

:
Crop production is one of the major interactions between humans and the natural environment, in the process, carbon is translocated cyclically inside the ecosystem. Data assimilation algorithm has advantages in mechanism and robustness in yield estimation, however, the computational efficiency is still a major obstacle for widespread application. To address the issue, a novel hybrid method based on the combination of the Crop Biomass Algorithm of Wheat (CBA-Wheat) to the Simple Algorithm For Yield (SAFY) model and the transfer learning method was proposed in this paper, which enables winter wheat yield estimation with acceptable accuracy and calculation efficiency. The transfer learning techniques learn the knowledge from the SAFY model and then use the knowledge to predict wheat yield. The main results showed that: (1) The comparison using CBA-Wheat between measured AGB and predicted AGB all reveal a good correlation with R2 of 0.83 and RMSE of 1.91 t ha−1, respectively; (2) The performance of yield prediction was as follows: transfer learning method (R2 of 0.64, RMSE of 1.05 t ha−1) and data assimilation (R2 of 0.64, RMSE of 1.01 t ha−1). At the farm scale, the two yield estimation models are still similar in performance with RMSE of 1.33 t ha−1 for data assimilation and 1.13 t ha−1 for transfer learning; (3) The time consumption of transfer learning with complete simulation data set is significantly lower than that of the other two yield estimation tests. The number of pixels to be simulated was about 16,000, and the computational efficiency of the data assimilation algorithm and transfer learning without complete simulation datasets. The transfer learning model shows great potential in improving the efficiency of production estimation.

1. Introduction

Crop yield is of great importance to food security and directly relates to the economic benefits of farmers [1]. It is therefore an imperative task for robust large-scale yield estimation in agricultural research and application [2].
With the advantages of non-destructive, high throughput and spatially continuous observation, remote sensing (RS) is a popular method for monitoring bio-physical parameters [3,4,5,6]. Over the past decades, the application of RS information obtained by optical and radar sensors on the terrestrial, aerial and satellite platforms in estimating crop production have considerably increased, including empirical, semi-empirical and mechanism models. Compared with traditional image processing tools, the Google Earth Engine (GEE) platform facilitates RS data acquisition, RS image processing, and analysis [7]. The GEE platform provides favorable conditions for crop yield prediction, crop growth monitoring and crop classification on a large regional scale [8,9,10]. For simplicity and computational efficiency, the majority of studies on crop yield prediction have been applied for crop yield estimation using vegetation indices (VIs) [11]. Some of the vegetation indices, i.e., the Normalized Difference Vegetation Index (NDVI) [12], the Enhanced Vegetation Index (EVI) [13] and the Optimized soil-adjusted vegetation indices (OSAVI) [14], are based on the Red-NIR isolines of the electromagnetic spectrum, prompted by the motivation to predict bio-physical parameters using bands available from different remote sensing platforms. With the future application of RS data in agriculture, multi-band VI, multi-calculation form and radar VIs are constructed [15]. The simplest method to estimate crop yield is to find the statistical relationships between only VIs and crop yield at specific phenological stages or multiple phenological stages, which could express acceptable crop yield estimation accuracy and performance under specific regions and years. The booting stage and flowering stage are often considered be the most suitable period to predict crop yield using optical remote sensing data [6,16]. Poor Spatial-temporal extrapolation without recalibration and revalidation confines the further application of these methods. The yield formation process is affected by many factors, such as soil, meteorology and genes [17,18,19]. Typically, this is now done by using different methods, e.g., hierarchical linear model [6] and machine learning methods [5,20], to fuse RS data with information related to crop growth. By integrating hyperspectral data and meteorological data, Li et al. [6] developed hierarchical linear modeling (HLM) to solve the interannual expandable problem of the wheat yield prediction model. To better explain the response of crop yield to changes in sub-seasonal environmental factors across the entire growing season, many machine learning methods and phenological information were considered for yield estimation models [5,21,22]. Despite the excellent performance of a semi-empirical, a well-trained model depends on limited mechanisms of crop growth that still have defects in expressing the process of crop yield formation, and their generality for multiple regions has not been thoroughly validated.
The process-oriented crop growth model (CGM) is a verified tool to simulate crop development and crop yield formulation at a farm scale with detailed input data, e.g., crop varieties, agronomic management practices, soil, and climate factors [22,23,24]. As a bridge between CGMs and RS data, data assimilation technology combined with the advantages of RS data and CGMs has been recognized as a reliable way to improve crop yield simulations over a large area. Crop growth variables are selected for assimilation to adjust the CGM simulation and thus obtain updated crop yield predictions. Among all the variables available for assimilation, leaf area index, aboveground dry biomass, canopy cover and soil moisture have been popularly elected in the assimilation system [25,26,27,28]. Aboveground dry biomass (AGB) is strongly related to crop yield, consequently, statistical regression, machine learning method and radiation transfer model have all been used to monitor AGB. Active and passive RS often has the ability to better predict AGB at low coverage, but all face the difficulty of saturation at high biomass [29]. Multi-source input variables are used in AGB prediction models, including optical VIs and SAR indicators, optical VIs and imagery textures, and optical VIs and phenological information [30,31,32]. Therefore, an in-depth study on agronomic parameters prediction using RS data is of great significance to further enhance crop yield prediction performance in the data assimilation system. In general, the estimation of agronomic parameters using the combination of multi-source data is more effective than that of only single-source data [33]. Different CGMs differ in their structure, parameterization and complexity, and may obtain different results even if the same observation data is applied [22,23,24]. Therefore, the data assimilation combined only one crop growth model, a yield estimation assimilation system using CGMs have also been constructed and achieved acceptable performance [34,35]. Although crop yield estimation using data assimilation has made great progress, problems remained. Primarily, for the periods of state variables extraction, they were not well represented, previously used state variables cannot represent the environmental factors resulting in the lack of universality in crop yield estimation under different environments and conditions. Second, for the assimilation algorithms, the computational intensity was always a problem that could not be ignored when estimating production in a large area. Summarizing the existing issues above, the existing problem should be using sufficient data which can be mined from simulation results of CGMs and RS data to complete the estimation of crop pre-harvest yield with less calculation.
Recently, machine learning methods have been used to learn abstract features from high-spatial-resolution RS images for yield estimation [20]. Despite the excellent performance of machine learning methods, a trained network depends on the large number of labeled samples that are difficult and expensive to collect in agricultural applications. Therefore, an economical approach such as transfer learning techniques would be desirable for crop yield prediction for different environmental conditions. To date, there has been limited work done to construct transfer learning techniques to improve local crop yield predictions. Therefore, the major objectives were as follows in this study: (1) to explore dynamic patterns of aboveground dry biomass (AGB) observations and simulated crop yield to present a transfer learning yield inversion model; (2) to examine the performance and transferability analysis of transfer learning yield inversion model cross-year yield estimation; (3) test the computational efficiency and accuracy of the data assimilation system and the proposed method.

2. Materials and Methods

2.1. Study Location and Field Data Collection

Field experiments (40.17°N, 116.43°E) were conducted in 2017~2018 and 2018~2019 at National Precision Agriculture Research Center near Beijing, China (see Figure 1). The historical experiments involve different cultivars, nitrogen management and growing seasons [32]. Two experiments (Calibration sets), 2017–2018 (Exp. 1) and 2018–2019 (Exp. 2) were designed as completely randomized designs with cultivars (Jingdong 18 and Lunxuan 167) and four nitrogen levels (N1: 0 kg N ha−1; N2: 90 kg N ha−1; N3: 180 kg N ha−1, and N4: 270 kg N ha−1). There were 32 plots per growing season, with 128 biomass samples (32 plots × 4 periods) and 32 yields collected annually, see Li et al. [32] for details. At each sampling site, Zadok’s growth stage [36] was investigated. Additionally, a yield survey was performed inside the farm in 2018 and 2019, 32 and 71 yield samples were randomly collected respectively for verification of the yield estimation model (Figure 1b).

2.2. Data Collection

2.2.1. Satellite Imagery Data Acquisition and Preprocessing

Planet satellite (www.planet.com, accessed on 1 September 2022) is an emerging remote sensing satellite in 2018~2019, with high frequency, good image quality and high data coverage efficiency. There are more than 170 satellites in the Planet’s small satellite constellation. Hundreds of satellites independently photograph the world every day, which can achieve daily global coverage. The Planet satellite imagery parameters are shown in Table 1. The images used in this paper are 3B-level data products after sensor calibration, radiometric calibration, orthorectification and atmospheric correction. The detailed process of satellite data processing is shown in Bai et al. [37].
A commonly used vegetation index from literature, the enhanced vegetation index 2 (EVI2) [38], was used to predict AGB in this study. The formula is as follows:
EVI 2 = 2.5 × ( NIR R ) / ( NIR + 2.4 × R + 1 )
where 2.5, 2.4 and 1 are the adjustment factors to reduce the impact of soil, atmosphere and saturation on the prediction of agricultural parameters. NIR represents the reflectance of the near-infrared band and R represents the reflectance of the red band.

2.2.2. Meteorological Data

The meteorological data from the ERA5 dataset (ECMWF, http://www.ecmwf.int, accessed on 1 September 2022) was used in this study, including daily minimum temperatures (Tmin), maximum temperatures (Tmax) and solar radiation, were acquired with a resolution of 0.125° for the experimental site. The daily meteorological data in this study are shown in Figure 2.

2.2.3. Field Measured AGB

In each growth period, the wheat tillers number per unit area was investigated in each experimental plot. Destructive sampling of 20 randomly selected tillers were taken to a laboratory and 20 tillers samples in each sampling point were then oven-dried to a constant weight. The wheat AGB (t ha−1) was calculated as follows:
AGB = W × T S × A × 10 2
where S represents the number of tillers in the destructive samples, W represents the weight (g) of the destructive samples, A represents the sampled area (m2), and T represents the total number of counted tillers in the sampled area (1 m2 in this study). The 10−2 represent the coefficient transformed from g m−1 to t ha−1.

2.2.4. Winter Wheat Yield Observations

A 1 m2 area of the winter wheat was destructive sampled to measure the harvested yield. Each sample point was in the center area of the experimental plot to avoid edge effects. The harvested grain had a 14% moisture content was threshed, air dried, and weighed on an electronic balance.

2.3. Data Analysis

In this study, one transfer learning method for yield prediction based on the Simple Algorithm For Yield (SAFY) model (see Section 2.3.1) [39], and the deep neural network (DNN) [40] method was constructed. Data assimilation using Shuffled complex evolution with PCA (SP-UCI) [41,42] by integrating the remote sensed AGB and SAFY model was tested in comparison to the capacity of the constructed model (transfer learning method). A detailed flowchart for the yield estimation method is presented in Figure 3.

2.3.1. Principle of SAFY Model

SAFY, a simple crop growth model, is used for simulating AGB on a daily step and is easy to integrate RS data for estimating AGB and crop yield [3,43,44]. During the AGB increased period, daily AGB production is based on light-use efficiency (LUE) theory, depending on effective LUE (ELUE), photosynthetically active radiation absorbed by crop canopy (APAR), daily average temperature (Ta) and temperature-stress-function (FT) [40]. Which uses the following equation:
Δ AGB = ELUE × F T ( Ta ) × APAR
APAR = ε 1 × ε c × Rg
where, ε 1 ,   ε c and Rg represent light-interception efficiency, climatic efficiency and daily incoming global variation. Where, ε 1 depends on light-interception coefficient (k) and green leaf area index (GLAI) based on Beer’s law:
ε 1 = 1 e k × GLAI
The daily maximum and minimum temperature directly affect the daily AGB accumulation rate, which is accounted for by Ta in the SAFY model determined by minimal temperature (Tmin), maximal temperature (Tmax), and optical temperature (Topt) for growth. Which uses the following equation:
F T ( Ta ) = { 1 [ Topt Ta Topt Tmin ] 2 ,   if   Tmin < Ta < Topt 1 [ Ta Topt Tmax Topt ] 2 ,   if   Topt < Ta < Tmax 0 ,   if   Ta < Tmin   or   Ta > Tmax
The dynamics of GLAI can be split into two phenological phases, leaf extent expansion and leaf disappearance during crop senescence, which is determined by the sum of temperature (STT) after crop emergence. In the leaf expansion stage, both AGB and GLAI are in the growth stage, a portion of the daily accumulated AGB by allocation to leaf function (Pl) was partitioned to dry leaf biomass, and then the increase of daily dry leaf biomass is converted in increase of daily leaf area (ΔGAI+) based on specific leaf area (SLA). If Pl > 0, ΔGAI+ leads to:
Δ GAI + = Δ AGB Pl SLA
Pl = 1 Pl a e Pl b SMT
where Pla and Plb represent two partition-to-leaf parameters and SMT represents temperature accumulation. When the SMT reached the given threshold STT, the leaves start to senesce at a constant rate (Rs). If SMT > STT, it leads to:
Δ GAI = GLAI × (   Ta STT ) / Rs
During the grain filling phase, the crop yield is proportional to the maximal AGB, with a constant harvest index (HI) partitioned to grains. It leads to:
Yield = HI × AGB max
Parameters of the SAFY model such as GLAI and ELUE in the model indirectly characterize the simulated yield differences caused by other agro-environmental stress. The description and range of parameters of the SAFY model are as follows (see Table 2):

2.3.2. Predicted AGB from Planet Imageries Using the CBA-Wheat Model

Most existing AGB monitoring models perform well only at a specific growth stage at a single location, with poor transfer during different growing stages during the season [28,43]. AGB is an important feature of the SAFY model, and excellent AGB estimation is of great significance for transfer learning to learn the characteristics of the SAFY model for yield estimation. The newly developed crop biomass algorithm for wheat (CBA-Wheat) [28] was chosen to predict winter AGB over the growing season, with a two-level piecewise model, such as:
AGB = α × EVI 2 + β
α = 1.27 e p 1 × ZS
β = p 2 × ZS p 3
where α and β represent the CBA-Wheat coefficients, ZS represent Zadok’s growth stage [37]. p1, p2 and p3 are the CBA-Wheat models constructed based on Li et al. [32]. Although this paper uses the same sets of ground data as Li et al. [32], considering the differences between the Planet satellite and ground hyperspectral data, the coefficients are fine-tuned.

2.3.3. Transfer Learning Method

The general framework of the transfer learning method (see Figure 3a) in this study is to summarize a crop growth model (SAFY) to use the information in the simulated labeled AGB-yield datasets to estimate the wheat yield and improve the efficiency of yield estimation. The transferability of features extracted by deep neural networks has been demonstrated and applied to natural language processing, computer vision and other fields [47,48,49]. This study then chose root mean squared error (RMSE) and loss value to evaluate the model accuracy. The implementation steps of the transfer learning method are as follows:
(I).
SAFY parameter sets construction: 200,000 sets of parameter combinations of SAFY were generated based on Monte Carlo (MC) algorithm.
(II).
Construction of possible AGB datasets and yield datasets based on the SAFY model: The parameter set constructed in step I was input into the SAFY model to obtain the possible AGB datasets and yield datasets. The time efficiency test of transfer learning is divided into two types with or without simulated data sets: no simulated datasets and with simulated datasets.
(III).
Train DNN model: A four-layer fully-connected network is constructed to train simulated AGB-yield from the SAFY model. The model was pre-trained using the AGB datasets and yield datasets simulated in step II, and fine-tuned by the measured data using transfer learning.
(IV).
Forecast yield based on transfer learning method: The AGB predicted from the CBA-Wheat model is utilized as the input layer of the trained DNN to predict winter wheat yield.

2.3.4. Data Assimilation

In this section, the Shuffled complex evolution with PCA (SP-UCI) algorithm is improved based on the SCE-UA algorithm and combines compound evolutionary algorithm, simplex algorithm and polynomial resampling. SP-UCI not only maintains the searchability of particle swarm in the whole parameter space but also efficiently realizes high-dimensional parameter optimization [41,50], and it is useful and effective global optimization. By using the SP-UCI algorithm, three sensitive parameters such as D0, ELUE and STT were optimized, and then the dynamic growth simulation was performed for the entire growth period of winter wheat using the SAFY model. The implementation steps of the data assimilation method (Figure 3b) are as follows:
(I).
AGB predicted model construction: the AGB prediction results based on the CBA-Wheat model are chosen as the state variable to estimate the yield in the assimilation system.
(II).
Run SAFY: SAFY model is run based on initialized model parameters and meteorological data.
(III).
Cost function calculation: The cost function is built on the basis of the relationship between the measured AGB and the model simulated AGB.
RMSE = 1 n i = 1 n ( Y i Y i ) 2
where, n, Yi′, Yi, and p are the number of samples, predicted AGB value, measured AGB value, and the number of independent variables, respectively.
(I).
Determine iteration termination conditions: When the objective function cannot be improved by 0.01% or the cost function is calculated more than 10,000 times to terminate the cycle.
(II).
Test the error between the model measured yield and the simulated yield.

3. Results

3.1. Validation of AGB Retrieved from CBA-Wheat Model and SAFY Model

Figure 4a shows that the trend of simulated AGB could reflect almost field-measured AGB in all nitrogen treatments. Four phases including the jointing phase, booting phase, flowering phase and filling phase were picked out. The distribution of AGB in the measured dataset, possible datasets, CBA-Wheat and parameters data assimilation dataset were illustrated in Figure 4b from jointing phase to filling phase, which demonstrates the broad coverage of the possible datasets. Moreover, the shapes of the four kinds of AGB datasets were similar. In this paper, the range and mean value of measured AGB were 0.46~19.50 t ha−1 and 6.96 t ha−1, respectively. The AGB ranges from the possible datasets of the jointing phase, booting phase, anthesis phase and filling phase were 0.10~7.72 t ha−1, 2.13~12.79 t ha−1, 3.20~19.11 t ha−1 and 5.24~26.36 t ha−1, respectively. The range of these values included the AGB values of other data sets in the corresponding period, indicating that the datasets input the transfer learning model include all possible situations of field production.

3.2. AGB Distribution for Different Datasets in Different Stages

The predicted AGB values were compared to 384 destructive measured AGB during four growth phases. The comparison between the measured AGB and the predicted AGB using CBA-Wheat (Figure 5) reveals a good correlation with R2 of 0.83 and RMSE of 1.91 t ha−1. CBA-Wheat and data assimilation have similar AGB estimation accuracy at different growth stages, which provides a solid foundation for the transfer learning method to replace the crop growth model for yield estimation.

3.3. Winter Wheat Yield Prediction

The simulated dataset with 60,000 samples covering the yield with 0.06~2.49 t ha−1 generated from the SAFY model was used to validate the effectiveness in yield inversions. The RMSE and loss value of the simulated dataset versus epochs are shown in Figure 6. When the training Epoch was more than 52, RMSE and loss value was less than 1.44 t ha−1 and 0.03, respectively. Meanwhile, the values converged to stable as the training epochs increased. Figure 7 demonstrates the high consistency between the measured and the predicted yield derived from the data assimilation and transfer learning method. The wheat yield prediction performance was as follows: transfer learning method (Figure 7a: R2 of 0.64, RMSE of 1.05 t ha−1) and data assimilation (Figure 7b: R2 of 0.64, RMSE of 1.01 t ha−1).
To compare the computational efficiency of the transfer learning method and data assimilation algorithm, three supplementary experiments are designed (Figure 8).The time consumption of transfer learning with complete simulation datasets was significantly lower than that of the other two yield estimation tests. When the number of pixels was about 16,000, the calculation efficiency of the data assimilation system was the same. When the number of pixels is more than 16,000, the transfer learning without complete simulation datasets had better computational efficiency compared with data assimilation system.

3.4. Farm-Land Verification of Transfer Learning Method

The SAFY and transfer learning method were verified on the farmland. The results of estimating yield with the data assimilation and transfer learning method are shown in Figure 9. At the farm scale, the two yield estimation models were still similar in performance for data assimilation (RMSE of 1.33 t ha−1 and R2 of 0.46) and for transfer learning (RMSE of 1.13 t ha−1 and R2 of 0.47). The farm-scale yield results (Figure 10) were consistent with the actual farm production. The above-mentioned results illustrated that, the transfer learning technique had good performance in calibration datasets and validation datasets. The transfer learning method, i.e., the DNN model, could exploit full merit of the prior knowledge extracted from the crop growth model to reach favorable yield predicting performance. Moreover, the generalization of the transfer learning method could also provide a potential insight for timely monitoring of wheat growth in the field.

4. Discussion

Previous studies have attempted to predict yield using optical remote sensing technology [5,20,21]. The present study proposed a novel methodology for crop yield estimation based on satellite imagery. The advantages and limitations of the transfer learning method would be discussed from three aspects: the advantage of applying CBA-Wheat to predict AGB, the comparison between the transfer learning method and data assimilation and potential extension and limitation.

4.1. Advantage of Applying CBA-Wheat to Predict AGB

The AGB level at different growth stages has an obvious impact on the increase or decrease of harvest yield because crop yield is part of AGB. Different relationships of AGB-EVI2 at various growing phases have previously been reported in Roth., [51] and Li et al., [32], hence this was a widespread problem. On the one hand, the difference in organic matter accumulation rate and vegetation spectral response with the growth period resulted in different AGB changes corresponding to the change of unit vegetation index in different growth periods [32,52]. On the other hand, the limited canopy information that can be detected by the optical vegetation index limits the application of AGB estimation in multiple growth periods [53,54]. To eliminate the phenomenon of appeal, plant height, phenological information, vertical distribution and multi-angle optical information have been used in the modified AGB unified model of the multi-growth period. CBA-Wheat is a model aiming at the fact that VI cannot be directly used for AGB estimation in different growth periods. Compared with the method of integrating other information with the vegetation index, CBA-Wheat is simple, efficient and accurate, and can meet the accuracy of transfer learning. In this paper, EVI2 is used to estimate biomass and yield. EVI2 has the advantages of sensitivity and anti-saturation at a high biomass level, and the band of EVI2 is included in the existing satellite data [38]. However, EVI2 is built based on optical satellites and is vulnerable to clouds, rain and other meteorological factors. When using the model in a larger area, it is necessary to consider the fusion of multi-source data (such as combining active remote sensing data, etc.).

4.2. Comparison between Transfer Learning Method and Data Assimilation

Previous studies demonstrated the capacity of using data assimilation to improve crop productivity estimation for a large area [3,40,45]. In this research, instead of many complex crop growth models such as the Decision Support System of Agrotechnology Transfer modeling system (DSSAT), the Agricultural Production Systems sIMulator (APSIM), we chose only one mechanically simpler model (SAFY model) [40,55,56]. The reason for using the SAFY model is its advantages of requiring fewer detailed inputs, and capable of producing similar outputs compared with the complex models [57,58]. With optimization of tuning cultivar parameters, the SAFY model was able to predict the yield interannually, with a high R2 of 0.64 to measured yield (Figure 7). The main reason for limiting its application can be attributed to low computational efficiency. Although the statistical method is highly efficient, it is not considered in this paper due to poor mechanism. The DNN network was trained on the AGB and yield generated from SAFY, then transferred to the field winter wheat yield prediction tasks. The proposed method has relatively simple architecture than data assimilation system. The yield prediction accuracies reached satisfactory results in the validation under various circumstances, including different locations and across years. The accuracy of the transfer model proposed in this study is similar to that of the data assimilation algorithm in yield prediction. Crop growth models are built according to the crop growth and development principles, this fact makes it possible to use SAFY simulated data as a base knowledge for winter wheat yield prediction [48]. However, the crop growth model is generality used to simulate crop growth and development under normal climate and management conditions, without considering extreme climatic conditions and pests [40,55,56]. Therefore, learning more crop growth models with different structures will help to predict yield and other parameters under different growth conditions. On the premise that the simulation data set has been trained well, the computational efficiency of transfer learning with simulated datasets was much higher than that of a data assimilation algorithm (Figure 8). It can be explained that the transfer learning method can express more efficient results with fewer parameters [40]. At the image pixel number is about 16,000, the time-consuming of the data assimilation algorithm and transfer learning without complete simulation datasets is the same. Building usable simulation data sets is the most time-consuming step in transfer learning without complete simulation datasets. The turning point of computing efficiency will be different in different computer equipment and running software. Figure 10 showed higher yield was observed in the area close to the edge of the fields. On the one hand, the pixel at the edge of the field is a mixed pixel of grass and wheat, which may overestimate the yield; on the other hand, the wheat growing at the edge of the field had abundant in light, temperature and other environmental conditions, so it has a high organic matter accumulation. The remote sensing images used in this paper are mainly optical data, whose data quality is affected by clouds and rain, and then affects the yield prediction of the farm [23]. In the future, it is necessary to use multi-source data for large-scale regional yield estimation.

4.3. Potential Extension and Limitation

The results provided reliable evidence that crop yield could be predicted using the newly proposed transfer learning method with good accuracy and expansibility similar to data assimilation. Some potential extensions can be considered although more development is still needed. Massive simulation data based on the crop growth model is an important foundation of the hybrid method and provide representative information for later yield estimation. This study has only focused on winter wheat yield, other output parameters prediction of other crop growth models for different crops such as LAI, plant nitrogen concentration and crop quality, can be applied to enhance the estimation performance. The proposed hybrid method could enhance the use of simulation data generated by crop growth models, which would reduce the difficulty of manual sampling and data acquisition in the application of yield production in large regions. If the constructed sample set is poorly representative, the prediction performance of the migration model is still poor [59]. It worth highlighting that high-quality dataset is necessary, in order to apply transfer learning technology in practical agriculture and reach high yield estimation accuracy. Therefore, it is strongly required to formulate more scientific and standardized field sampling standards and establish a representative database, which will have a far-reaching impact on the development of yield prediction. Yield estimation from a farm was used to verify the effectiveness of this method. Although the transfer learning method was verified in different N management conditions, the verification of other treatments or other crop species was incomplete. Future research remains indispensable to evaluate the feasibility of such a method, for the extension when upscaling to different ecological areas. Whether a set of simulation datasets can be shared for transfer learning under different production conditions also needs to be further studied.

5. Conclusions

In this paper, to effectively utilize the satellite data to predict wheat yield, a transfer learning strategy was proposed, which included CBA-Wheat algorithm and the transfer learning method with SAFY simulation dataset and satellite-based observations. The method was fully analyzed, validated against field observations and compared with the data assimilation method, with the main conclusions as follows: (1) the comparison between measured AGB and predicted AGB using CBA-Wheat all reveal a good correlation with R2 and RMSE of 0.83 and 1.91 t ha−1, respectively. (2) the performance of yield prediction was as follows: transfer learning method (R2 of 0.64, RMSE of 1.05 t ha−1) and data assimilation (R2 of 0.64, RMSE of 1.01 t ha−1). At the farm scale, the two yield estimation models are still similar in performance with RMSE of 1.33 t ha−1 for data assimilation and 1.13 t ha−1 for transfer and learning. (3) The time consumption of transfer learning with complete simulation data set is significantly lower than that of the other two yield estimation tests. While the number of pixels to be simulated is about 16,000, the computational efficiency of the data assimilation algorithm and transfer learning without complete simulation datasets.

Author Contributions

Conceptualization, Y.Z. (Yu Zhao), S.H. and Z.L.; methodology, Y.Z. (Yu Zhao), S.H. and Z.L.; software, Y.Z. (Yu Zhao); validation, Y.Z., S.H. and Y.M.; formal analysis, X.S. and G.Y.; investigation, G.Y. and Y.Z. (Yan Zhu); data curation, J.C., S.H. and Y.Z. (Yu Zhao); writing—original draft preparation, Y.Z. (Yu Zhao); writing—review and editing, Y.Z. (Yan Zhu); visualization, H.F.; funding acquisition, X.S., Z.L. and G.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the earmarked fund for China Agriculture Research System (CARS-03), the National Natural Science Foundation of China (Grant No. 42271396) and the National Key Research and Development Program (2019YFE0125300).

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like thank Weiguo Li, Hanzhong Liu and Hong Chang for acquiring data in the field experiments of this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lobell, D.B.; Asseng, S. Comparing estimates of climate change impacts from process-based and statistical crop models. Environ. Res. Lett. 2017, 12, 015001. [Google Scholar] [CrossRef]
  2. Zaks, D.P.M.; Kucharik, C.J. Data and monitoring needs for a more ecological agriculture. Environ. Res. Lett. 2011, 6, 014017. [Google Scholar] [CrossRef]
  3. Battude, M.; Al Bitar, A.; Morin, D.; Cros, J.; Huc, M.; Marais Sicre, C.; Dantec, V.; Demarez, V. Estimating maize biomass and yield over large areas using high spatial and temporal resolution sentinel-2 like remote sensing data. Remote Sens. Environ. 2016, 184, 668–681. [Google Scholar] [CrossRef]
  4. Huang, J.; Ma, H.; Sedano, F.; Lewis, P.; Liang, S.; Wu, Q.; Su, W.; Zhang, X.; Zhu, D. Evaluation of regional estimates of winter wheat yield by assimilating three remotely sensed reflectance datasets into the coupled WOFOST–PROSAIL model. Eur. J. Agron. 2019, 102, 1–13. [Google Scholar] [CrossRef]
  5. Jiang, H.; Hu, H.; Zhong, R.; Xu, J.; Xu, J.; Huang, J.; Wang, S.; Ying, Y.; Lin, T. A deep learning approach to conflating heterogeneous geospatial data for corn yield estimation: A case study of the USA corn belt at the county level. Glob. Chang. Biol. 2019, 26, 1–13. [Google Scholar] [CrossRef]
  6. Li, Z.; Taylor, J.; Yang, H.; Casa, R.; Jin, X.; Li, Z.; Song, X.; Yang, G. A hierarchical interannual wheat yield and grain protein prediction model using spectral vegetative indices and meteorological data. Field Crop. Res. 2020, 248, 107711. [Google Scholar] [CrossRef]
  7. Gorelick, N.; Hancher, M.; Dixon, M.; Ilyushchenko, S.; Thau, D.; Moore, R. Google earth engine: Planetary-scale geospatial analysis for everyone. Remote Sens. Environ. 2017, 202, 18–27. [Google Scholar] [CrossRef]
  8. Li, W.; Ciais, P.; Stehfest, E.; Vuuren, D.; Popp, A.; Arneth, A.; Fulvio, F.; Doelman, J.; Humpenoder, F.; Harper, A.B.; et al. Mapping the yields of lignocellulosic bioenergy crops from observations at the global scale. Earth Syst. Sci. Data 2020, 12, 789–804. [Google Scholar] [CrossRef] [Green Version]
  9. Huang, H.; Huang, J.; Li, X.; Zhuo, W.; Wu, Y.; Niu, Q.; Su, W.; Yuan, W. A dataset of winter wheat aboveground biomass in China during 2007–2015 based on data assimilation. Sci. Data 2022, 9, 200. [Google Scholar] [CrossRef]
  10. Descals, A.; Wich, S.; Meijaard, E.; Gaveau, A.L.A.; Pee dell, S.; Zantoi, Z. High-resolution global map of smallholder and industrial closed-canopy oil palm plantations. Earth Syst. Sci. Data 2021, 13, 1211–1231. [Google Scholar] [CrossRef]
  11. Tennakoon, S.; Murty, V.; Eiumnoh, A. Estimation of cropped area and grain yield of rice using remote sensing data. Int. Remote Sens. 1992, 13, 427–439. [Google Scholar] [CrossRef]
  12. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring Vegetation Systems in the Great Plains with ERTS; NASA Special Publication: Washington, DC, USA, 1974; Volum 1, pp. 48–62. [Google Scholar]
  13. Sims, D.A.; Rahman, A.F.; Cordova, V.D.; El-Masri, B.Z.; Baldocchi, D.D.; Bolstad, P.V.; Flanagan, L.B.; Goldstein, A.H.; Hollinger, D.Y.; Misson, L.; et al. A new model of gross primary productivity for North American ecosystems based solely on the enhanced vegetation index and land surface temperature from MODIS. Remote Sens. Environ. 2008, 112, 1633–1646. [Google Scholar] [CrossRef]
  14. Rondeaux, G.; Steven, M.; Baret, F. Optimization of soil-adjusted vegetation indices. Remote Sens. Environ. 1996, 55, 95–107. [Google Scholar] [CrossRef]
  15. Jin, X.; Yang, G.; Xu, X.; Yang, H.; Feng, H.; Li, Z.; Shen, J.; Zhao, C. Combined multi-temporal optical and radar parameters for estimating lai and biomass in winter wheat using HJ and RADARSAR-2 data. Remote Sens. 2015, 7, 13251–13272. [Google Scholar] [CrossRef] [Green Version]
  16. Zhou, X.; Zheng, H.; Xu, X.; He, J.; Ge, X.; Yao, X.; Cheng, T.; Zhu, Y.; Cao, W.; Tian, Y. Predicting grain yield in rice using multi-temporal vegetation indices from UAV-based multispectral and digital imagery. ISPRS J. Photogramm. Remote Sens. 2017, 130, 246–255. [Google Scholar] [CrossRef]
  17. Zhang, L.; Zhang, Z.; Luo, Y.; Cao, J.; Li, Z. Optimizing genotype-environment-management interactions for maize farmers to adapt to climate change in different agro-ecological zones across China. Sci. Total Environ. 2020, 728, 138614. [Google Scholar] [CrossRef]
  18. Lobell, D.B.; Asner, G.P. Climate and management contributions to recent trends in US agricultural yields. Science 2003, 299, 1032. [Google Scholar] [CrossRef]
  19. Bryan, B.; King, D.; Zhao, G. Influence of management and environment on Australian wheat: Information for sustainable intensification and closing yield gaps. Environ. Res. Lett. 2014, 9, 044005. [Google Scholar] [CrossRef]
  20. Li, L.; Wang, B.; Feng, P.; Wang, H.; He, Q.; Wang, Y.; Liu, D.L.; Li, Y.; He, J.; Feng, H.; et al. Crop yield forecasting and associated optimum lead time analysis based on multi-source environmental data across china. Agric. For. Meteorol. 2021, 308–309, 108558. [Google Scholar] [CrossRef]
  21. Cao, J.; Zhang, Z.; Tao, F.; Zhang, L.; Luo, Y.; Zhang, J.; Han, J.; Xie, J. Integrating multi-source data for rice yield prediction across china using machine learning and deep learning approaches. Agric. For. Meteorol. 2021, 297, 108275. [Google Scholar] [CrossRef]
  22. Moulin, S.; Bondeau, A.; Delecolle, R. Combining agricultural crop models and satellite observations: From field to regional scales. Int. J. Remote Sens. 1998, 19, 1021–1036. [Google Scholar] [CrossRef]
  23. Dorigo, W.A.; Zurita-Milla, R.; Wit, A.; Brazile, J.; Singh, R.; Schaepman, M.E. A review on reflective remote sensing and data assimilation techniques for enhanced agroecosystem modeling. Int. J. Appl. Earth Observ. Geoinform. 2007, 9, 165–193. [Google Scholar] [CrossRef]
  24. Jin, X.; Kumar, L.; Li, Z.; Feng, H.; Xu, X.; Yang, G.; Wang, J. A review of data assimilation of remote sensing and crop models. Eur. J. Agron. 2018, 92, 141–152. [Google Scholar] [CrossRef]
  25. Huang, J.; Tian, L.; Liang, S.; Ma, H.; Becker-Reshef, I.; Huang, Y.; Su, W.; Zhang, X.; Zhu, D.; Wu, W. Improving winter wheat yield estimation by assimilation of the leaf area index from Landsat TM and MODIS data into the WOFOST model. Agric. For. Meteorol. 2015, 204, 106–121. [Google Scholar] [CrossRef] [Green Version]
  26. Ines, A.V.M.; Das, N.N.; Hansen, J.W.; Njoku, E.G. Assimilation of remotely sensed soil moisture and vegetation with a crop simulation model for maize yield prediction. Remote Sens. Environ. 2013, 138, 149–164. [Google Scholar] [CrossRef] [Green Version]
  27. Hu, S.; Shi, L.; Zha, Y.; Williams, M.; Lin, L. Simultaneous state-parameter estimation supports the evaluation of data assimilation performance and measurement design for soil-water atmosphere-plant system. J. Hydrol. 2017, 555, 812–831. [Google Scholar] [CrossRef]
  28. Huang, J.; Gómez-Dans, J.L.; Huang, H.; Ma, H.; Wu, Q.; Lewis, P.E.; Liang, S.; Chen, Z.; Xue, J.H.; Wu, Y.; et al. Assimilation of remote sensing into crop growth models: Current status and perspectives. Agric. For. Meteorol. 2019, 276–277, 107609. [Google Scholar] [CrossRef]
  29. Li, P.; Zhang, X.; Wang, W.; Zheng, H.; Yao, X.; Tian, Y.; Zhu, Y.; Cao, W.; Chen, Q.; Cheng, T. Estimating aboveground and organ biomass of plant canopies across the entire season of rice growth with terrestrial laser scanning. Int. J. Appl. Earth Obs. Geoinf. 2020, 91, 102132. [Google Scholar] [CrossRef]
  30. Yang, H.; Yang, G.; Gaulton, R.; Zhao, C.; Li, Z.; Taylor, J.; Wicks, D.; Minchella, A.; Chen, E.; Yang, X. In-season biomass estimation of oilseed rape (Brassica napus L.) using fully polarimetric SAR imagery. Precis. Agric. 2019, 20, 630–648. [Google Scholar] [CrossRef] [Green Version]
  31. Beamish, A.; Raynolds, M.K.; Epstein, H.; Frost, G.V.; Macander, M.J.; Bergstedt, H.; Bartsch, A.; Kruse, S.; Miles, V.; Tanis, C.M.; et al. Recent trends and remaining challenges for optical remote sensing of Arctic tundra vegetation: A review and outlook. Remote Sens. Environ. 2020, 246, 111872. [Google Scholar] [CrossRef]
  32. Li, Z.; Zhao, Y.; Taylor, J.; Gaulton, R.; Jin, X.; Song, X.; Li, Z.; Meng, Y.; Chen, P.; Feng, H.; et al. Comparison and transferability of thermal, temporal and phenological-based in-season predictions of above-ground biomass in wheat crops from proximal crop reflectance data. Remote Sens. Environ. 2022, 273, 112967. [Google Scholar] [CrossRef]
  33. Guo, Y.; Fu, Y.; Hao, F.; Zhang, X.; Wu, W.; Jin, X.; Bryant, C.R.; Senthilnath, J. Integrated phenology and climate in rice yields prediction using machine learning methods. Ecol. Indic. 2020, 120, 106935. [Google Scholar] [CrossRef]
  34. Guo, L.; Wang, J.; Xiao, Z.; Zhou, H.; Song, J. Data-based mechanistic modelling and validation for leaf area index estimation using multi-angular remote-sensing observation time series. Int. J. Remote Sens. 2014, 35, 4655–4672. [Google Scholar] [CrossRef]
  35. Dhillon, M.S.; Dahms, T.; Kuebert-Flock, C.; Borg, E.; Conrad, C.; Ullmann, T. Modelling crop biomass from synthetic remote sensing time series: Example for the DEMMIN test site, Germany. Remote Sens. 2020, 12, 1819. [Google Scholar] [CrossRef]
  36. Zadoks, J.C.; Chang, T.T.; Konzak, C.F. A decimal code for the growth stages of cereals. Weed Res. 1974, 14, 415–421. [Google Scholar] [CrossRef]
  37. Bai, X.; Li, Z.; Li, W.; Zhao, Y.; Li, M.; Chen, H.; Wei, S.; Jiang, Y.; Yang, G.; Zhu, X. Comparison of Machine-Learning and CASA Models for Predicting Apple Fruit Yieldsfrom Time-Series Planet Imageries. Remote Sens. 2021, 13, 3073. [Google Scholar] [CrossRef]
  38. Jiang, Z.; Huete, A.R.; Didan, K.; Miura, T. Development of a two-band enhanced vegetation index without a blue band. Remote Sens. Environ. 2008, 112, 3833–3845. [Google Scholar] [CrossRef]
  39. Duchemin, B.; Maisongrande, P.; Boulet, G.; Benhadj, I. A simple algorithm for yield estimates: Evaluation for semi-arid irrigated winter wheat monitored with green leaf area index. Environ. Model. Softw. 2008, 23, 876–892. [Google Scholar] [CrossRef] [Green Version]
  40. Bengio, Y. Learning Deep Architectures for AI; Now Publishers Inc.: Delft, The Netherlands, 2009. [Google Scholar]
  41. Duan, Q.Y.; Gupta, V.K.; Sorooshian, S. Shuffled complex evolution approach for effective and efficient global minimization. J. Optim. Theory Appl. 1993, 76, 501–521. [Google Scholar] [CrossRef]
  42. Duan, Q.; Sorooshian, S.; Gupta, V.K. Optimal use of the SCE-UA global optimization method for calibrating watershed models. J. Hydrol. 1994, 158, 265–284. [Google Scholar] [CrossRef]
  43. Dong, T.; Liu, J.; Qian, B.; Zhao, T.; Jing, Q.; Geng, X.; Wang, J.; Huffman, T.; Shang, J. Estimating winter wheat biomass by assimilating leaf area index derived from fusion of Landsat-8 and MODIS data. Int. J. Appl. Earth. Obs. 2016, 49, 63–74. [Google Scholar] [CrossRef]
  44. Kang, Y.; Ozdogan, M. Field-level crop yield mapping with Landsat using a hierarchical data assimilation approach. Remote Sens. Environ. 2019, 228, 144–163. [Google Scholar] [CrossRef]
  45. Wang, E.; Martres, P.; Zhao, Z.; Ewert, F.; Maiorano, A.; Rötter, R.; Kimball, B.A.; Ottman, M.J.; Wall, G.W.; White, J.W.; et al. The uncertainty of crop yield projections is reduced by improved temperature response functions. Nat. Plants 2017, 3, 17102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Claverie, M.; Demarez, V.; Duchemin, B.; Hagolle, O.; Keravec, P.; Marciel, B.; Ceschia, E.; Dejoux, J.F.; Dedieu, G. Spatialization of crop leaf area index and biomass by combining a simple crop model SAFY and high spatial and temporal resolutions remote sensing data. In Proceedings of the 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, Africa, 12–17 July 2009. [Google Scholar] [CrossRef]
  47. Tzeng, E.; Hoffman, J.; Zhang, N.; Saenko, K.; Darrell, T. Deep Domain Confusion: Maximizing for Domain Invariance. arXiv 2014, arXiv:1412.3474. [Google Scholar]
  48. Yosinski, J.; Clune, J.; Bengio, Y.; Lipson, H. How Transferable Are Features in Deep Neural Networks? In Proceedings of the Advances in Neural Information Processing Systems 27 (NIPS 2014), Montreal, QC, Canada, 8–13 December 2014. [Google Scholar]
  49. Tan, C.; Sun, F.; Kong, T.; Zhang, W.; Yang, C.; Liu, C. A Survey on Deep Transfer Learning. In Proceedings of the 27th International Conference on Artificial Neural Networks, Rhodes, Greece, 4–7 October 2018. [Google Scholar]
  50. Chu, W.; Gao, X.; Sorooshian, S. A new evolutionary search strategy for global optimization of high-dimensional problems. Inf. Sci. 2011, 181, 4909–4927. [Google Scholar] [CrossRef] [Green Version]
  51. Roth, L.; Streit, B. Predicting cover crop biomass by lightweight UAS-based RGB and NIR photography: An applied photogrammetric approach. Precis. Agric. 2018, 19, 93–114. [Google Scholar] [CrossRef] [Green Version]
  52. Zhao, Y.; Meng, Y.; Feng, H.; Han, S.; Yang, G.; Li, Z. Should phenological information be applied to predict agronomic traits across growth stages of winter wheat? Crop J. 2022, 10, 1346–1352. [Google Scholar] [CrossRef]
  53. Yue, J.; Yang, G.; Li, C.; Li, Z.; Wang, Y.; Feng, H.; Xu, B. Estimation of winter wheat above-ground biomass using unmanned aerial vehicle-based snapshot hyperspectral sensor and crop height improved models. Remote Sens. 2017, 9, 708. [Google Scholar] [CrossRef] [Green Version]
  54. Duan, D.; Zhao, C.; Li, Z.; Yang, G. Estimating total leaf nitrogen concentration in winter wheat by canopy hyperspectral data and nitrogen vertical distribution. J. Integr. Agric. 2019, 18, 1562–1570. [Google Scholar] [CrossRef]
  55. Jones, J.W.; Hoogenboom, G.; Porter, C.H.; Boote, K.J.; Batchelor, W.D.; Hunt, L.A.; Wilkens, P.W.; Singh, U.; Gijsman, A.J.; Ritchie, J.T. The DSSAT cropping system model. Eur. J. Agron. 2003, 18, 235–265. [Google Scholar] [CrossRef]
  56. Keating, B.A.; Carberry, P.S.; Hammer, G.L.; Probert, M.E.; Robertson, M.J.; Holzworth, D.; Huth, N.I.; Hargreaves, J.N.G.; Meinke, H.; Hochman, Z.; et al. An overview of APSIM, a model designed for farming systems simulation. Eur. J. Agron. 2003, 18, 267–288. [Google Scholar] [CrossRef] [Green Version]
  57. Dzotsi, K.A.; Jones, J.W.; Adiku, S.G.K.; Naab, J.B.; Singh, U.; Porter, C.H.; Gijsman, A.J. Modeling soil and plant phosphorus within DSSAT. Ecol. Model. 2010, 221, 2839–2849. [Google Scholar] [CrossRef]
  58. Della Peruta, R.; Keller, A.; Schulin, R. Sensitivity analysis, calibration and validation of EPIC for modelling soil phosphorus dynamics in Swiss agro-ecosystems. Environ. Model. Softw. 2014, 62, 97–111. [Google Scholar] [CrossRef]
  59. Zhang, Y.; Hui, J.; Qin, Q.; Sun, Y.; Zhang, T.; Sun, H.; Li, M. Transfer-learning-based approach for leaf chlorophyll content estimation of winter wheat from hyperspectral data. Remote Sens. Environ. 2021, 267, 112724. [Google Scholar] [CrossRef]
Figure 1. Experimental farm location for field observations: (a) test site; (b) validation set area; (c) calibration set area. Orange area in Figure 1a represent the Changping district.
Figure 1. Experimental farm location for field observations: (a) test site; (b) validation set area; (c) calibration set area. Orange area in Figure 1a represent the Changping district.
Remotesensing 14 05474 g001
Figure 2. Daily temperature (a) and radiation (b) during 2017~2019 in field experiments.
Figure 2. Daily temperature (a) and radiation (b) during 2017~2019 in field experiments.
Remotesensing 14 05474 g002
Figure 3. The flowchart of transfer learning method (a) and data assimilation (b) for winter wheat. SAFY, CBA_Wheat, DNN, SPU-CI, AGB and EVI2 represent the simple algorithm for yield model, crop biomass algorithm for wheat, the deep neural network, the Shuffled complex evolution with PCA, aboveground dry biomass and the enhanced vegetation index 2, respectively.
Figure 3. The flowchart of transfer learning method (a) and data assimilation (b) for winter wheat. SAFY, CBA_Wheat, DNN, SPU-CI, AGB and EVI2 represent the simple algorithm for yield model, crop biomass algorithm for wheat, the deep neural network, the Shuffled complex evolution with PCA, aboveground dry biomass and the enhanced vegetation index 2, respectively.
Remotesensing 14 05474 g003
Figure 4. Winter AGB distribution in different stages. (a) the measured and the simulated AGB in four different stages. (b) the AGB distribution in the measured dataset, the possible dataset from the SAFY model (Section 2.3.3), CBA-wheat dataset and data assimilation of the four growth phases.
Figure 4. Winter AGB distribution in different stages. (a) the measured and the simulated AGB in four different stages. (b) the AGB distribution in the measured dataset, the possible dataset from the SAFY model (Section 2.3.3), CBA-wheat dataset and data assimilation of the four growth phases.
Remotesensing 14 05474 g004
Figure 5. Comparison of measured and predicted AGB with CBA-wheat model.
Figure 5. Comparison of measured and predicted AGB with CBA-wheat model.
Remotesensing 14 05474 g005
Figure 6. Performances of yield estimation model built using possible dataset and DNN: (a) RMSE, (b) loss value. RMSE, DNN and epochs represent root mean square error, the deep neural network and DNN iteration number, respectively.
Figure 6. Performances of yield estimation model built using possible dataset and DNN: (a) RMSE, (b) loss value. RMSE, DNN and epochs represent root mean square error, the deep neural network and DNN iteration number, respectively.
Remotesensing 14 05474 g006
Figure 7. Scatter plots of measured yield versus predicted yield (tons per hectare): transfer learning method (a) and data assimilation (b).
Figure 7. Scatter plots of measured yield versus predicted yield (tons per hectare): transfer learning method (a) and data assimilation (b).
Remotesensing 14 05474 g007
Figure 8. Comparison of computational efficiency between data assimilation and transfer learning algorithms. Computer configuration: Lenovo Legion R7000 2021 CPU R7-5800H @ 3.20 Hz and 16G random-access memory.
Figure 8. Comparison of computational efficiency between data assimilation and transfer learning algorithms. Computer configuration: Lenovo Legion R7000 2021 CPU R7-5800H @ 3.20 Hz and 16G random-access memory.
Remotesensing 14 05474 g008
Figure 9. Predicted yield (tons per hectare) derived from transfer learning models (a) and data assimilation (b).
Figure 9. Predicted yield (tons per hectare) derived from transfer learning models (a) and data assimilation (b).
Remotesensing 14 05474 g009
Figure 10. Winter wheat yield maps over experimental seasons: (a) 2018 and (b) 2019 at Xiaotangshan derived using Planet satellite imagery with the transfer learning models.
Figure 10. Winter wheat yield maps over experimental seasons: (a) 2018 and (b) 2019 at Xiaotangshan derived using Planet satellite imagery with the transfer learning models.
Remotesensing 14 05474 g010
Table 1. Planet satellite remote sensing data parameters.
Table 1. Planet satellite remote sensing data parameters.
BandsWavelength (nm)Spatial Resolution (m)
Blue455–5153
Green500–5903
Red560–6703
NIR780–8603
Table 2. Parameter descriptions in the SAFY model.
Table 2. Parameter descriptions in the SAFY model.
ParameterAbbreviationUnitValueReferences
FixedClimatic efficiency ε c -0.48Duchemin et al. [39]
Light-interception coefficient ε 1 -0.5Duchemin et al. [39]
The optimal temperatureTopt°C21Wang et al. [45]
Minimum temperaturTmin°C0Wang et al. [45]
Maximum temperaturTmax°C37Wang et al. [45]
Specific leaf areaSLAm2/g0.022Claverie et al. [46]
Initial aboveground biomassAGB0g/m24.5Calibrated
Leaf Partitioning Coefficient aPla-0.16Calibrated
Leaf Partitioning Coefficient bPlb-1.4Calibrated
Senesce rateRs°C/d10.8Calibrated
CalibratedDay of emergenceD0d0–15This study
Sum of temperatureSTT°C1200–1600This study
Effective Light Use EfficiencyELUEg/MJ1.3–2.5This study
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhao, Y.; Han, S.; Meng, Y.; Feng, H.; Li, Z.; Chen, J.; Song, X.; Zhu, Y.; Yang, G. Transfer-Learning-Based Approach for Yield Prediction of Winter Wheat from Planet Data and SAFY Model. Remote Sens. 2022, 14, 5474. https://doi.org/10.3390/rs14215474

AMA Style

Zhao Y, Han S, Meng Y, Feng H, Li Z, Chen J, Song X, Zhu Y, Yang G. Transfer-Learning-Based Approach for Yield Prediction of Winter Wheat from Planet Data and SAFY Model. Remote Sensing. 2022; 14(21):5474. https://doi.org/10.3390/rs14215474

Chicago/Turabian Style

Zhao, Yu, Shaoyu Han, Yang Meng, Haikuan Feng, Zhenhai Li, Jingli Chen, Xiaoyu Song, Yan Zhu, and Guijun Yang. 2022. "Transfer-Learning-Based Approach for Yield Prediction of Winter Wheat from Planet Data and SAFY Model" Remote Sensing 14, no. 21: 5474. https://doi.org/10.3390/rs14215474

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop