Next Article in Journal
Multi-National Banknote Classification Based on Visible-light Line Sensor and Convolutional Neural Network
Next Article in Special Issue
Phenoliner: A New Field Phenotyping Platform for Grapevine Research
Previous Article in Journal
A Study of Vicon System Positioning Performance
Previous Article in Special Issue
Optical Sensing to Determine Tomato Plant Spacing for Precise Agrochemical Application: Two Scenarios
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sensors
Volume 17
Issue 7
10.3390/s17071593
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

Leaf Area Index Estimation Using Chinese GF-1 Wide Field View Data in an Agriculture Region

by
Xiangqin Wei
1,2,3,
Xingfa Gu
1,2,3,*,
Qingyan Meng
1,3,
Tao Yu
1,3,
Xiang Zhou
1,3,
Zheng Wei
3,
Kun Jia
4,* and
Chunmei Wang
1,3
1
Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Application Technology Center of China High-Resolution Earth Observation System, Beijing 100101, China
4
State Key Laboratory of Remote Sensing Science, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
*
Authors to whom correspondence should be addressed.
Sensors 2017, 17(7), 1593; https://doi.org/10.3390/s17071593
Submission received: 26 April 2017 / Revised: 20 June 2017 / Accepted: 5 July 2017 / Published: 8 July 2017
(This article belongs to the Special Issue Sensors in Agriculture)
Browse Figures
Versions Notes

Abstract

:
Leaf area index (LAI) is an important vegetation parameter that characterizes leaf density and canopy structure, and plays an important role in global change study, land surface process simulation and agriculture monitoring. The wide field view (WFV) sensor on board the Chinese GF-1 satellite can acquire multi-spectral data with decametric spatial resolution, high temporal resolution and wide coverage, which are valuable data sources for dynamic monitoring of LAI. Therefore, an automatic LAI estimation algorithm for GF-1 WFV data was developed based on the radiative transfer model and LAI estimation accuracy of the developed algorithm was assessed in an agriculture region with maize as the dominated crop type. The radiative transfer model was firstly used to simulate the physical relationship between canopy reflectance and LAI under different soil and vegetation conditions, and then the training sample dataset was formed. Then, neural networks (NNs) were used to develop the LAI estimation algorithm using the training sample dataset. Green, red and near-infrared band reflectances of GF-1 WFV data were used as the input variables of the NNs, as well as the corresponding LAI was the output variable. The validation results using field LAI measurements in the agriculture region indicated that the LAI estimation algorithm could achieve satisfactory results (such as R2 = 0.818, RMSE = 0.50). In addition, the developed LAI estimation algorithm had potential to operationally generate LAI datasets using GF-1 WFV land surface reflectance data, which could provide high spatial and temporal resolution LAI data for agriculture, ecosystem and environmental management researches.

1. Introduction

Vegetation is the basic component of the terrestrial ecosystem and plays an important role in energy exchange, carbon cycling and hydrological cycling process on the earth surface [1,2,3]. Therefore, timely and accurate land surface vegetation information is of great significance for earth system science, ecological environment assessment and climate change related studies [4]. Leaf area index (LAI), generally defined as one half of the total green leaf area per unit of horizontal ground surface area [5], is an important parameter for characterizing land surface vegetation conditions [6,7,8]. LAI is an essential parameter that characterizes the density of leaves and canopy structure, which can reflect the vegetation’s ability for biophysical processes such as photosynthesis, respiration and transpiration, and is also an important input variable for the carbon cycle models, crop growth models and water cycle models [9,10,11,12,13]. Therefore, accurate estimation of LAI on regional and global scales is of great significance for ecosystem modelling, biogeochemical cycle modelling, agriculture monitoring and related applications.
Remote sensing provides the only effective means to estimate LAI at the regional and global scales for its ability of continuous observation and to provide broad and impartial earth observation data [7,8,14,15]. Currently, several global LAI products have been generated using medium to low spatial resolution satellite data, such as the MODIS LAI products [6], GEOV1 LAI Products [7], MERIS LAI products [16], GLASS LAI products [17] and GLOBCARBON LAI products [18]. However, the currently existed LAI products generally have low spatial resolutions at the kilometric level, whereas decametric spatial resolution LAI data will be better suited for applications related to agriculture monitoring, as compared to kilometric resolution LAI data, which are usually larger than the typical scales of most croplands. The wide field view (WFV) sensor on board GF-1, the first satellite of the China High-resolution Earth Observation System, can acquire multi-spectral data with high spatial and temporal resolutions [19], which are valuable data sources for dynamic LAI monitoring at regional scale. Therefore, exploring the application potential of GF-1 WFV data on land surface LAI estimation and developing the specific LAI estimation algorithm for GF-1 WFV data are urgently needed.
The key of LAI estimation using remote sensing data is how to establish the relationship between LAI and land surface reflectance data according to the radiative transfer process of photon in the vegetation canopy and its spectral response characteristics [20,21,22,23]. Many algorithms have been developed to retrieve LAI using satellite remote sensing data and generally two types of algorithms can be distinguished including empirical methods and physical model based methods [6,24,25,26,27]. The empirical LAI estimation methods are based on the statistical relationships between LAI and vegetation indices (VIs), which are calibrated using field LAI measurements and remote sensing data or simulated data from canopy radiative transfer models [28,29]. The empirical methods simplify the complex radiative transfer process in the canopy, which are simple and can obtain satisfactory LAI estimation accuracy in small regions. Furthermore, when using VIs as independent variables, the VIs can highlight vegetation information and weaken the influences of canopy shadows, soil backgrounds, atmospheric conditions and angle effects [29]. However, the empirical methods only use the VIs calculated from several band reflectances, which cannot make full use of the multi-band spectral information of remote sensing data. Furthermore, the reduction of the VIs from multi-band reflectances to an index also reduces the LAI estimation constraints and increases the uncertainty of LAI estimation results, resulting in the empirical relationship between LAI and VIs changing with sensor types, vegetation types, time and geographical areas. Therefore, the application of empirical LAI estimation methods to a large scale is very difficult, due to the complexity of land surface.
The physical model based methods are mainly based on the simulation of radiative transfer in vegetation canopy and establishing the physical relationship between canopy reflectances and LAI, and then the canopy radiative transfer models are inversed [30,31]. The direct inversion of radiative transfer models is very difficult due to the complexity of the models, and iterative optimization (OPT) method [32], lookup table (LUT) method [6] and machine learning methods [26] are usually used for indirect inversion of physical models to achieve LAI estimation. The OPT method is based on the iterative minimization of a cost function which requires hundreds of runs of the canopy radiative transfer model for each pixel and therefore computationally too demanding. For practical applications, LUT method and machine learning methods are popular LAI estimation methods, which are based on the database simulated by physical models. The LUT method is conceptually the simplest technique by finding the solution for a given set of reflectance measurements, which consisted of selecting the closest cases in the database according to a cost function, and then based on extracting the corresponding set of LAI [6,33]. However, the LUT method usually requires a fixed number of input variables, unless there are very large lookup tables, which are difficult to manipulate. Machine learning methods can efficiently and accurately approximate the complex nonlinear functions, and train the algorithm parameters through training samples to realize the efficient and accurate estimation of LAI using remote sensing data [26,34]. Commonly used machine learning methods mainly include neural networks (NNs), as well as support vector machines and decision trees [35]. NNs trained over radiative transfer model simulations have been applied, with success, to estimate LAI from several sensors’ data, leading to several operational LAI production algorithms, such as the CYCLOPES and MERIS LAI products [26,36]. Therefore, based on the reality of work in the field of LAI estimation using remote sensing data, the NNs inversion of physical models is a potential and reasonable choice for LAI estimation from GF-1WFV data.
Therefore, the objective of this study is to develop a general LAI estimation algorithm for GF-1 WFV reflectance data under various land surface conditions using NNs trained over radiative transfer model simulations. Meanwhile, the field LAI measurements in an agriculture region with maize as the dominated crop type are used to assess the performance of GF-1 WFV data on LAI estimation using the proposed algorithm.

