Next Article in Journal
Linear Displacement Calibration System Integrated with a Novel Auto-Alignment Module for Optical Axes
Next Article in Special Issue
SimTalk: Simulation of IoT Applications
Previous Article in Journal
Mechanical and Electronic Video Stabilization Strategy of Mortars with Trajectory Correction Fuze Based on Infrared Image Sensor
Previous Article in Special Issue
Best Water Vapor Information Layer of Himawari-8-Based Water Vapor Bands over East Asia
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 20
Issue 9
10.3390/s20092460
sensors-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Multi-Species Leaf Area Index Based on Chinese GF-1 Satellite Data Using Look-Up Table and Gaussian Process Regression Methods

by
Yangyang Zhang
1,2,3,
Jian Yang
1,2,*,
Xiuguo Liu
1,
Lin Du
1,
Shuo Shi
4,
Jia Sun
1 and
Biwu Chen
4
1
School of Geography and Information Engineering, China University of Geosciences, Wuhan 430074, China
2
Artificial Intelligence School, Wuchang University of Technology, Wuhan 430223, China
3
College of Resources and Environmental Sciences, Gansu Agricultural University, Lanzhou 730070, China
4
State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan 430079, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(9), 2460; https://doi.org/10.3390/s20092460
Submission received: 30 March 2020 / Revised: 22 April 2020 / Accepted: 24 April 2020 / Published: 26 April 2020
(This article belongs to the Special Issue Simulation Studies on Remote Sensing Scenarios)
Browse Figures
Versions Notes

Abstract

:
Leaf area index (LAI) is an important biophysical parameter, which can be effectively applied in the estimation of vegetation growth status. At present, amounts of studies just focused on the LAI estimation of a single plant type, while plant types are usually mixed rather than single distribution. In this study, the suitability of GF-1 data for multi-species LAI estimation was evaluated by using Gaussian process regression (GPR), and a look-up table (LUT) combined with a PROSAIL radiative transfer model. Then, the performance of the LUT and GPR for multi-species LAI estimation was analyzed in term of 15 different band combinations and 10 published vegetation indices (VIs). Lastly, the effect of the different band combinations and published VIs on the accuracy of LAI estimation was discussed. The results indicated that GF-1 data exhibited a good potential for multi-species LAI retrieval. Then, GPR exhibited better performance than that of LUT for multi-species LAI estimation. What is more, modified soil adjusted vegetation index (MSAVI) was selected based on the GPR algorithm for multi-species LAI estimation with a lower root mean squared error (RMSE = 0.6448 m2/m2) compared to other band combinations and VIs. Then, this study can provide guidance for multi-species LAI estimation.

1. Introduction

The leaf area index (LAI) is usually defined as one half of the total green leaf area per unit of horizontal ground surface area [1]. As the key vegetation biophysical parameter, LAI can be efficiently applied in monitoring vegetation photosynthesis, respiration, transpiration, soil respiration, and the energy exchange of the canopy-atmosphere [2,3]. Hence, LAI has been generally applied in many fields, such as precision agriculture, vegetation health monitoring, and crop yield estimation [4,5,6,7]. Therefore, it is significant to accurately invert vegetation LAI in real time. At present, the field measurement methods of LAI mainly contain the direct destructive sampling method, the indirect digital hemispherical photography method, plant canopy analyzers, and smartphones application named PocketLAI [8]. However, the classical methods are difficult to obtain LAI for large areas. In addition, crop LAI is also difficult to obtain by field measurements due to the high spatial and temporal heterogeneity.
Remote sensing techniques provide an effective method for quantitative assessment of LAI due to the advantage of fast, non-destructive and large scales. Many studies have successfully retrieved LAI through using remote sensing data [9,10,11,12,13,14,15]. The LAI inversion methods by remote sensing data are generally divided into two main aspects: statistical and physical approaches [16,17,18]. The statistical approaches are usually in term of the linear or nonlinear relationship between ground-based LAI measurements and spectral vegetation indices [19,20], which are simple and fast. However, a certain physical mechanism and universality is lacking [21,22]. The physical inversion strategies are built based on the radiative transfer models (RTM) [23,24]. There are three different inversion strategies for LAI estimation: iterative optimization, look-up table (LUT), and hybrid inversion [25,26]. The iterative optimization is a classical inversion approach, but it needs high computational cost and exhibits poor convergence. Additionally, many studies have shown that the LUT and hybrid inversion exhibited better performance for LAI retrieval [27,28].
The LUT relies on a large multiple-solutions database (rather than the single best solution) for the inversion of vegetation parameters. Measurements or model errors have little influence on this inversion method [29], which improved the robustness of LAI inversion [30], the LUT method has been widely used for LAI inversion in many studies due to the advantage of stability and robustness [31,32,33]. Hybrid inversion, which was combined the RTM with machine learning regression algorithms (MLRAs), uses RTM to generate simulated datasets as the training set of MLRAs. The hybrid inversion has the advantages of the statistical approaches’ simplicity and the physical method’s universality, which can accurately and rapidly invert vegetation LAI [34,35]. All kind of MLRAs are applied for the inversion of RTM, for instance, artificial neural network (ANN), support vector regression (SVR), random forest (RF), and Gaussian process regression (GPR) [36,37,38,39]. Several studies have shown that GPR exhibits more performance for LAI estimation among these MLRAs [36,40,41].
Using a high-resolution satellite is an effective way to extract LAI information, which can better meet the requirements for monitoring geographical conditions. Thus, the China High Resolution Earth Observation System (CHEOS) affords high spatial–temporal and near-real-time (NRT) observation information for vegetation estimation. Chinese satellites of High-Resolution Satellite 1 (GF-1) with a swath width of 800 km, coverage repetitive cycle of 4 days, spatial resolution of 16 m, the first satellite of CHEOS, was launched on 26 April 2013. In this study, GF-1 data was used to invert vegetation LAI and generate the LAI map.
With the development of remote sensing techniques, many studies have been conducted on LAI inversion based on remote sensing data. However, most of these studies may be only focused on the single species LAI estimation, such as wheat [32,33,42,43], maize [44,45], and sugar beet [26,46]. Relatively few research has been carried out on LAI inversion of multi-species. However, plant types are mostly mixed rather than single distribution. It is worth exploring whether the LAI of multi-species can be rapidly and efficiently retrieved based on remote sensing data.
Therefore, the objectives of this study are to (1) evaluate the suitability of GF-1 data for multi-species LAI estimation by using the LUT and GPR algorithms combined with the PROSAIL RTM; (2) explore the performance of the LUT and GPR for estimating multi-species LAI in term of 15 different band combinations and ten published vegetation indices (VIs); and (3) analyze the effect of the different band combinations and VIs on the accuracy of LAI estimation.

2. Study Area and Data

2.1. Study Area