2. Study Area

The Shenzhou county of China (Figure 1) is selected as the study area to investigate the performance of LAI estimation algorithm for GF-1 WFV data. The study area is located on the North China Plain (centred at 115°35′ E, 37°53′ N) covering approximately 20 km × 25 km. The study area belongs to the temperate climate zone and is a typical upland field agriculture area in the North China Plain. The annual average temperature is approximately 13.4 °C, and the annual average precipitation is approximately 486 mm, which is mainly concentrated in July and August. It is relatively flat farmland with an average altitude of about 20 m above the sea level, so that uncertainties of LAI estimation caused by topographical facts will be reduced to a minimum. The main autumn grain crops planted in the study area are maize. Maize season begins in mid-June and harvests in early October. Though the study area is not big, it has the representative characteristics of crop type distribution in the North China Plain.

3. Data and Pre-Processes

3.1. Field LAI Measurements

In order to validate the LAI estimation accuracy using GF-1 WFV data, the field campaigns were completed covering the main growing seasons of maize, and the specific field observation dates were 27 June, 21 July, 14 August and 5 September 2014. Twenty-three representative square sample sites (30 m × 30 m) within the study area (Figure 1) were selected based on the crop growth conditions and the high spatial resolution remote sensing data. The sample sites were located in relatively homogeneous regions with approximately 50 m around the sample sites having similar crop growth conditions, thus the uncertainty of LAI estimation caused by the pixel matching error between GF-1 WFV data and the sample site would be minimized. The center of each sample site was determined using the handheld global positioning system (GPS) receiver which had a positioning accuracy of approximately ±3 m. LAI was measured using an LAI-2000 plant canopy analyzer (Li-Cor, Inc., Lincoln, NE, USA), and there were three measuring plots in each sample site and the average LAI value of the three LAI measurements was regarded as the LAI of the sample site. The number of LAI measurement sites on 27 June was four because maize in this period was very small, and on 21 July, 14 August and 5 September were all twenty-three. These field LAI measurements were used to validate the LAI estimates using the GF-1 WFV data.

3.2. GF-1 WFV Data and Pre-Processes

In order to validate the proposed LAI estimation algorithm, GF-1 WFV data synchronized with the field LAI measurements were required. The spatial resolution of GF-1 WFV data was 16 m and the temporal resolution was four days for the four WFV sensors combined. According to the field survey dates and the image quality of GF-1 WFV data, the high-quality GF-1 WFV data close to the field survey dates were selected as much as possible, and finally 5 GF-1 WFV data were collected (Table 1). The GF-1 WFV data acquired on 24 June, 18 July and 15 August 2014 were close to the corresponding field survey dates (3 days or less), which could be approximated that LAI in this close time did not change and the field LAI measurements could be directly used to validate the LAI estimates using GF-1 WFV data. However, corresponding to the field LAI measurements on 5 September 2014, there were no suitable GF-1 WFV data close to the field survey date due to the effect of clouds, thus the GF-1 WFV data acquired on 24 August and 18 September 2014 were selected to estimate LAI, which would be interpolated to estimate LAI on 5 September using the spline interpolation method and then the LAI estimation accuracy would be validated using the filed LAI measurements on 5 September.
The pre-process of GF-1 WFV data contained radiance calibration, atmospheric correction and geometric correction. The radiance calibration was converting the DN values of the raw data to radiances using the following equation:
L e   =   Gain   × DN   +   Offset
where Le was the radiance, and Gain and Offset were calibration coefficients obtained from the China Centre for Resources Satellite Data and Application, which were listed in Table 2.
The Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) model was used to atmospheric correction of GF-1 WFV data [37]. The input parameters for FLAASH model were determined based on the imaging time and imaging parameters of each GF-1 WFV data. After the atmospheric correction, the GF-1 WFV land surface reflectance data were obtained and the reflectances of green, red and near-infrared (NIR) bands would be used to estimate LAI. The geometric correction of GF-1 WFV data was conducted using two-order polynomial transformation method with bilinear interpolation resampling. High quality Landsat-8 Operational Land Imager data [38] were selected as the base map, and ground control points were selected from the images manually. The resulted geometric co-registration error for each GF-1 WFV data was less than one pixel of the Landsat data (30 m). Finally, the subset images of the atmospherically and geometrically corrected GF-1 WFV data covered the sample sites were extracted to further investigate the performance of the LAI estimation algorithm for GF-1 WFV data (Figure 1).

4. Methods

A flowchart of the LAI estimation algorithm development for GF-1 WFV data was presented in Figure 2. The training sample dataset was firstly generated using the spectral reflectance simulations based on the radiative transfer model and the relative spectral response profiles of GF-1 WFV sensors. Then, the NNs were trained over the training sample dataset to develop the LAI estimation algorithm for GF-1 WFV data. Finally, LAI could be estimated using the pre-processed GF-1 WFV land surface reflectance data using the trained NNs or set to zero when the pixel was non-vegetation which was judged by the normalized difference vegetation index (NDVI) value.

4.1. Generating LAI Training Sample Dataset from Radiative Transfer Model Simulations

The canopy radiative transfer model quantitatively described the physical relationship between LAI and canopy spectral reflectances. The widely used coupled PROSPECT and SAIL (PROSAIL) model, which was easy to use, had general robustness and also had consistent performance in validation practices [22,39,40] was selected to simulate the satellite observations of canopy reflectance based on the relative spectral response profiles of the GF-1 WFV sensors. The PROSPECT model was a plate model based radiative transfer model, which simulated the hemispherical reflectance and transmittance of leaves from spectral wavelength from 0.4 to 2.5 μm based on the leaf biochemical and biophysical parameters [41,42]. The input parameters for PROSPECT model included leaf chlorophyll a + b concentration (Cab), water content (Cw), dry matter content (Cm), brown pigment content (Cbrown), carotenoid content (Car) and leaf structure parameter (N), and the output parameters were hemispherical reflectance and transmittance of leaves, which were also the input parameters for SAIL model. The SAIL model with hot-spot correction, which assumed the canopy as a turbid medium, was selected in this study [22,26]. The canopy structure in the SAIL model was characterized by LAI, the average leaf angle inclination (ALA) assuming an ellipsoidal distribution and the hot-spot parameter.
The underlying soil reflectances were also input variables for the PROSAIL model. The soil reflectances from a globally distributed soil spectral library released by the International Soil Reference and Information Centre (access at: http://www.isric.org), which contained various soil types with various properties and had representative of various soil types [19,43], were selected for the inputs of PROSAIL model. To remove data redundancy generated by similar soil reflectances and to avoid huge computations in PROSAIL simulations, several representative soil reflectances should be determined from the original data. The spectral angle mapper method [44,45] was used to evaluate the similarity of different soil reflectances and further determined the representative soil reflectances for PROSAIL model. The soil reflectances having spectral angle value smaller than 0.05 would be considered as similar soil reflectances which would be averaged as a representative soil reflectance. Finally, 13 soil reflectances were determined to represent the possible range of soil reflectances in the PROSAIL model (Figure 3).
It had been demonstrated that reasonable error of the input variables of radiative transfer model was permitted which did not lead to obvious loss of parameter inversion accuracy, therefore some input variables in the PROSAIL model could be fixed in the simulations [46,47]. Considering the objective of this study was to develop a general LAI estimation algorithm using GF-1 WFV data under various land surface conditions, the input variables of PROSAIL model were given reasonable ranges or fixed values (Table 3) based on previous studies, such as the Leaf Optical Properties Experiment 93 and the algorithms of CYCLOPES LAI product [22,26,48,49]. Therefore, the PROSAIL model could simulate the vegetation canopy reflectances based on the vegetation physicochemical parameters, geometric parameters and soil reflectances, meanwhile the remote sensing data could also obtain vegetation canopy reflectances by atmospheric correction, thus the remote sensing reflectance data were associated with LAI through the physical process of radiation transmission.
For any combination of the input variables, the canopy reflectance was computed for each wavelength from 0.4 to 2.5 μm at a spectral resolution of 1 nm, and then the canopy reflectance was resampled to simulate GF-1 WFV observations using the relative spectral response profiles. Because the relative spectral response profiles of the 4 WFV sensors on board GF-1 satellite were slightly different, the simulations of GF-1 WFV reflectance were separately generated for each WFV sensor. The canopy reflectance simulations using PROSAIL model resulted in 2,021,760 cases of matched reflectances and LAI values for each WFV sensor. To account for uncertainties in the satellite observations and model simulations, a white Gaussian noise with signal to noise ratio of 100 was added to the simulated reflectances, which were used as the learning dataset for the NNs.

4.2. Neural Networks

The NNs mimic human learning to build relationships between variables, which can approximate multivariate non-linear relationships and are robust to noisy data [50,51], thus have been widely used for estimating land surface parameters from remote sensing data [7,52,53]. The popular back propagation NNs (BPNNs) were selected for LAI estimation algorithm development using GF-1 WFV data [4,7]. The BPNNs could learn from the training sample dataset and built relationships between reflectances and LAI, then the trained BPNNs could generate the optimal LAI estimates based on the actual reflectances from GF-1 WFV data. The architecture of the BPNNs used for LAI estimation from GF-1 WFV data was shown in Figure 4. The input variables of the BPNNs included the reflectances of the green, red and NIR bands, and the output variable was the corresponding LAI. The number of nodes in the hidden layer was set to six. The BPNNs activation functions in the hidden layer and output nodes were set to “sigmoid” and “linear”, respectively. The Levenberg–Marquardt minimization algorithm was used to calibrate the synaptic coefficients because of its efficient convergence capacity [19,54]. The training sample dataset made of pairs of reflectances and LAI, which were generated from the PROSAIL model simulations, was randomly split into two parts: 90% of the cases were used for training BPNNs, meanwhile the rest 10% of the cases were used to test the hyper-specialization during the training process. The BPNNs models were generated for each of the four GF-1 WFV sensors based on the specific training sample dataset for each WFV sensor.

4.3. LAI Estimating Procedure for GF-1 WFV Data

The NDVI generated from red and NIR band reflectances was an important indicator for vegetation growth conditions, and usually used to estimate LAI based on empirical relationships between LAI and NDVI [55,56,57]. Therefore, to avoid the presentation of abnormal NNs inversion values caused by the non-vegetation pixels, it was a good strategy to remove non-vegetation pixels using NDVI threshold value before LAI estimation using the NNs. In this study, the NDVI threshold value of 0.05 was used to distinguish vegetation and non-vegetation pixels. The pixels with NDVI value smaller than 0.05 was identified as non-vegetation pixels and LAI would be set to zero, whereas the LAI for other pixels would be estimated using the trained NNs, from which the GF-1 WFV data should be judged from that sensor and then adopt the corresponding trained NNs.

4.4. LAI Estimation Accuracy Assessment

Direct comparison of field LAI measurements with LAI estimates using GF-1 WFV land surface reflectance data is a reliable way to assess the accuracy of the proposed LAI estimation algorithm. Comparison between the field survey and predicted LAI will be performed by using five statistical indices: R2 of linear regression, root mean square error (RMSE), relative RMSE (RRMSE), Nash–Sutcliffe model efficiency coefficient (EF) [58] and coefficient of residual mass (CRM) [59]. The RMSE values show how much the predicted LAI values under- or over-estimate the field observations. The RRMSE is calculated by dividing the RMSE by the mean field observations. The EF value compares the predicted LAI values to the average value of the field observations. A negative EF value indicates that the average value of the field observations gives a better estimate than the predicted values. The CRM is a measure of the tendency of the developed LAI estimation method to overestimate or underestimate the field observations. A positive CRM value indicates that the method underestimates the field observations and the negative value indicates a tendency to overestimate. For a perfect fit between field observed and predicted LAI data, values of R2, RMSE, RRMSE, CRM, and EF should equal 1.0, 0.0, 0.0, 0.0 and 1.0, respectively.

5. Results

The LAI estimation results using the proposed algorithm and the pre-processed GF-1 WFV land surface reflectance data are shown in Figure 5. In the visual aspect, the spatial and temporal distributions of LAI estimates are reasonable. As for the spatial distribution, LAI estimates using GF-1 WFV data acquired on 27 June generally show small values, because maize at this period is at the early emergence stage with small LAI, while larger LAI estimates only present at an orchard on the upper left corner of the study area. After July, with the rapid growth of maize, the large LAI estimates are mainly distributed in the area of farmland, and the LAI estimates at other areas are smaller. Meanwhile, LAI estimates in the cloud and cloud shadow regions are all zero. Regarding the temporal variations, LAI estimates of maize show low values in late June and continuous increases in July and mid-August, reaching the peak LAI values in the late August and a downward trend in mid-September. The temporal variation characteristics of LAI estimates exactly reflect the growth processes of maize from germination, rapid growth to peak value, and then gradually decline, which is in accord with the growth characteristics of maize. Therefore, the visual observation of the spatial and temporal characteristics of LAI estimates using GF-1 WFV data could preliminarily indicate the reasonability of the proposed LAI estimation algorithm.
The field LAI measurements located in the cloud and cloud shadow regions were deleted based on the visual observation, and the remaining field LAI measurements, which contained 43 LAI data, were used to directly validate the LAI estimates from GF-1 WFV data using the proposed algorithm (Figure 6). It could be seen that there was a good linear relationship between field LAI measurements and LAI estimates from GF-1 WFV data. The LAI estimation accuracy was satisfactory with values of R2, RMSE, RRMSE, CRM, and EF equaling 0.818%, 0.5%, 19.0%, 0.04% and 0.797%, respectively. In addition, some larger differences between the field survey LAI and the predicted LAI were observed in the region of large LAI values. The main reason might be that the field LAI values measured on 5 September were matched with LAI estimates interpolated from LAI estimates on 24 August and 18 September, whereas maize in this period was close to the mature stage with great LAI change rate, thus the interpolated LAI results might have some uncertainties. The other reason might be that the clouds influenced the atmospheric conditions of the study area which further influenced the accuracy of obtaining GF-1 WFV land surface reflectances and caused the uncertainties of LAI estimates.
In summary, the difference between the field LAI measurements and LAI estimates from GF-1 WFV data was small and the LAI estimation results were satisfactory. Therefore, the results indicated that the proposed LAI estimation algorithm for GF-1 WFV data was reliable and GF-1 WFV data could achieve acceptable performance on LAI estimation, which had the potential for providing high temporal and spatial resolution LAI dataset for related applications, such as agriculture monitoring.

6. Discussion and Conclusions

This study proposed a LAI estimation algorithm for GF-1 WFV land surface reflectance data based on BPNNs with training samples generated from the radiative transfer model simulations under different soil and vegetation conditions, thus the algorithm could adapt to a variety of different underlying conditions. The validation results using field LAI measurements in an agriculture region with maize as the dominated crop type showed that the proposed algorithm could achieve satisfactory LAI estimation accuracy (e.g., RMSE = 0.50), which indicated GF-1 WFV data had good performance on LAI estimation and the algorithm had potential to operationally estimate LAI from GF-1 WFV land surface reflectances. The LAI estimation algorithm was automatically operated without prior knowledge on land cover, and no human interaction and no empirical model parameters were needed. Therefore, the proposed LAI estimation algorithm for GF-1 WFV data could overcome the difficulties in determination of model parameters for empirical LAI estimation methods, which were generally changed with time, region and vegetation types. In addition, the studies of LAI estimation using remote sensing data based on NNs trained over radiative transfer model simulations in the past mainly focused on the kilometric spatial resolution remote sensing data, and several global LAI products were generated from SPOT VEGETATION, MERIS and MODIS data. Therefore, the developed LAI estimation algorithm for decametric spatial resolution remote sensing data using NNs based on physical model was a new try in this study, and LAI with decametric spatial resolution would be more useful for agriculture, ecosystem and environment management than kilometric spatial resolution data which were usually larger than the typical scales of most landscapes.
There were also some limitations about this study, even though it achieved satisfactory LAI estimation performance. Firstly, the field LAI measurements for validation were collected from an agriculture region covering the main growing season of maize, though they had temporal–spatial representativeness of cropland and the agriculture monitoring was the main application field of LAI, exploratory investigations should be strengthened using much more validation data from various vegetation types in the future work. Moreover, the pre-processes of GF-1 WFV data were conducted manually. If an automated data preprocessing approach for radiation correction and geometric correction could be developed, a streamlined LAI estimation workflow would be formed from the original DN values to LAI estimates, which would make the LAI estimation from GF-1 WFV data much simpler and quicker. Finally, the accuracy assessment of LAI estimation algorithm was difficult to achieve due to the difficulties in obtaining real LAI values from field measurements to match with the pixel level LAI estimates. Currently, using plant canopy analyzer to obtain field LAI measurements was the commonly used method, which could obtain relatively accurate LAI data. However, some errors might be made when matching the field LAI measurements and pixel scale LAI estimates due to the non-absolute uniformity of the land surface, though averaging the multi-point measurements could reduce this error to a certain extent. Perhaps using an unmanned aircraft system flying in the lower altitude to take photos from the nadir and developing corresponding LAI extraction algorithm from photos was a considerable way for real LAI data collection. Although the field LAI measurements had many uncertainties, they still played an important role in assessing LAI estimation algorithms using remote sensing data, and more effective real LAI obtaining approaches would be expected in the future.
In conclusion, the proposed LAI estimation algorithm for GF-1 WFV data was reliable and GF-1 WFV data were confirmed as having a good performance on LAI estimation, which could provide reliable high spatial and temporal resolution LAI data for related applications. Further work would focus on validating the proposed LAI estimation algorithm using field LAI measurements with less uncertainties under various land cover type conditions.

Acknowledgments

This study was partially supported by the Major Special Project of the China High-Resolution Earth Observation System and the Science and Technology Subsidies to Leading Scientists by National Administration of Surveying, Mapping and Geoinformation.

Author Contributions

X.W., X.G. and K.J. conceived and designed the experiments; X.W. performed the experiments; X.W. and K.J. analyzed the data; X.W. wrote the paper; and X.W., X.G., Q.M., T.Y., X.Z., Z.W., K.J. and C.W. contributed to revise the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiang, B.; Liang, S.L.; Yuan, W.P. Observational evidence for impacts of vegetation change on local surface climate over northern China using the Granger causality test. J. Geophys. Res. Biogeosci. 2015, 120, 1–12. [Google Scholar] [CrossRef]
  2. Jia, K.; Liang, S.; Liu, S.; Li, Y.; Xiao, Z.; Yao, Y.; Jiang, B.; Zhao, X.; Wang, X.; Xu, S.; et al. Global land surface fractional vegetation cover estimation using general regression neural networks from MODIS surface reflectance. IEEE Trans. Geosci. Remote Sens. 2015, 53, 4787–4796. [Google Scholar] [CrossRef]
  3. Bonan, G.B. Forests and climate change: Forcings, feedbacks, and the climate benefits of forests. Science 2008, 320, 1444–1449. [Google Scholar] [CrossRef] [PubMed]
  4. Liang, S. Advances in Land Remote Sensing System, Modeling Inversion and Application; Springer: Dordrecht, The Netherlands, 2008. [Google Scholar]
  5. Chen, J.M.; Black, T.A. Defining leaf area index for non-flat leaves. Plant Cell Environ. 1992, 15, 421–429. [Google Scholar] [CrossRef]
  6. Myneni, R.B.; Hoffman, S.; Knyazikhin, Y.; Privette, J.L.; Glassy, J.; Tian, Y.; Wang, Y.; Song, X.; Zhang, Y.; Smith, G.R.; et al. Global products of vegetation leaf area and fraction absorbed PAR from year one of MODIS data. Remote Sens. Environ. 2002, 83, 214–231. [Google Scholar] [CrossRef]
  7. Baret, F.; Weiss, M.; Lacaze, R.; Camacho, F.; Makhmara, H.; Pacholcyzk, P.; Smets, B. GEOV1: LAI and FAPAR essential climate variables and FCOVER global time series capitalizing over existing products. Part 1: Principles of development and production. Remote Sens. Environ. 2013, 137, 299–309. [Google Scholar] [CrossRef]
  8. Xiao, Z.Q.; Liang, S.L.; Wang, J.D.; Chen, P.; Yin, X.J.; Zhang, L.Q.; Song, J.L. Use of General Regression Neural Networks for Generating the GLASS Leaf Area Index Product from Time-Series MODIS Surface Reflectance. IEEE Trans. Geosci. Remote Sens. 2014, 52, 209–223. [Google Scholar] [CrossRef]
  9. Myneni, R.B.; Ross, J.; Asrar, G. A review on the theory of photon transport in leaf canopies. Agric. For. Meteorol. 1989, 45, 1–153. [Google Scholar] [CrossRef]
  10. Sellers, P.J.; Dickinson, R.E.; Randall, D.A.; Betts, A.K.; Hall, F.G.; Berry, J.A.; Collatz, G.J.; Denning, A.S.; Mooney, H.A.; Nobre, C.A.; et al. Modeling the exchanges of energy, water, and carbon between continents and the atmosphere. Science 1997, 275, 502–509. [Google Scholar] [CrossRef] [PubMed]
  11. Arora, V. Modeling vegetation as a dynamic component in soil-vegetation-atmosphere transfer schemes and hydrological models. Rev. Geophys. 2002, 40. [Google Scholar] [CrossRef]
  12. Bu, H.; Sharma, L.K.; Denton, A.; Franzen, D.W. Sugar Beet Yield and Quality Prediction at Multiple Harvest Dates Using Active-Optical Sensors. Agron. J. 2016, 108, 273–284. [Google Scholar] [CrossRef]
  13. Bu, H.; Sharma, L.K.; Denton, A.; Franzen, D.W. Comparison of Satellite Imagery and Ground-Based Active Optical Sensors as Yield Predictors in Sugar Beet, Spring Wheat, Corn, and Sunflower. Agron. J. 2017, 109, 299–308. [Google Scholar] [CrossRef]
  14. Sharma, L.K.; Bu, H.; Denton, A.; Franzen, D.W. Active-Optical Sensors Using Red NDVI Compared to Red Edge NDVI for Prediction of Corn Grain Yield in North Dakota, U.S.A. Sensors 2015, 15, 27832–27853. [Google Scholar] [CrossRef] [PubMed]
  15. Sharma, L.K.; Bu, H.; Franzen, D.W. Comparison of two ground-based active-optical sensors for in-season estimation of corn (Zea mays L.) yield. J. Plant Nutr. 2016, 39, 957–966. [Google Scholar] [CrossRef]
  16. Gobron, N.; Pinty, B.; Verstraete, M.; Govaerts, Y. The MERIS Global Vegetation Index (MGVI): Description and preliminary application. Int. J. Remote Sens. 1999, 20, 1917–1927. [Google Scholar] [CrossRef]
  17. Liang, S.L.; Zhao, X.; Liu, S.H.; Yuan, W.P.; Cheng, X.; Xiao, Z.Q.; Zhang, X.T.; Liu, Q.; Cheng, J.; Tang, H.R.; et al. A long-term Global LAnd Surface Satellite (GLASS) data-set for environmental studies. Int. J. Digit. Earth 2013, 6, 5–33. [Google Scholar] [CrossRef]
  18. Plummer, S.; Arino, O.; Simon, M.; Steffen, W. Establishing a earth observation product service for the terrestrial carbon community: The globcarbon initiative. Mitig. Adapt. Strateg. Glob. Chang. 2006, 11, 97–111. [Google Scholar] [CrossRef]
  19. Jia, K.; Liang, S.; Gu, X.; Baret, F.; Wei, X.; Wang, X.; Yao, Y.; Yang, L.; Li, Y. Fractional vegetation cover estimation algorithm for Chinese GF-1 wide field view data. Remote Sens. Environ. 2016, 177, 184–191. [Google Scholar] [CrossRef]
  20. Chen, J.M.; Leblanc, S.G. A four-scale bidirectional reflectance model based on canopy architecture. IEEE Trans. Geosci. Remote Sens. 1997, 35, 1316–1337. [Google Scholar] [CrossRef]
  21. 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]
  22. Jacquemoud, S.; Verhoef, W.; Baret, F.; Bacour, C.; Zarco-Tejada, P.J.; Asner, G.P.; Francois, C.; Ustin, S.L. PROSPECT + SAIL models: A review of use for vegetation characterization. Remote Sens. Environ. 2009, 113, S56–S66. [Google Scholar] [CrossRef]
  23. Sharma, L.K.; Franzen, D.W. Use of corn height to improve the relationship between active optical sensor readings and yield estimates. Precis. Agric. 2014, 15, 331–345. [Google Scholar] [CrossRef]
  24. Liang, S. Recent developments in estimating land surface biogeophysical variables from optical remote sensing. Prog. Phys. Geogr. 2007, 31, 501–516. [Google Scholar] [CrossRef]
  25. Liang, S.; Li, X.; Wang, J. Advanced Remote Sensing: Terrestrial Information Extraction and Applications; Academic Press: Oxford, UK, 2012; p. 800. [Google Scholar]
  26. Baret, F.; Hagolle, O.; Geiger, B.; Bicheron, P.; Miras, B.; Huc, M.; Berthelot, B.; Nino, F.; Weiss, M.; Samain, O.; et al. LAI, fAPAR and fCover CYCLOPES global products derived from VEGETATION—Part 1: Principles of the algorithm. Remote Sens. Environ. 2007, 110, 275–286. [Google Scholar] [CrossRef]
  27. Sharma, L.; Bali, S.; Dwyer, J.; Plant, A.; Bhowmik, A. A Case Study of Improving Yield Prediction and Sulfur Deficiency Detection Using Optical Sensors and Relationship of Historical Potato Yield with Weather Data in Maine. Sensors 2017, 17, 1095. [Google Scholar] [CrossRef] [PubMed]
  28. Wang, F.-M.; Huang, J.-F.; Tang, Y.-L.; Wang, X.-Z. New Vegetation Index and Its Application in Estimating Leaf Area Index of Rice. Rice Sci. 2007, 14, 195–203. [Google Scholar] [CrossRef]
  29. Chen, J.M. Evaluation of Vegetation Indices and a Modified Simple Ratio for Boreal Applications. Can. J. Remote Sens. 1996, 22, 229–242. [Google Scholar] [CrossRef]
  30. 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]
  31. Gascon, F.; Gastellu-Etchegorry, J.P.; Lefevre-Fonollosa, M.J.; Dufrene, E. Retrieval of forest biophysical variables by inverting a 3-D radiative transfer model and using high and very high resolution imagery. Int. J. Remote Sens. 2004, 25, 5601–5616. [Google Scholar] [CrossRef]
  32. Bicheron, P.; Leroy, M. A Method of Biophysical Parameter Retrieval at Global Scale by Inversion of a Vegetation Reflectance Model. Remote Sens. Environ. 1999, 67, 251–266. [Google Scholar] [CrossRef]
  33. Weiss, M.; Baret, F.; Myneni, R.B.; Pragnère, 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]
  34. Roujean, J.L.; Lacaze, R. Global mapping of vegetation parameters from POLDER multiangular measurements for studies of surface-atmosphere interactions: A pragmatic method and its validation. J. Geophys. Res. Atmos. 2002, 107. [Google Scholar] [CrossRef]
  35. Yang, L.; Jia, K.; Liang, S.; Liu, J.; Wang, X. Comparison of Four Machine Learning Methods for Generating the GLASS Fractional Vegetation Cover Product from MODIS Data. Remote Sens. 2016, 8, 682. [Google Scholar] [CrossRef]
  36. García-Haro, F.J.; Camacho, F.; Meliá, J. Inter-comparison of SEVIRI/MSG and MERIS/ENVISAT biophysical products over Europe and Africa. In Proceedings of the 2nd MERIS/(A)ATSR User Workshop, ESA SP-666, Frascati, Italy, 22–26 September 2008. [Google Scholar]
  37. Cooley, T.; Anderson, G.P.; Felde, G.W.; Hoke, M.L.; Ratkowski, A.J.; Chetwynd, J.H.; Gardner, J.A.; Adler-Golden, S.M.; Matthew, M.W.; Berk, A.; et al. FLAASH, a MODTRAN4-based atmospheric correction algorithm, its application and validation. In Proceedings of the 2002 IEEE International Geoscience and Remote Sensing Symposium & 24th Canadian Symposium on Remote Sensing, Toronto, Canada, 24–28 June 2002; pp. 1414–1418. [Google Scholar]
  38. Jia, K.; Wei, X.Q.; Gu, X.F.; Yao, Y.J.; Xie, X.H.; Li, B. Land cover classification using Landsat 8 Operational Land Imager data in Beijing, China. Geocarto Int. 2014, 29, 941–951. [Google Scholar] [CrossRef]
  39. Kuusk, A. The hot spot effect in plant canopy reflectance. In Photon-Vegetation Interactions; Myneni, R.B., Ross, J., Eds.; Springer: Berlin/Heidelberg, Germany, 1991; pp. 139–159. [Google Scholar]
  40. Verhoef, W.; Jia, L.; Xiao, Q.; Su, Z. Unified optical-thermal four-stream radiative transfer theory for homogeneous vegetation canopies. IEEE Trans. Geosci. Remote Sens. 2007, 45, 1808–1822. [Google Scholar] [CrossRef]
  41. Jacquemoud, S.; Baret, F. PROSPECT—A model of leaf optical-properties spectra. Remote Sens. Environ. 1990, 34, 75–91. [Google Scholar] [CrossRef]
  42. Feret, J.-B.; François, C.; Asner, G.P.; Gitelson, A.A.; Martin, R.E.; Bidel, L.P.R.; Ustin, S.L.; le Maire, G.; Jacquemoud, S. PROSPECT-4 and 5: Advances in the leaf optical properties model separating photosynthetic pigments. Remote Sens. Environ. 2008, 112, 3030–3043. [Google Scholar] [CrossRef]
  43. Shepherd, K.D.; Palm, C.A.; Gachengo, C.N.; Vanlauwe, B. Rapid characterization of organic resource quality for soil and livestock management in tropical agroecosystems using near-infrared spectroscopy. Agron. J. 2003, 95, 1314–1322. [Google Scholar] [CrossRef]
  44. Dennison, P.E.; Halligan, K.Q.; Roberts, D.A. A comparison of error metrics and constraints for multiple endmember spectral mixture analysis and spectral angle mapper. Remote Sens. Environ. 2004, 93, 359–367. [Google Scholar] [CrossRef]
  45. Jia, K.; Li, Q.Z.; Tian, Y.C.; Wu, B.F.; Zhang, F.F.; Meng, J.H. Accuracy Improvement of Spectral Classification of Crop Using Microwave Backscatter Data. Spectrosc. Spectr. Anal. 2011, 31, 483–487. [Google Scholar]
  46. Goel, N.S.; Strebel, D.E. Inversion of vegetation canopy reflectance models for estimating agronomic variables. I. Problem definition and initial results using the suits model. Remote Sens. Environ. 1983, 13, 487–507. [Google Scholar] [CrossRef]
  47. Qu, Y.H.; Zhang, Y.Z.; Wang, J.D. A dynamic Bayesian network data fusion algorithm for estimating leaf area index using time-series data from in situ measurement to remote sensing observations. Int. J. Remote Sens. 2012, 33, 1106–1125. [Google Scholar] [CrossRef]
  48. Hosgood, B.; Jacquemoud, S.; Andreoli, G.; Verdebout, J.; Pedrini, G.; Schmuck, G. Leaf Optical Properties EXperiment 90 (LOPEX93); Report EUR 16095 EN; European Commission, Joint Research Center, Institute for Remote Sensing Applications: Brussels, Italy, 1990. [Google Scholar]
  49. Qu, Y.; Wang, J.; Wan, H.; Li, X.; Zhou, G. A Bayesian network algorithm for retrieving the characterization of land surface vegetation. Remote Sens. Environ. 2008, 112, 613–622. [Google Scholar] [CrossRef]
  50. Murthy, C.S.; Raju, P.V.; Badrinath, K.V.S. Classification of wheat crop with multi-temporal images: Performance of maximum likelihood and artificial neural networks. Int. J. Remote Sens. 2003, 24, 4871–4890. [Google Scholar] [CrossRef]
  51. Ahmad, S.; Kalra, A.; Stephen, H. Estimating soil moisture using remote sensing data: A machine learning approach. Adv. Water Resour. 2010, 33, 69–80. [Google Scholar] [CrossRef]
  52. 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]
  53. Jia, K.; Liang, S.; Wei, X.Q.; Li, Q.; Du, X.; Jiang, B.; Yao, Y.J.; Zhao, X.; Li, Y. Fractional forest cover changes in Northeast China from 1982 to 2011 and its relationship with climatic variations. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 775–783. [Google Scholar] [CrossRef]
  54. Ngia, L.S.H.; Sjoberg, J. Efficient training of neural nets for nonlinear adaptive filtering using a recursive Levenberg-Marquardt algorithm. IEEE Trans. Signal Process. 2000, 48, 1915–1927. [Google Scholar] [CrossRef]
  55. Curran, P.J.; Dungan, J.L.; Gholz, H.L. Seasonal LAI in slash pine estimated with landsat TM. Remote Sens. Environ. 1992, 39, 3–13. [Google Scholar] [CrossRef]
  56. Duan, H.; Yan, C.; Tsunekawa, A.; Song, X.; Li, S.; Xie, J. Assessing vegetation dynamics in the Three-North Shelter Forest region of China using AVHRR NDVI data. Environ. Earth Sci. 2011, 64, 1011–1020. [Google Scholar] [CrossRef]
  57. Masson, V.; Champeaux, J.L.; Chauvin, F.; Meriguet, C.; Lacaze, R. A global database of land surface parameters at 1-km resolution in meteorological and climate models. J. Clim. 2003, 16, 1261–1282. [Google Scholar] [CrossRef]
  58. Nash, J.E.; Sutcliffe, J.V. River flow forecasting through conceptual models part I—A discussion of principles. J. Hydrol. 1970, 10, 282–290. [Google Scholar] [CrossRef]
  59. Loague, K.; Green, R.E. Statistical and graphical methods for evaluating solute transport models: Overview and application. J. Contam. Hydrol. 1991, 7, 51–73. [Google Scholar] [CrossRef]