Heihe River Basin of China is located in the middle of the Hexi Corridor (96°42′–102°00′E, 37°41′–42°42′N), which is the second largest inland river basin in the northwest China. Heihe River Basin is located in the center of Eurasia, which is far away from the ocean and surrounded by mountains. The climate of the basin is mainly controlled by westerly circulation in mid-high latitude and is influenced by cold air mass in the polar regions. Thus, it has a typical continental climate, the precipitation is rare and concentrated, the sunshine is sufficient and the temperature difference between day and night is large. In this study, the lower reaches of the Heihe River Basin are selected as the study area (Figure 1). The study area includes a large of the Gobi Desert, which contains desert riparian forests, shrub, and meadow vegetation [47].

2.2. Field Measurement Data

As the study area is located in the northwest of China, the average temperature and precipitation in July will reach the maximum value, and vegetation grows vigorously, which is conducive to the measurement and evaluation of LAI. Then, field LAI observation started from 22 July and ended on 1 August in 2014. Considering the characteristics of downstream vegetation distribution, uniform vegetation with large area distribution (over 100 m × 100 m) was selected as the observation sample plot. A total of 29 observation sample plots were selected, and the plant species included Hami melon, Tamarix chinensis, reed, weed, flower wood, bitter bean, and so on (Figure 1). LAI field observations were non-destructively measured by using the LAI-2200 plant canopy analyzer (LI-COR, Inc., Lincoln, NE, USA). The average LAI was obtained based on one above-canopy measurement and four below-canopy measurements of the incoming radiation in each sample plot, every sample plot was repeatedly measured at least 2 times, and the measurement time of all samples was from 8:00 to 12:00 [48]. Due to the experimental conditions limited, the data of the field measurements were downloaded from the Cold and Arid Regions Science Data Center. Then, every plot was matched with GF-1 pixel by using the ENVI 5.3 software according to the longitude and latitude coordinates.

2.3. Remote Sensing Data

GF-1 Wild Field Camera (WFV) data was used to estimate LAI of the study area. GF-1 was launched by China on 26 April 2013, and the GF-1 WFV camera can obtain 16 m multi-spectral color images (Blue, Green, Red, and Near-infrared bands). GF-1 satellite is equipped with four WFV camera sensors, namely WFV1, WFV2, WFV3, and WFV4, which simultaneously image and have different absolute radiometric calibration coefficients in orbit. The imaging swath width of each WFV sensor is about 200 km, and the imaging swath width of the four WFV sensors can reach about 800 km. The key parameters of the WFV sensor were listed in Table 1. Due to the existence of clouds or haze, a high-quality GF-1 WFV image data on 29 July in 2014, which was covered in the study area, was downloaded.
To convert the digital number (DN) of the raw image to land surface reflectance, the ENVI 5.3 software was used to preprocess the GF-1 image data in this study, including radiometric calibration, atmospheric correction, and geometric correction.

3. Methods

3.1. Inversion Schemes

Figure 2 shows the inversion flow chart of this study. Firstly, LAI was evaluated based on the reflectance of GF-1 WFV and VIs. In order to fully assess the performance of the spectral information for LAI estimation, 15 combinations of four bands were used (Table 2). In addition, LAI and VIs had a good correlation according to previous studies, so 10 types of VIs were selected, which were listed in Table 3. These VIs, which included simple difference indices, a normalized difference ratio, triangular vegetation indices, and modified versions, were commonly used to estimate LAI [19,31,33]. The relative wavelengths of VIs were located in the visible and near-infrared bands and in the GF-1 band ranges. Then, LUT and GPR methods were adopted to realize LAI inversion based on remote sensing spectral bands and VIs by using PROSAIL model. Finally, the optimal strategies were achieved through the verification of ground measurements data.

3.2. LUT Inversion

The LUT inversion is a conceptually simple approach to retrieve LAI based on a large multi-solution database. Thus, the database determines the implementation of LUT. Generally, the simulated dataset has been generated by PROSAIL RTM, which is a combination of the SAIL (scattering by arbitrarily inclined leaves) [60] canopy reflectance model and PROSPECT leaf optical properties model [26]. In this study, PROSAIL5B, a combination of PROSPECT5 and 4SAIL, was used to simulate canopy reflectance with many input canopy biophysical parameters through the output leaf reflectance and transmittance of PROSPECT with some input leaf biochemical parameters to SAIL [61]. To build LUT, the PROSAIL model was used to simulate the broadband canopy reflectance of GF-1 WFV through randomly generated input parameters. The key input parameters of PROSAIL are summarized by many literatures, which was listed in Table 4 [33,62,63,64,65]. Due to the lack of understanding of the actual growth status of vegetation in the study area, the simulation dataset should contain enough vegetation types and states to improve the robust of LUT and GPR for LAI estimation. In addition, the brown pigment of the PROSPECT model was also set. To get appropriate inversion results, Weiss et al. found that the size of LUT needs 100,000 reflectance/VIs, which can well balance the contradiction between saving computing resources and improving inversion accuracy [27].
To select the solution of the inverse problem, the LUT solution can be calculated by averaging similar parameter combinations with the smallest differences between the reflectance/VIs and the simulated reflectance/VIs [61,66,67]. The relative root mean square error (RRMSE) was used as the cost function to assess the performance of model for LAI estimation in this study, and was defined as follows:
R R M S E R = 1 m ∑ i = 1 m ( R R S , i − R S i m u l a t e d , i R R S , i ) 2
R R M S E V I = ( V I R S , i − V I S i m u l a t e d , i V I R S , i ) 2
where RRS is the measured band reflectance; RSimulated is the simulated reflectance; VIRS is the measured VI; VISimulated is the simulated VI, which was calculated from the simulated reflectance; and m is the number of spectral bands.
A set of reflectance/VIs variables corresponding to the minimization of the RRMSE value in the LUT are taken as the solution. However, the solution may not be the only and optimal result due to measurement errors and model deficiencies. The solution selected a threshold 10% of the ordered minimum RRMSE and these parameter combinations were averaged with the LUT size of 100,000, which can solve the ill-posed inversion problem. It was feasible for the accuracy of LAI estimation with the 10% threshold, and was consistent with the previous studies [61,67].

3.3. GPR Inversion

GPR is a nonparametric model for regression analysis of data using the Gaussian process prior. The model hypothesis of GPR includes two parts: noise (regression residual) and Gaussian process prior, and the essence of its solution is Bayesian inference. Without limiting the form of kernel function, GPR is theoretically a universal approximation of any continuous function on a compact space. GPR uses the Gaussian process as a priority, which assumes that the training samples are the sampling of Gaussian process, which resulted in estimation results are closely related to the kernel function. The practical significance of the kernel function in GPR is a covariance function, which describes the correlation among training samples. In this study, a squared exponential kernel function was applied in the GPR model, which was a widely used standard covariance function. Moreover, the squared exponential kernel function can capture sample similarity well in most problems and just one hyperparameter needed to be tuned. The detailed introduction of GPR and kernel functions can be referred to in [36,68,69,70].
In this paper, GPR inversion uses the PROSAIL model to build a simulated dataset that including simulated spectral reflectance/VIs and corresponding simulated LAI as the training dataset, to obtain the regression relationship between spectral reflectance/VIs and LAI. This inversion method has the advantages of simple empirical methods and strong university of physical model, and can accurately and rapidly retrieve LAI. The number of GPR training datasets would be determined in the next part.