Sensors 17 01593 g001 550
Figure 1. Square region in the left image shows the geo-location of the Shenzhou study area in Hebei Province, and the right image is the GF-1 WFV data acquired on 15 August 2014. Generally, the red patches are farmland in the right image and the blue patches are non-vegetated regions, such as residential areas and roads. The green triangles on the right image indicate the field survey locations.
Figure 1. Square region in the left image shows the geo-location of the Shenzhou study area in Hebei Province, and the right image is the GF-1 WFV data acquired on 15 August 2014. Generally, the red patches are farmland in the right image and the blue patches are non-vegetated regions, such as residential areas and roads. The green triangles on the right image indicate the field survey locations.
Sensors 17 01593 g001
Sensors 17 01593 g002 550
Figure 2. Flowchart of the leaf area index (LAI) estimation algorithm for GF-1 WFV data.
Figure 2. Flowchart of the leaf area index (LAI) estimation algorithm for GF-1 WFV data.
Sensors 17 01593 g002
Sensors 17 01593 g003 550
Figure 3. The 13 soil reflectances used to represent the possible range of spectral shapes for the PROSAIL model.
Figure 3. The 13 soil reflectances used to represent the possible range of spectral shapes for the PROSAIL model.
Sensors 17 01593 g003
Sensors 17 01593 g004 550
Figure 4. The architecture of the back propagation neural networks (BPNNs) used for LAI estimation from GF-1 WFV reflectance data.
Figure 4. The architecture of the back propagation neural networks (BPNNs) used for LAI estimation from GF-1 WFV reflectance data.
Sensors 17 01593 g004
Sensors 17 01593 g005 550
Figure 5. GF-1 WFV land surface reflectance data ((a) 27 June 2014, (b) 18 July 2014, (c) 15 August 2014, (d) 24 August 2014, and (e) 18 September 2014) and their corresponding LAI estimates ((f): 27 June 2014, (g) 18 July 2014, (h) 15 August 2014, (i) 24 August 2014, and (j) 18 September 2014).
Figure 5. GF-1 WFV land surface reflectance data ((a) 27 June 2014, (b) 18 July 2014, (c) 15 August 2014, (d) 24 August 2014, and (e) 18 September 2014) and their corresponding LAI estimates ((f): 27 June 2014, (g) 18 July 2014, (h) 15 August 2014, (i) 24 August 2014, and (j) 18 September 2014).
Sensors 17 01593 g005
Sensors 17 01593 g006 550
Figure 6. Scatter plots between field survey LAI and GF-1 WFV data predicated LAI.
Figure 6. Scatter plots between field survey LAI and GF-1 WFV data predicated LAI.
Sensors 17 01593 g006
Table
Table 1. The main characteristics of GF-1 wide field view (WFV) data used in this study.
Table 1. The main characteristics of GF-1 wide field view (WFV) data used in this study.
WFV SensorDate (dd/mm/yy)Matched Field Survey DateData Quality
WFV127 June 201427 June 2014Good
WFV218 July 201421 July 2014With cloud
WFV115 August 201414 August 2014Good
WFV324 August 20145 September 2014Good
WFV418 September 20145 September 2014With cloud
Table
Table 2. The calibration coefficients of GF-1 WFV data in 2014.
Table 2. The calibration coefficients of GF-1 WFV data in 2014.
WFV SensorBandsGainOffset
WFV1Blue band (Band1)0.20040
Green band (Band2)0.16480
Red band (Band3)0.12430
NIR band (Band4)0.15630
WFV2Blue band (Band1)0.17330
Green band (Band2)0.13830
Red band (Band3)0.11220
NIR band (Band4)0.13910
WFV3Blue band (Band1)0.17450
Green band (Band2)0.15140
Red band (Band3)0.12570
NIR band (Band4)0.14620
WFV4Blue band (Band1)0.17130
Green band (Band2)0.16000
Red band (Band3)0.14970
NIR band (Band4)0.14350
Table
Table 3. The input variables of PROSAIL model for LAI estimation algorithm development.
Table 3. The input variables of PROSAIL model for LAI estimation algorithm development.
ParametersUnitsValue RangeStep
LAIm2/m20–70.2
ALA°30–7010
N-1–20.5
Cabμg/cm230–6010
Cmg/cm20.005–0.0150.005
Carμg/cm20-
Cwcm0.005–0.0150.005
Cbrown-0–0.50.5
Hot-0.1-
Solar zenith angle°25–5510