3.4. Statistical Evaluation

To evaluate the performance of LUT and GPR for LAI estimation, the coefficient of determination (R2) and root-mean-square error (RMSE) were used. Sets y i and y i ’ denote the measured and predicted values respectively, y ¯ refers to the average of y i , and n is the number of measurements:
R 2 = 1 − ∑ i = 1 n ( y i − y i ′ ) 2 ∑ i = 1 n ( y i − y ¯ ) 2
R M S E = ∑ i = 1 n ( y i − y i ′ ) 2 n

4. Results

4.1. Determination of the Number of GPR Training Datasets

The synthetic dataset generated by the PROSAIL model was used as the training dataset of GPR model, which solved the demand of the model for the measured dataset. The training dataset directly influences the inversion accuracy of the GPR model. However, the training dataset is not positively correlated with the inversion accuracy of the model. This possible interpretation is that the accumulation of the random Gaussian noise and the occurred of the overfitting of the training model with the increasing of the training dataset number may be influenced by the accuracy of the LAI estimation. In this study, a different number of training datasets was used to analyze the accuracy of LAI inversion based on GPR.
As it can be seen from Figure 3, when the number of training datasets was 3000, the inversion accuracy of GPR was optimal. García-Haro, F. J. also introduced that when the number of training samples by PROSAIL model generation was around 3000, the RMSE values of LAI estimation becomes stable [65]. Therefore, 3000 simulated datasets were selected as training datasets of GPR model in this paper.

4.2. Different Band Combinations on LAI Inversion

The performance of LUT and GPR for LAI inversion with different band combinations (Table 2) was evaluated and was shown in Figure 4. Theoretically, the best inversion strategy is to have the highest R2 and the lowest RMSE, but in most cases, R2 and RMSE are not always consistent. For the LAI estimation based on the LUT-Band strategies with the band reflectance of GF-1 data (Figure 4a), the B2B4 (R2 = 0.4258, RMSE = 0.6700 m2/m2) exhibits the optimal performance for LAI inversion by the comparison of R2 and RMSE values, which was selected to estimate LAI in the following study. In addition, the B3B4 (red, Nir; R2 = 0.4260, RMSE = 0.7786 m2/m2) and B2B3B4 (green, red, Nir; R2 = 0.4290, RMSE = 0.7848 m2/m2) also displayed the potential for predicting the LAI based on GF-1 data. The remaining LUT-Band had a relatively lower RMSE values (RMSE < 1.2 m2/m2), but the R2 values of LAI inversion were lower than 0.3. However, the R2 and RMSE of the B4 (Nir; R2 = 0.4026, RMSE = 1.3767 m2/m2) performed very inconsistently.
The GPR-Band strategies based on GF-1 data (Figure 4b) showed an inconsistency in R2 and RMSE. The B2 (green) and B3 (red) strategies had the lowest RMSE values (0.8470 m2/m2 and 0.8963 m2/m2) but also had the worst R2 values (0.0096 and 0.0102). The B2B4 (Green, Nir) and B3B4 (Red, Nir) of GPR performed consistently with a lower RMSE values (0.9856 m2/m2 and 0.9975 m2/m2) and a higher R2 values (0.4042 and 0.3906). The B4 (Nir) of GPR on LAI estimation appeared to be similar to the LUT-B4, with the highest R2 value of 0.4212 and a higher RMSE value of 1.3792 m2/m2. All the other GPR-Band strategies on LAI inversion performed poor accuracy with all the R2 lower than 0.3.
These results indicated that it is important to select the appropriate band combination strategy to improve the accuracy of LAI estimation by using LUT and GPR methods. Thus, the optimal inversion strategy should be selected by comprehensive consideration of R2 and RMSE. For LUT and GPR, the B2B4 (green, Nir) strategy exhibited the optimal performance for LAI estimation. However, the GPR took only one-third of the time of LUT. Thus, the B2B4 was selected as the spectral characteristics parameter to analyze LAI by using the GPR model.

4.3. Different Published VIs on LAI Estimation

The performance of LUT and GPR for LAI inversion with different VIs (Table 3) was evaluated and was shown in Figure 5. For LUT and GPR, ten published VIs exhibited the basically similar performance for LAI estimation. Modified soil adjusted vegetation index (MSAVI) strategy got the best accuracy of LAI inversion based on the GF-1 data in both LUT and GPR with the low RMSE values (0.6338 m2/m2 and 0.6448 m2/m2) and a high R2 (0.4303 and 0.4257). In particular, the green normalized difference vegetation index (GNDVI), green ratio vegetation index (GRVI), and normalized difference vegetation index (NDVI) strategies all achieved a poor accuracy with a lower R2 and a higher RMSE for both LUT and GPR. Apart from the above VIs, the remaining VIs of LUT and GPR performed moderately.
Figure 5 shows that different VIs could achieve different inversion results for LUT and GPR methods. For the comparison of the R2 and RMSE values of LAI estimation based on different published VIs, MSAVI was considered as the optimal strategy to estimate multi-species LAI based on the GF-1 satellite data and selected for LAI estimation in the following study.

4.4. Comparison of LUT and GPR with Different Inversion Strategies on LAI Estimation

In this study, the suitability of GF-1 data for multi-species LAI estimation was discussed by using LUT and GPR algorithms based on the PROSAIL RTM. The accuracy of VIs-based inversion was higher than that of band-based inversion by using LUT. MSAVI exhibited the optimal performance for LAI estimation with the low RMSE value of 0.6338 m2/m2 and the relatively high R2 value of 0.4304 (Figure 4a and Figure 5a). Similar to the LUT inversion results, the VIs strategies also performed better performance for the LAI estimation than that of the band-based by using GPR. Additionally, MSAVI still achieved the optimal performance based on GPR with the low RMSE value of 0.6448 m2/m2 and the relatively high R2 value of 0.4257 (Figure 4b and Figure 5b). The reason of why the VIs inversion results was superior to that of the band-based may be that the VIs not only included sensitive bands of LAI, and also strengthened the difference among the reflectivity of bands, and weakened the influences of some factors (soil background). In addition, it is significant to choose an appropriate inversion strategy than an inversion method in order to improve the accuracy of LAI inversion. According to the accuracy of LAI inversion, the MSAVI is the optimal VIs by using LUT and GPR. However, GPR takes only one-third of the time of LUT regardless of the inversion strategy. Thus, considering the factors of precision and efficiency, MSAVI based on GPR was select as the optimal inversion method to estimate multi-species LAI from GF-1 satellite data in this study.

4.5. Mapping of Multi-Species LAI

Based on the above comparative analysis results, MSAVI based on GPR algorithm was selected for multi-species LAI retrieval by the GF-1 WFV data.
Figure 6 shows the LAI map generated from MSAVI strategy of GPR based on GF-1 data on 29 July, 2014. The LAI inversion values of the study area ranged from 0 to 4. Since the study area included a large area of the Gobi Desert, LAI values of the whole region were all low. The low LAI values were distributed in the west of the study area, and the high LAI values mostly appeared in the middle region. In general, the LAI retrieval map was basically consistent with the actual situations, and there was good consistency between the inversion and measured LAI values with the RMSE value of 0.6448 m2/m2.