Share and Cite

MDPI and ACS Style

Wei, X.; Gu, X.; Meng, Q.; Yu, T.; Zhou, X.; Wei, Z.; Jia, K.; Wang, C. Leaf Area Index Estimation Using Chinese GF-1 Wide Field View Data in an Agriculture Region. Sensors 2017, 17, 1593. https://doi.org/10.3390/s17071593

AMA Style

Wei X, Gu X, Meng Q, Yu T, Zhou X, Wei Z, Jia K, Wang C. Leaf Area Index Estimation Using Chinese GF-1 Wide Field View Data in an Agriculture Region. Sensors. 2017; 17(7):1593. https://doi.org/10.3390/s17071593

Chicago/Turabian Style

Wei, Xiangqin, Xingfa Gu, Qingyan Meng, Tao Yu, Xiang Zhou, Zheng Wei, Kun Jia, and Chunmei Wang. 2017. "Leaf Area Index Estimation Using Chinese GF-1 Wide Field View Data in an Agriculture Region" Sensors 17, no. 7: 1593. https://doi.org/10.3390/s17071593

APA Style

Wei, X., Gu, X., Meng, Q., Yu, T., Zhou, X., Wei, Z., Jia, K., & Wang, C. (2017). Leaf Area Index Estimation Using Chinese GF-1 Wide Field View Data in an Agriculture Region. Sensors, 17(7), 1593. https://doi.org/10.3390/s17071593

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