5. Discussion

Multi-species LAI (e.g., Hami melon, Tamarix chinensis, reed, weed, flower wood, bitter bean, and so on) obtained from the GF-1 WFV sensor data was evaluated by comparing the different bands and VIs strategies based on the LUT and GPR methods. The optimal inversion strategy and method for multi-species LAI retrieval were determined. Additionally, the LAI map with the RMSE of 0.6448 m2/m2 was finally acquired through the optimal inversion strategy. Compared with the previous studies [71,72,73], the optimal inversion strategy exhibited better performance for multi-species LAI estimation. Thus, the GF-1 WFV data has good potential being used for multi-species LAI retrieval by choosing the appropriate inversion strategy. However, there are still some issues that need to be solved.
The inversion accuracy of LAI by using LUT is generally similar to that of GPR, which is in agreement with the previous studies [27,28]. However, the calculative process of LUT is slower and more complex than that of GPR, and some initial values of parameters also need to be set. Additionally, GPR is a hybrid inversion approach based on the RTM, which is different from the empirical method and does not need lots of measured data. The hybrid approach only requires measured data to verify the inversion accuracy. Compared with other inversion methods based on RTM, the hybrid inversion model has advantages in simple, fast, and accurate calculation in estimating vegetation parameters [35,74].
Different inversion strategies have a great important influence on the accuracy of LAI estimation. Fifteen band combinations and 10 VIs strategies of LUT and GPR displayed different performance for LAI estimation. In particular, B2B4 and MSAVI are the optimal methods on LAI inversion by GF-1 data. Both MSAVI and B2B4 strategies contain a near-infrared wavelength, which indicates that the near-infrared band of the GF-1 WFV sensor is sensitive to LAI estimation with multi-species on the underlying surface. It provides a reference for LAI inversion under similar underlying surface using the GF-1 WFV sensor data in the future study.
Due to the variety of plant types that was included in the study area, the improvement of LAI inversion accuracy is limited to a certain extent. The possible interpretation is that LAI of different plant types presents different special characteristics [75,76], as well as the application of inversion strategies. At present, many studies’ main focus on estimating LAI for a single plant type [3,14,32,33,43,53,67], while relatively few studies also discussed LAI estimation from multiple plant types [71,72,73]. These studies estimated global LAI by using different inversion strategies based on Sentinel-2 satellite data. The performance of Sentinel-2-derived for LAI estimation was quantified using global ground observation with consistent measurement criteria [73], and the RMSE value of Sentinel-2 LAI was 1.09 m2/m2 through comparison to ground LAI. The Simplified Level 2 Product Processor (SL2P) was used to estimate multi-crop LAI by Sentinel-2/MSI and Landsat-8/OLI data and was validated by using an agricultural area [71]. Additionally, the RMSE for MSI/PL2P and OLI/SL2P were 0.98 m2/m2 and 1.63 m2/m2 for LAI estimation, respectively. A new simple Sentinel-2 LAI index (SeLI) was proposed and used for multi-crop green LAI estimation [72]. SeLI performed well and had the potential to be used for LAI estimation of a multi-crop dataset (RMSE = 0.69 m2/m2). Compared with previous studies, it is feasible to use GF-1 data for multi-species LAI estimation. To improve the LAI inversion accuracy of multiple plant types, an efficient use of prior knowledge and accurate classification of different plant types and other effective methods will be studied in future work. In addition, LAI estimation of different species is also need to be further discussed by using classification and clustering in the following work.

6. Conclusions

In this study, the performance of GF-1 data for multi-species LAI estimation by using LUT and GPR algorithms combined with the PROSAIL RTM was analyzed. The effect of the 15 different bands combinations and 10 published VIs on the accuracy of LAI estimation was also assessed. The results demonstrated that the accuracy of VIs-based inversion was higher than that of band-based inversion for LAI estimation by using LUT and GPR. Moreover, the MSAVI strategy got the optimal performance of LAI inversion based on GF-1 data in both LUT and GPR with a high R2 (0.4303 and 0.4257) and low RMSE values (0.6338 m2/m2 and 0.6448 m2/m2). Nevertheless, GNDVI based on LUT and GPR had the worst accuracy of LAI estimation with the RMSE of 0.7790 m2/m2 and 0.8070 m2/m2, and the R2 of 0.2788 and 0.2863, respectively. These results indicated that different inversion methods had a great influence on LAI estimation. Selecting an appropriate inversion model can effectively improve the accuracy of the LAI estimation. If just considering the R2 values, the performance of the LUT (R2 = 0.4303) may be similar to that of GPR (R2 = 0.4257) for LAI estimation based on MSAVI. However, the running time of GPR was only one-third of that of LUT. Thus, GPR combined with MSAVI was optimal and was selected for multi-species LAI estimation in term of the GF-1 data. This study suggested that GF-1 satellite data performed well and had the potential to be used for multi-species LAI estimation. In addition, it is of great significance and value to research multi-species LAI inversion by using remote sensing satellite data.

Author Contributions

Conceptualization, Y.Z. and J.Y.; methodology, Y.Z., X.L. and J.Y.; validation, Y.Z. and J.Y.; formal analysis, Y.Z. and J.Y.; writing—original draft preparation, Y.Z. and J.Y.; writing—review and editing, X.L., L.D., S.S., J.S. and B.C.; supervision, X.L. and J.Y.; funding acquisition, J.Y. All authors have commented and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key R&D Program of China (2018YFB0504500); National Natural Science Foundation of China (NSFC; 41801268); Natural Science Foundation of Hubei Province (2018CFB272), and Fundamental Research Funds for the Central Universities, China University of Geosciences (Wuhan; CUG170661, CUGQY1947).

Acknowledgments

The data set is provided by “Cold and Arid Region Science Data Center at Lanzhou” (http://westdc.westgis.ac.cn).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, J.M.; Black, T.A. Defining Leaf-Area Index for Non-Flat Leaves. Plant Cell Environ. 1992, 15, 421–429. [Google Scholar] [CrossRef]
  2. Soudani, K.; François, C.; Maire, G.L.; Dantec, V.L.; Dufrêne, E. Comparative analysis of IKONOS, SPOT, and ETM+ data for leaf area index estimation in temperate coniferous and deciduous forest stands. Remote Sens. Environ. 2006, 102, 161–175. [Google Scholar] [CrossRef] [Green Version]
  3. Campos-Taberner, M.; Garcia-Haro, F.J.; Camps-Valls, G.; Grau-Muedra, G.; Nutini, F.; Crema, A.; Boschetti, M. Multitemporal and multiresolution leaf area index retrieval for operational local rice crop monitoring. Remote Sens. Environ. 2016, 187, 102–118. [Google Scholar] [CrossRef]
  4. Jay, S.; Maupas, F.; Bendoula, R.; Gorretta, N. Retrieving LAI, chlorophyll and nitrogen contents in sugar beet crops from multi-angular optical remote sensing: Comparison of vegetation indices and PROSAIL inversion for field phenotyping. Field Crops Res. 2017, 210, 33–46. [Google Scholar] [CrossRef] [Green Version]
  5. 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]
  6. Mora, M.; Avila, F.; Carrasco-Benavides, M.; Maldonado, G.; Olguín-Cáceres, J.; Fuentes, F. Automated computation of leaf area index from fruit trees using improved image processing algorithms applied to canopy cover digital photograpies. Comput. Electron. Agric. 2016, 123, 195–202. [Google Scholar] [CrossRef] [Green Version]
  7. Cotter, M.; Asch, F.; Hilger, T.; Rajaona, A.; Schappert, A.; Stuerz, S.; Yang, X. Measuring leaf area index in rubber plantations—A challenge. Ecol. Indic. 2017, 82, 357–366. [Google Scholar] [CrossRef]
  8. Bréda, N.J. Ground-based measurements of leaf area index: A review of methods, instruments and current controversies. J. Exp. Bot. 2003, 54, 2403–2417. [Google Scholar] [CrossRef]
  9. Korhonen, L.; Packalen, P.; Rautiainen, M. Comparison of Sentinel-2 and Landsat 8 in the estimation of boreal forest canopy cover and leaf area index. Remote Sens. Environ. 2017, 195, 259–274. [Google Scholar] [CrossRef]
  10. Gonzalez-Sanpedro, M.C.; Le Toan, T.; Moreno, J.; Kergoat, L.; Rubio, E. Seasonal variations of leaf area index of agricultural fields retrieved from Landsat data. Remote Sens. Environ. 2008, 112, 810–824. [Google Scholar] [CrossRef] [Green Version]
  11. Yao, Y.J.; Liu, Q.H.; Liu, Q.; Li, X.W. LAI retrieval and uncertainty evaluations for typical row-planted crops at different growth stages. Remote Sens. Environ. 2008, 112, 94–106. [Google Scholar] [CrossRef]
  12. Shang, J.L.; Liu, J.G.; Huffman, T.; Qian, B.D.; Pattey, E.; Wang, J.F.; Zhao, T.; Geng, X.Y.; Kroetsch, D.; Dong, T.F.; et al. Estimating plant area index for monitoring crop growth dynamics using Landsat-8 and RapidEye images. J. Appl. Remote Sens. 2014, 8, 085196. [Google Scholar] [CrossRef] [Green Version]
  13. Fang, H.L.; Liang, S.L. Retrieving leaf area index with a neural network method: Simulation and validation. IEEE Trans. Geosci. Remote Sens. 2003, 41, 2052–2062. [Google Scholar] [CrossRef] [Green Version]
  14. Tian, J.Y.; Wang, L.; Li, X.J.; Gong, H.L.; Shi, C.; Zhong, R.F.; Liu, X.M. Comparison of UAV and WorldView-2 imagery for mapping leaf area index of mangrove forest. Int. J. Appl. Earth Obs. Geoinf. 2017, 61, 22–31. [Google Scholar] [CrossRef]
  15. Comba, L.; Biglia, A.; Aimonino, D.R.; Tortia, C.; Mania, E.; Guidoni, S.; Gay, P. Leaf Area Index evaluation in vineyards using 3D point clouds from UAV imagery. Precis. Agric. 2019, 1–16. [Google Scholar] [CrossRef] [Green Version]
  16. Kimes, D.S.; Knyazikhin, Y.; Privette, J.L.; Abuelgasim, A.A.; Gao, F. Inversion methods for physically-based models. Remote Sens. Rev. 2000, 18, 381–439. [Google Scholar] [CrossRef]
  17. Dorigo, W.A.; Zurita-Milla, R.; de Wit, A.J.W.; 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 Obs. Geoinf. 2007, 9, 165–193. [Google Scholar] [CrossRef]
  18. Baret, F.; Buis, S. Estimating Canopy Characteristics from Remote Sensing Observations: Review of Methods and Associated Problems. In Advances in Land Remote Sensing: System, Modeling, Inversion and Application; Liang, S., Ed.; Springer: Dordrecht, The Netherlands, 2008; pp. 173–201. [Google Scholar]
  19. Haboudane, D.; Miller, J.R.; Pattey, E.; Zarco-Tejada, P.J.; Strachan, I.B. Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: Modeling and validation in the context of precision agriculture. Remote Sens. Environ. 2004, 90, 337–352. [Google Scholar] [CrossRef]
  20. Darvishzadeh, R.; Skidmore, A.; Atzberger, C.; van Wieren, S. Estimation of vegetation LAI from hyperspectral reflectance data: Effects of soil type and plant architecture. Int. J. Appl. Earth Obs. Geoinf. 2008, 10, 358–373. [Google Scholar] [CrossRef]
  21. Colombo, R.; Merom, M.; Marchesi, A.; Busetto, L.; Rossini, M.; Giardino, C.; Panigada, C. Estimation of leaf and canopy water content in poplar plantations by means of hyperspectral indices and inverse modeling. Remote Sens. Environ. 2008, 112, 1820–1834. [Google Scholar] [CrossRef]
  22. Eitel, J.U.H.; Gessler, P.E.; Smith, A.M.S.; Robberecht, R. Suitability of existing and novel spectral indices to remotely detect water stress in Populus spp. For. Ecol. Manag. 2006, 229, 170–182. [Google Scholar] [CrossRef]
  23. Meroni, M.; Colombo, R.; Panigada, C. Inversion of a radiative transfer model with hyperspectral observations for LAI mapping in poplar plantations. Remote Sens. Environ. 2004, 92, 195–206. [Google Scholar] [CrossRef]
  24. Combal, B.; Baret, F.; Weiss, M.; Trubuil, A.; Mace, D.; Pragnere, A.; Myneni, R.; Knyazikhin, Y.; Wang, L. Retrieval of canopy biophysical variables from bidirectional reflectance - Using prior information to solve the ill-posed inverse problem. Remote Sens. Environ. 2003, 84, 1–15. [Google Scholar] [CrossRef]
  25. Schlerf, M.; Atzberger, C. Inversion of a forest reflectance model to estimate structural canopy variables from hyperspectral remote sensing data. Remote Sens. Environ. 2006, 100, 281–294. [Google Scholar] [CrossRef]
  26. Jacquemoud, S.; Baret, F.; Andrieu, B.; Danson, F.M.; Jaggard, K. Extraction of Vegetation Biophysical Parameters by Inversion of the Prospect Plus Sail Models on Sugar-Beet Canopy Reflectance Data—Application To Tm and Aviris Sensors. Remote Sens. Environ. 1995, 52, 163–172. [Google Scholar] [CrossRef]
  27. Weiss, M.; Baret, F.; Myneni, R.B.; Pragnere, A.; Knyazikhin, Y. Investigation of a model inversion technique to estimate canopy biophysical variables from spectral and directional reflectance data. Agronomie 2000, 20, 3–22. [Google Scholar] [CrossRef]
  28. Richter, K.; Atzberger, C.; Vuolo, F.; Weihs, P.; D’Urso, G. Experimental assessment of the Sentinel-2 band setting for RTM-based LAI retrieval of sugar beet and maize. Can. J. Remote Sens. 2009, 35, 230–247. [Google Scholar] [CrossRef]
  29. Chen, J.M.; Liu, J.; Leblanc, S.G.; Lacaze, R.; Roujean, J.L. Multi-angular optical remote sensing for assessing vegetation structure and carbon absorption. Remote Sens. Environ. 2003, 84, 516–525. [Google Scholar] [CrossRef]
  30. Darvishzadeh, R.; Atzberger, C.; Skidmore, A.; Schlerf, M. Mapping grassland leaf area index with airborne hyperspectral imagery: A comparison study of statistical approaches and inversion of radiative transfer models. ISPRS J. Photogramm. Remote Sens. 2011, 66, 894–906. [Google Scholar] [CrossRef]
  31. Duan, S.-B.; Li, Z.-L.; Wu, H.; Tang, B.-H.; Ma, L.; Zhao, E.; Li, C. Inversion of the PROSAIL model to estimate leaf area index of maize, potato, and sunflower fields from unmanned aerial vehicle hyperspectral data. Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 12–20. [Google Scholar] [CrossRef]
  32. Li, H.; Chen, Z.-X.; Jiang, Z.-W.; Wu, W.-B.; Ren, J.-Q.; Liu, B.; Tuya, H. Comparative analysis of GF-1, HJ-1, and Landsat-8 data for estimating the leaf area index of winter wheat. J. Integr. Agric. 2017, 16, 266–285. [Google Scholar] [CrossRef]
  33. Li, H.; Liu, G.; Liu, Q.; Chen, Z.; Huang, C. Retrieval of Winter Wheat Leaf Area Index from Chinese GF-1 Satellite Data Using the PROSAIL Model. Sensors (Basel) 2018, 18, 1120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Xiao, Z.Q.; Liang, S.L.; Wang, J.D.; Song, J.L.; Wu, X.Y. A Temporally Integrated Inversion Method for Estimating Leaf Area Index from MODIS Data. IEEE Trans. Geosci. Remote Sens. 2009, 47, 2536–2545. [Google Scholar] [CrossRef]
  35. Verger, A.; Baret, F.; Camacho, F. Optimal modalities for radiative transfer-neural network estimation of canopy biophysical characteristics: Evaluation over an agricultural area with CHRIS/PROBA observations. Remote Sens. Environ. 2011, 115, 415–426. [Google Scholar] [CrossRef]
  36. Verrelst, J.; Alonso, L.; Camps-Valls, G.; Delegido, J.; Moreno, J. Retrieval of Vegetation Biophysical Parameters Using Gaussian Process Techniques. IEEE Trans. Geosci. Remote Sens. 2012, 50, 1832–1843. [Google Scholar] [CrossRef]
  37. Karimi, S.; Sadraddini, A.A.; Nazemi, A.H.; Xu, T.; Fard, A.F. Generalizability of gene expression programming and random forest methodologies in estimating cropland and grassland leaf area index. Comput. Electron. Agric. 2018, 144, 232–240. [Google Scholar] [CrossRef]
  38. Durbha, S.S.; King, R.L.; Younan, N.H. Support vector machines regression for retrieval of leaf area index from multiangle imaging spectroradiometer. Remote Sens. Environ. 2007, 107, 348–361. [Google Scholar] [CrossRef]
  39. Bacour, C.; Baret, F.; Béal, D.; Weiss, M.; Pavageau, K. Neural network estimation of LAI, fAPAR, fCover and LAI × C ab, from top of canopy MERIS reflectance data: Principles and validation. Remote Sens. Environ. 2006, 105, 313–325. [Google Scholar] [CrossRef]
  40. Verrelst, J.; Rivera, J.P.; Gitelson, A.; Delegido, J.; Moreno, J.; Camps-Valls, G. Spectral band selection for vegetation properties retrieval using Gaussian processes regression. Int. J. Appl. Earth Obs. Geoinf. 2016, 52, 554–567. [Google Scholar] [CrossRef]
  41. Lazaro-Gredilla, M.; Titsias, M.K.; Verrelst, J.; Camps-Valls, G. Retrieval of Biophysical Parameters With Heteroscedastic Gaussian Processes. IEEE Geosci. Remote Sens Lett. 2014, 11, 838–842. [Google Scholar] [CrossRef]
  42. Yao, X.; Wang, N.; Liu, Y.; Cheng, T.; Tian, Y.; Chen, Q.; Zhu, Y. Estimation of Wheat LAI at Middle to High Levels Using Unmanned Aerial Vehicle Narrowband Multispectral Imagery. Remote Sens. 2017, 9, 1304. [Google Scholar] [CrossRef] [Green Version]
  43. Wang, L.; Dong, T.; Zhang, G.; Niu, Z. LAI Retrieval Using PROSAIL Model and Optimal Angle Combination of Multi-Angular Data in Wheat. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 1730–1736. [Google Scholar] [CrossRef]
  44. Yang, F.; Sun, J.L.; Fang, H.L.; Yao, Z.F.; Zhang, J.H.; Zhu, Y.Q.; Song, K.S.; Wang, Z.M.; Hu, M.G. Comparison of different methods for corn LAI estimation over northeastern China. Int. J. Appl. Earth Obs. Geoinf. 2012, 18, 462–471. [Google Scholar]
  45. Koetz, B.; Baret, F.; Poilve, H.; Hill, J. Use of coupled canopy structure dynamic and radiative transfer models to estimate biophysical canopy characteristics. Remote Sens. Environ. 2005, 95, 115–124. [Google Scholar] [CrossRef]
  46. Combal, B.; Baret, F.; Weiss, M. Improving canopy variables estimation from remote sensing data by exploiting ancillary information. Case study on sugar beet canopies. Agronomie 2002, 22, 205–215. [Google Scholar]
  47. Yun-Ling, L.I.; Yan, D.H.; Pei, Y.S.; Qin, D.Y. Dynamic variation of landscape in Heihe River Basin. J. Hehai Univ. Nat. Sci. 2005, 33, 6–10. [Google Scholar]
  48. Geng, L.M.M.; Song, Y.; Liang, J.; Li, Y.; Zhou, S. HiWATER: Dataset of Leaf Area Index by LAI2200 in the lower reaches of the Heihe River Basin. Cold Arid Reg. Sci. Data Center Lanzhou 2015. [Google Scholar] [CrossRef]
  49. Gitelson, A.; Merzlyak, M.N. Quantitative Estimation of Chlorophyll-a Using Reflectance Spectra—Experiments with Autumn Chestnut and Maple Leaves. J. Photochem. Photobiol. B-Biol. 1994, 22, 247–252. [Google Scholar] [CrossRef]
  50. Gitelson, A.; Merzlyak, M.N. Spectral Reflectance Changes Associated with Autumn Senescence of Aesculus-Hippocastanum L and Acer-Platanoides L Leaves—Spectral Features and Relation To Chlorophyll Estimation. J. Plant Physiol. 1994, 143, 286–292. [Google Scholar] [CrossRef]
  51. Rouse, J.W.; Haas, R.H. Monitoring Vegetation Systems in the Great Plain with ERTS. Third ERTS Symposium 1973, 1, 309–317. [Google Scholar]
  52. Richardson, A.J.; Wiegand, C.L. Distinguishing Vegetation from Soil Background Information. Photogramm. Eng. Remote Sens. 1977, 43, 1541–1552. [Google Scholar]
  53. Broge, N.H.; Leblanc, E. Comparing prediction power and stability of broadband and hyperspectral vegetation indices for estimation of green leaf area index and canopy chlorophyll density. Remote Sens. Environ. 2001, 76, 156–172. [Google Scholar] [CrossRef]
  54. Jiang, Z.Y.; 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]
  55. Gitelson, A.A.; Kaufman, Y.J.; Merzlyak, M.N. Use of a green channel in remote sensing of global vegetation from EOS-MODIS. Remote Sens. Environ. 1996, 58, 289–298. [Google Scholar] [CrossRef]
  56. Gitelson, A.A.; Gritz, Y.; Merzlyak, M.N. Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. J. Plant Physiol. 2003, 160, 271–282. [Google Scholar] [CrossRef]
  57. Daughtry, C.S.T.; Walthall, C.L.; Kim, M.S.; de Colstoun, E.B.; McMurtrey, J.E. Estimating corn leaf chlorophyll concentration from leaf and canopy reflectance. Remote Sens. Environ. 2000, 74, 229–239. [Google Scholar] [CrossRef]
  58. Yang, Z.; Willis, P.; Mueller, R. Impact of Band-Ratio Enhanced AWIFS Image on Crop Classification Accuracy. J. Am. Chem. Soc. 2008, 136, 355–359. [Google Scholar]
  59. Qi, J.; Chehbouni, A.; Huete, A.R.; Kerr, Y.H.; Sorooshian, S. A Modified Soil Adjusted Vegetation Index. Remote Sens. Environ. 1994, 48, 119–126. [Google Scholar] [CrossRef]
  60. Verhoef, W. Light-Scattering by Leaf Layers with Application to Canopy Reflectance Modeling—The Sail Model. Remote Sens. Environ. 1984, 16, 125–141. [Google Scholar] [CrossRef] [Green Version]
  61. Nigam, R.; Bhattacharya, B.K.; Vyas, S.; Oza, M.P. Retrieval of wheat leaf area index from AWiFS multispectral data using canopy radiative transfer simulation. Int. J. Appl. Earth Obs. Geoinf. 2014, 32, 173–185. [Google Scholar] [CrossRef]
  62. Liang, L.; Di, L.; Zhang, L. Estimation of crop LAI using hyperspectral vegetation indices and a hybrid inversion method. Remote Sens. Environ. 2015, 165, 123–134. [Google Scholar] [CrossRef]
  63. Liang, L.; Qin, Z.; Zhao, S. Estimating crop chlorophyll content with hyperspectral vegetation indices and the hybrid inversion method. Int. J. Remote Sens. 2016, 37, 2923–2949. [Google Scholar] [CrossRef]
  64. Rivera-Caicedo, J.P.; Verrelst, J.; Muñoz-Marí, J.; Camps-Valls, G.; Moreno, J. Hyperspectral dimensionality reduction for biophysical variable statistical retrieval. ISPRS J. Photogramm. Remote Sens. 2017, 132, 88–101. [Google Scholar] [CrossRef]
  65. Javier Garcia-Haro, F.; Campos-Taberner, M.; Munoz-Mari, J.; Laparra, V.; Camacho, F.; Sanchez-Zapero, J.; Camps-Valls, G. Derivation of global vegetation biophysical parameters from EUMETSAT Polar System. ISPRS J. Photogramm. Remote Sens. 2018, 139, 57–74. [Google Scholar] [CrossRef]
  66. Richter, K.; Atzberger, C.; Vuolo, F.; D’Urso, G. Evaluation of Sentinel-2 Spectral Sampling for Radiative Transfer Model Based LAI Estimation of Wheat, Sugar Beet, and Maize. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2011, 4, 458–464. [Google Scholar] [CrossRef]
  67. Tanaka, S.; Kawamura, K.; Maki, M.; Muramoto, Y.; Yoshida, K.; Akiyama, T. Spectral Index for Quantifying Leaf Area Index of Winter Wheat by Field Hyperspectral Measurements: A Case Study in Gifu Prefecture, Central Japan. Remote Sen. 2015, 7, 5329–5346. [Google Scholar] [CrossRef] [Green Version]
  68. Camps-Valls, G.; Gomez-Chova, L.; Muoz-Mari, J.; Vila-Frances, J.; Amoros, J.; Valle-Tascon, S.D.; Calpe-Maravilla, J. Biophysical parameter estimation with adaptive Gaussian Processes. In Proceedings of the 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 12–17 July 2009. [Google Scholar]
  69. Verrelst, J.; Rivera, J.P.; Moreno, J.; Camps-Valls, G. Gaussian processes uncertainty estimates in experimental Sentinel-2 LAI and leaf chlorophyll content retrieval. ISPRS J. Photogramm. Remote Sens. 2013, 86, 157–167. [Google Scholar] [CrossRef]
  70. Svendsen, D.H.; Martino, L.; Campos-Taberner, M.; García-Haro, F.J.; Camps-Valls, G. Joint Gaussian processes for biophysical parameter retrieval. IEEE Trans. Geosci. Remote Sens. 2017, 56, 1718–1727. [Google Scholar] [CrossRef] [Green Version]
  71. Djamai, N.; Fernandes, R.; Weiss, M.; McNairn, H.; Goïta, K. Validation of the Sentinel Simplified Level 2 Product Prototype Processor (SL2P) for mapping cropland biophysical variables using Sentinel-2/MSI and Landsat-8/OLI data. Remote Sens. Environ. 2019, 225, 416–430. [Google Scholar] [CrossRef]
  72. Pasqualotto, N.; Delegido, J.; Van Wittenberghe, S.; Rinaldi, M.; Moreno, J. Multi-Crop Green LAI Estimation with a New Simple Sentinel-2 LAI Index (SeLI). Sensors 2019, 19, 904. [Google Scholar] [CrossRef] [Green Version]
  73. Hu, Q.; Yang, J.; Xu, B.; Huang, J.; Memon, M.S.; Yin, G.; Zeng, Y.; Zhao, J.; Liu, K. Evaluation of Global Decametric-Resolution LAI, FAPAR and FVC Estimates Derived from Sentinel-2 Imagery. Remote Sens. 2020, 12, 912. [Google Scholar] [CrossRef] [Green Version]
  74. Fang, H.L.; Liang, S.L. A hybrid inversion method for mapping leaf area index from MODIS data: Experiments and application to broadleaf and needleleaf canopies. Remote Sens. Environ. 2005, 94, 405–424. [Google Scholar] [CrossRef]
  75. Zhao, J.; Li, J.; Liu, Q.H.; Fan, W.J.; Zhong, B.; Wu, S.L.; Yang, L.; Zeng, Y.L.; Xu, B.D.; Yin, G.F. Leaf Area Index Retrieval Combining HJ1/CCD and Landsat8/OLI Data in the Heihe River Basin, China. Remote Sens. 2015, 7, 6862–6885. [Google Scholar] [CrossRef] [Green Version]
  76. Zhang, W.; Chen, Y.; Hu, S. Retrieving LAI in the Heihe and the Hanjiang river basins using landsat images for accuracy evaluation on MODIS LAI product. In Proceedings of the 2007 IEEE International Geoscience and Remote Sensing Symposium, Barcelona, Spain, 23–28 July 2007. [Google Scholar]
Sensors 20 02460 g001 550
Figure 1. Location of the study area and spatial distribution of the sample plots.
Figure 1. Location of the study area and spatial distribution of the sample plots.
Sensors 20 02460 g001
Sensors 20 02460 g002 550
Figure 2. Flow chart of leaf area index (LAI) inversion for this study.
Figure 2. Flow chart of leaf area index (LAI) inversion for this study.
Sensors 20 02460 g002
Sensors 20 02460 g003 550
Figure 3. Accuracy of LAI inversion accuracy with a different number of training datasets based on Gaussian process regression (GPR).
Figure 3. Accuracy of LAI inversion accuracy with a different number of training datasets based on Gaussian process regression (GPR).
Sensors 20 02460 g003
Sensors 20 02460 g004 550
Figure 4. Comparison of LUT and GPR with different bands of reflectance to predict LAI based on GF-1 data ((a) LUT and (b) GPR).
Figure 4. Comparison of LUT and GPR with different bands of reflectance to predict LAI based on GF-1 data ((a) LUT and (b) GPR).
Sensors 20 02460 g004
Sensors 20 02460 g005 550
Figure 5. Comparison of LUT and GPR with different VIs to estimate LAI based on GF-1 data ((a) LUT and (b) GPR).
Figure 5. Comparison of LUT and GPR with different VIs to estimate LAI based on GF-1 data ((a) LUT and (b) GPR).
Sensors 20 02460 g005
Sensors 20 02460 g006 550
Figure 6. LAI map of the Heihe study area estimated by the modified soil adjusted vegetation index (MSAVI) strategy of GPR based on the measured GF-1 data on 29 July, 2014.
Figure 6. LAI map of the Heihe study area estimated by the modified soil adjusted vegetation index (MSAVI) strategy of GPR based on the measured GF-1 data on 29 July, 2014.
Sensors 20 02460 g006
Table
Table 1. Key characteristics of GF-1 Wild Field Camera (WFV).
Table 1. Key characteristics of GF-1 Wild Field Camera (WFV).
BandWavelength Range (nm)Radiometric Resolution (bit)Spatial Resolution (m)Breadth (km)Revisit Period (d)
Blue (Band 1)450–5201016200 (1 CCD)
800 (4 CCD)		4
Green (Band 2)520–590
Red (Band 3)630–690
Near-infrared (Band 4)770–890
Table
Table 2. The 15 combinations of four bands of GF-1 WFV.
Table 2. The 15 combinations of four bands of GF-1 WFV.
No.BandNo.BandNo.Band
1B16B1, B311B1, B2, B3
2B27B1, B412B1, B2, B4
3B38B2, B313B1, B3, B4
4B49B2, B414B2, B3, B4
5B1, B210B3, B415B1, B2, B3, B4
Note: B1, B2, B3, and B4 are the wavelengths of blue, green, red, and near-infrared, respectively.
Table
Table 3. Selected ten vegetation indices (Vis) from previous studies.
Table 3. Selected ten vegetation indices (Vis) from previous studies.
No.VIFormulaReference
1NDVI(B4 − B3)/(B4 + B3)[49,50,51]
2DVIB4 − B3[52]
3TVI0.5 (120 (B4 − B2) − 200 (B3 − B2))[53]
4EVI22.5(B4 − B3)/[(B4 + 2.4B3) + 1][54]
5GNDVI(B4 − B2)/(B4 + B2)[55]
6GRVIB4/B2 − 1[56]
7MCARI 1.5 [ 2.5 ( B 4 − B 3 ) − 1.3 ( B 4 − B 2 ) ] ( ( 2 B 4 + 1 ) 2 − ( 6 B 4 − 5 B 3 ) − 0.5 ) [57]
8MNLI1.5 (B42 − B3)/(B42 + B3 +0.5)[58]
9MSAVI ( 2 B 4 + 1 ) − ( 2 B 4 + 1 ) 2 − 8 ( B 4 − B 3 ) 2 [59]
10MTVI2 1.5 [ ( 1.2 ( B 4 − B 2 ) − 2.5 ( B 3 − B 2 ) ) ] ( 2 B 4 + 1 ) 2 − ( 6 B 4 − 5 B 3 ) − 0.5 [19]
Table
Table 4. Ranges and distributions of PROSAIL input parameters for the look-up table (LUT) generation.
Table 4. Ranges and distributions of PROSAIL input parameters for the look-up table (LUT) generation.
ParameterVariablesUnitMaxMinAverageStd.Type
LeafN—2.511.51Gaussian
Cabμg.cm-29005040Gaussian
Carμg.cm-2200107Gaussian
Cbrown—1.500.20.8Gaussian
Cwcm0.0500.020.025Gaussian
Cmg.cm-20.0200.010.01Gaussian
CanopyLAIm2/m2703.52.5Gaussian
ALIAdegree80306020Gaussian
hspot—100.450.6Gaussian
Soilpsoil—100.50.5Gaussian
Solar and Sensorttsdegree7025——Fixed
ttodegree800——Fixed
psidegree120-120——Fixed

Share and Cite

MDPI and ACS Style

Zhang, Y.; Yang, J.; Liu, X.; Du, L.; Shi, S.; Sun, J.; Chen, B. Estimation of Multi-Species Leaf Area Index Based on Chinese GF-1 Satellite Data Using Look-Up Table and Gaussian Process Regression Methods. Sensors 2020, 20, 2460. https://doi.org/10.3390/s20092460

AMA Style

Zhang Y, Yang J, Liu X, Du L, Shi S, Sun J, Chen B. Estimation of Multi-Species Leaf Area Index Based on Chinese GF-1 Satellite Data Using Look-Up Table and Gaussian Process Regression Methods. Sensors. 2020; 20(9):2460. https://doi.org/10.3390/s20092460

Chicago/Turabian Style

Zhang, Yangyang, Jian Yang, Xiuguo Liu, Lin Du, Shuo Shi, Jia Sun, and Biwu Chen. 2020. "Estimation of Multi-Species Leaf Area Index Based on Chinese GF-1 Satellite Data Using Look-Up Table and Gaussian Process Regression Methods" Sensors 20, no. 9: 2460. https://doi.org/10.3390/s20092460

APA Style

Zhang, Y., Yang, J., Liu, X., Du, L., Shi, S., Sun, J., & Chen, B. (2020). Estimation of Multi-Species Leaf Area Index Based on Chinese GF-1 Satellite Data Using Look-Up Table and Gaussian Process Regression Methods. Sensors, 20(9), 2460. https://doi.org/10.3390/s20092460

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