Next Article in Journal
Impact of the Revisit of Thermal Infrared Remote Sensing Observations on Evapotranspiration Uncertainty—A Sensitivity Study Using AmeriFlux Data
Next Article in Special Issue
A New CNN-Bayesian Model for Extracting Improved Winter Wheat Spatial Distribution from GF-2 imagery
Previous Article in Journal
Seeking Optimal GNSS Radio Occultation Constellations Using Evolutionary Algorithms
Previous Article in Special Issue
Assimilating Remote Sensing Phenological Information into the WOFOST Model for Rice Growth Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Retrieving Corn Canopy Leaf Area Index from Multitemporal Landsat Imagery and Terrestrial LiDAR Data

1
College of Land Science and Technology, China Agricultural University, No. 17 Qinghua East Road, Haidian District, Beijing 100083, China
2
Key Laboratory of Remote Sensing for Agri-Hazards, Ministry of Agriculture, Beijing 100083, China
3
Department of Geography, The Ohio State University, Columbus, OH 43210, USA
*
Author to whom correspondence should be addressed.
Submission received: 24 January 2019 / Revised: 24 February 2019 / Accepted: 4 March 2019 / Published: 8 March 2019

Abstract

:
Leaf angle is a critical structural parameter for retrieving canopy leaf area index (LAI) using the PROSAIL model. However, the traditional method using default leaf angle distribution in the PROSAIL model does not capture the phenological dynamics of canopy growth. This study presents a LAI retrieval method for corn canopies using PROSAIL model with leaf angle distribution functions referred from terrestrial laser scanning points at four phenological stages during the growing season. Specifically, four inferred maximum-probability leaf angles were used in the Campbell ellipsoid leaf angle distribution function of PROSAIL. A Lookup table (LUT) is generated by running the PROSAIL model with inferred leaf angles, and the cost function is minimized to retrieve LAI. The results show that the leaf angle distribution functions are different for the corn plants at different phenological growing stages, and the incorporation of derived specific corn leaf angle distribution functions distribute the improvement of LAI retrieval using the PROSAIL model. This validation is done using in-situ LAI measurements and MODIS LAI in Baoding City, Hebei Province, China, and compared with the LAI retrieved using default leaf angle distribution function at the same time. The root-mean-square error (RMSE) between the retrieved LAI on 4 September 2014, using the modified PROSAIL model and the in-situ measured LAI was 0.31 m2/m2, with a strong and significant correlation (R2 = 0.82, residual range = 0 to 0.6 m2/m2, p < 0.001). Comparatively, the accuracy of LAI retrieved results using default leaf angle distribution is lower, the RMSE of which is 0.56 with R2 = 0.76 and residual range = 0 to 1.0 m2/m2, p < 0.001. This validation reveals that the introduction of inferred leaf angle distributions from TLS data points can improve the LAI retrieval accuracy using the PROSAIL model. Moreover, the comparations of LAI retrieval results on 10 July, 26 July, 19 August and 4 September with default and inferred corn leaf angle distribution functions are all compared with MODIS LAI products in the whole study area. This validation reveals that improvement exists in a wide spatial range and temporal range. All the comparisons demonstrate the potential of the modified PROSAIL model for retrieving corn canopy LAI from Landsat imagery by inferring leaf orientation from terrestrial laser scanning data.

Graphical Abstract

1. Introduction

Leaf area index (LAI) is an important parameter that controls many physical and biological processes in vegetation canopies and therefore controls the resulting productivity [1]. By definition, LAI is measured as the total one-sided leaf area per unit ground surface area for flat broad leaves; for non-flat leaves; it is half of the total light-intercepting area per unit ground surface area [1]. In precision agriculture, LAI has been used to monitor crop growth and development [2], estimate crop yield [3,4], map and classify crop vitality and yield [5], and to detect early crop stress [6]. Corn (Zea mays) is a productive C4 crop and LAI is a very important indicator for corn yield formation. Thus, the reliable estimation of LAI is an important tool for monitoring corn growth, estimating yield, cultivating corn, and breeding new cultivars.
There are currently two kinds of LAI estimation methods: direct measurements and indirect methods [6]. Compared with direct destructive sampling measurements, indirect estimation methods that use optical instruments and remote sensing images offer the advantage of rapid analysis for large areas, thereby overcoming the drawbacks of labor-intensive and time-consuming “direct” methods. The optical instruments include commercial canopy analysers (i.e., SunSCAN, AccuPAR, LAI-2000/LAI-2200, DEMON, PASTIS-57, digital hemispherical photography, etc.) and smartphone applications (i.e., PocketLAI and LAISmart) [7,8]. These methods are non-destructive, and are based on the statistical and probabilistic approach of foliar element distribution and arrangement within the crop canopy [9].
For the LAI retrieval using remote sensing images, there are two methods: empirical methods and physical models [10]. Empirical methods depend on the location, season, and plant attributes such as species, age, and density; therefore, they do not generalize well over large areas. In contrast, methods based on physical models are more attractive because they are robust over large areas and can be adjusted to account for a wide range of situations.
PROSAIL is a popular physical radiative transfer model for LAI estimation [11]. It was developed to couple the SAIL bidirectional canopy-reflectance model [12] and the PROSPECT model of leaf optical properties [13,14] in the 1990s. Many researchers have estimated canopy biophysical variables using the PROSAIL model in agriculture, forestry, environmental science, and ecology application [11]. However, there are three major problems for retrieving canopy biophysical parameters using PROSAIL: (1) the unknown foliage angle distribution, (2) the error because of the nonrandom spatial distribution of foliage, and (3) the contribution of the supporting non-photosynthetic material to radiation interception [15]. The non-random nature of foliage can be solved by applying simplifications such as a clumping index [6]. This article focuses on the problem of corn leaf angle estimation and identifying non-photosynthetic material of a corn canopy.
Leaf angle is important for estimating the bidirectional reflectance of plant canopies in remote sensing applications [16,17] and PROSAIL output is sensitive to leaf angle input [11]. Despite the importance of measuring the leaf angle distribution, only the Campbell ellipsoidal leaf angle distribution function [18] is commonly used to approximate this distribution by means of empirical expressions that have been parameterized for many species. For many species, including corn, the leaf angle distribution changes throughout the growth process. For corn, new leaves are more erect during the early growth stages but become more horizontal later. Therefore, to accurately retrieve corn LAI, it is necessary to describe the leaf angle distribution throughout the growth process. In the present study, our goal was to develop a model of corn leaf angle distribution suitable for use in the PROSAIL model. Light detection and ranging (LiDAR) technologies have enabled researchers to characterize the 3D canopy structure of corn and other crops, especially based on the high density of the point cloud produced by the LiDAR (Light Detection And Ranging) pulses during terrestrial laser scanning (TLS) [19,20]. Currently, there are two kinds of methods for deriving leaf angle distributions: computing normal vectors of neighboring TLS points [19,21,22] and computing the angle of leaf clusters classified from TLS data [23,24,25]. Zheng and Moskal developed the least square fitting method to compute the normal vectors of artificial tree leaves from TLS points for leaf orientation retrieval [19]. Bailey and Mahaffee calculate the normal vectors of broad-leafed black cottonwood tree and grapevine canopy for leaf angle probability density function estimation [21]. Vicari et al. calculate the normal vectors of leaves TLS points for four broadleaf tree species and estimated their leaf angle distributions [22]. Li et al. clustered the simulated tree leaves using density-based segmentation method for retrieving leaf angle distributions [24]. Xu et al. segmented single-leaf point cloud for extracting leaf angle distributions for three tree species (Ehretia macrophylla, crape myrtle and Fatsia japonica) [25]. These researches all focus on the tree leaf angle distributions using TLS points, especially the broadleaf tree species. Furthermore, this promising new technology of TLS is used to estimate the corn leaf angle distribution in this study, which offers the possibility of improving accuracy of corn canopy LAI retrieval. Therefore, our goal was to retrieve the corn canopy LAI using the inferred leaf angle distributions in a form suitable for incorporation in the PROSAIL model, with data provided by terrestrial LiDAR. Specifically, we aim to (1) develop an approach for inferring the leaf angle distribution functions for corn canopy at four phenological stages, and (2) retrieve corn canopy LAI using the inferred leaf angle distributions obtained from terrestrial LiDAR for these four phenological stages.

2. Materials and Methods

2.1. Study Area

The study area is located in Baoding City, Hebei Province, China, ranging from 39.08°N, 115.48°E to 39.58°N, 116.23°E (Figure 1a). The field is cultivated in a rotation with corn and winter wheat (Triticum aestivum); soybean (Glycine max), cotton (Gossypium hirsutum), and sweet potato (Ipomoea batatas) are grown occasionally. Corn is generally sown at the beginning of June and harvested in the middle of September, that is, a 3-month growing period. The soil is fertile and the summer is hot and rainy. Thus, the study area is the main planted area for corn in North China.

2.2. Data Collection

We measure the corn canopy LAI, chlorophyll content, leaf angles, leaf reflectance, canopy reflectance and background soil reflectance in sampling quadrats (Figure 1b). At each sampling site, we established eight study plots, each 30 m × 30 m in size (to match the resolution of Landsat pixels), with each plot divided into nine 10 m × 10 m quadrats for LAI measurement. Locations were determined using the scanner’s onboard GPS receiver, with a 2D accuracy of 10m, and were used to register the study plots to the corresponding Landsat pixels.

2.2.1. TLS Data Acquisition and Pre-Processing

Laser scanning was conducted with a Focus3D X330 scanner (FARO Technologies, Rugby, U.K., Figure 2a) at the same time as Landsat 7 ETM+ and Landsat 8 OLI image acquisition dates. The TLS scanner employs a 1550-nm laser within a maximum field of view of 300° vertically by 360° horizontally (Table 1, Figure 3b). Its scanning speed is 122,000 to 976,000 points/sec, with a scanning distance of 0.6 m to 330 m. In each quadrat, plot-level TLS scans were performed from three positions outside each plot (Figure 3a). There is one scanning position located inside the plot. Four white balls were placed inside each plot and scanned as registration points; their positions were subsequently used to permit the co-registration of the point clouds generated from different positions because the data obtained from the different scan positions had different coordinate systems. Figure 2 provides photographs of the fieldwork.
In pre-processing, data points representing “noise” associated with real returns from the corn canopy were removed. The noise points are mainly “flying points” for corn plants, because most leaf surfaces are inclined and its normal vectors are not parallel to the laser beam, which will result in the deviation between the received laser pulse and the emitted pulse. The points with this deviation greater than 30 will be classified as “flying points” [26]. The remaining points were co-registered and merged into a single composite point cloud per plot using the FARO SCENE software (www.faro.com). Figure 3c illustrates a sample result of the scanning.
To alleviate the contribution of non-photosynthetic materials to radiation interception, the difference of normals (DoN) method [27,28] is used to identify the points of corn leaves and stalk. Open-source code for the DoN operator was obtained from the Point Cloud Library (PCL; http://pointclouds.org) [29]. This identification is only done on field samplings and only for corn leaf angle distribution function computation.

2.2.2. Acquisition of Remote Sensing Images

Referring to corn’s growing season, we collected two Landsat 7 ETM+ images on 10 July (the stem elongation stage) and 26 July (the heading stage), and two Landsat 8 OLI images on 19 August (the flowering stage) and 4 September (the grain-filling stage) in 2014. All images were downloaded from the United States Geological Survey gateway (https://www.usgs.gov/). The bands used in this study were blue, green, red, and near-infrared (NIR). Table 2 provides details of these four images, the sun elevation and azimuth angle, and the viewing zenith and azimuth angle, which are the inputs of PROSAIL model. We pre-processed the remote sensing images in the whole study area for geo-registration, radiometric calibration, and atmospheric correction to improve the accuracy of the LAI retrieval. The corn planted area is extracted using a decision tree and a mixed-pixel un-mixing classifier in our previous work [30]. The decision tree is built using the NDVI difference calculated from Landsat 8 OLI images on 19 August and 4 September. The pixels with NDVI19 August > 0.6 and NDVI4 September > 0.55 are classified as the corn planted area by testing the NDVI difference of the corn planted area and the non-corn planted area in the whole study area.

2.2.3. Field Data Collection

In order to validate if the inferred leaf angle distributions can be used to improve LAI retrieval, an extensive field campaign was carried out on 4 September, 2014 to measure LAI using a LAI-2200 Plant Canopy Analyzer (LI-COR, Lincoln, NE, USA) (Figure 2c). We used the cover cap with a 45° field angle to eliminate the effect of non-plant objects within the range of the sensor’s field of view. The skylight was measured one time and the light under leaves was measured four times, respectively, for every sampling corn plant. We obtained six measurement values in each quadrat, following a zigzag pattern (Figure 1c). The measured LAI is the effective LAI.
In addition to LAI, we also measured the leaf chlorophyll content using a SPAD-502 leaf chlorophyll meter on 4 September 2014. For this measurement, we used five corn plants and six leaves per plant to measure the chlorophyll content. For each leaf, we measured the chlorophyll content near the leaf tip, leaf bottom, and middle of the leaf, and used the average value to estimate the chlorophyll content of the leaf. The measured relative SPAD values are transformed to absolute chlorophyll concentrations using the exponential fit model proposed by Markwell et al. [31]. The mean of all transformed chlorophyll values for each quadrat represented the actual chlorophyll content.
Finally, we measured the leaf angle by photographing the plants (Figure 2b) and manually analyzing the leaf angles in the images. The leaf reflectance, canopy reflectance and background soil reflectance are measured using SVC HR-1024 spectrometer (Spectra Vista Corporation, Poughkeepsie, NY, USA). The location of all the samplings was provided using a Huace i80 real-time kinematic (RTK) GPS receiver (Huace Ltd., Shanghai, China).

2.3. LAI Retrieval

We use the PROSAIL model to retrieve corn canopy LAI based on the remote sensing images acquired on 10 July, 26 July, 19 August and 4 September, 2014. We focus on testing if the introduction of the leaf angle distribution estimated from TLS data on each date will improve the accuracy of LAI retrieval in this study. Figure 4 summarizes the workflow of LAI retrieval in this study.

2.3.1. PROSAIL Model and Sensitivity Analysis

PROSAIL is the coupling of the PROSPECT model and SAIL model, which simulates the bidirectional reflectance at the top of the canopy in the spectral range from 400 nm to 2500 nm as a function of input variables that relate to the structure of the canopy, the leaf optical properties, the background soil reflectance, and the solar geometry [7,11]. It assumes that the canopy is a turbid medium in which leaves are randomly distributed [32]. The forward simulation expression of PROSAIL model is as followed:
ρ = PROSAIL (N, Cab, Car, Cw, Cm, LIDFa, LAI, hspot, tts, tto, psi, ρsoil)
From this formula, we can see that there are four kinds of input variables for PROSAIL model: leaf optical properties, canopy structure, background soil reflectance and sun-view geometry (see Table 3). Leaf optical properties is described by the mesophyll structural parameter (N), leaf chlorophyll content (Cab), leaf carotenoid content (Cab), dry matter (Cm), and equivalent water thickness (Cw). Leaf chlorophyll (Cab) is measured with a SPAD-502 chlorophyll meter in this study. Equivalent water thickness (Cw) is tied to the difference of fresh leaf weight and dry leaf weight (Cw = (Cfresh leaf − Cdry leaf)/LAI). Canopy structure is characterized by LAI, leaf angle distribution function (LIDFa) and the hot-spot parameter (hspot). LAI comes from the measured LAI value in field work using LAI-2200 Plant Canopy Analyzer. Leaf angle distribution function (LIDFa) comes from the calculated function using TLS points data of corn canopy in this study. The effect background soil (ρsoil) is descripted by the soil reflectance. The sun-view geometry is described by the solar zenith angle (tts), view zenith angle (tto), and the relative azimuth angle between sun and satellite sensor (psi).
In order to constrain the behavior of the PROSAIL model and mitigate the ill-posed problem of LAI retrieval [33], the sensitivity analysis of the model is done firstly during LAI retrieval process. In other words, the sensitivity analysis is done to evaluate the response of canopy reflectance at any wavelength by changing input variables of PROSAIL model. Sensitive input variables can strongly affect model outputs, whereas, insensitive variables have much less impact on the results. The sensitivity analysis is done by running the PROSAIL model using a given range of one input variables within a certain interval (see Table 3), and other input variables are set as unique values. For an example, the PROSAIL model will be run using a range value [0, 8] with 0.01 interval for LAI, and other inputs such as LIDFa, hspot, N, Cab, Car etc. are set as default values when we analyze the sensitivity of LAI input variable. If the simulated reflectance is obviously different using a range of LAI and default values of other variables (i.e., LIDFa, hspot, N, Cab, Car etc.), we can conclude that LAI is a sensitive input parameter.

2.3.2. Leaf Angle Distribution Function Inferred from the TLS Data

An accurate leaf angle distribution function is needed to calculate the radiation flux densities at the leaf surfaces [16,34]. Thus, we focused on improving only the inclination angle density function in this study. One of the most popular mathematical angular distribution models was developed by Campbell [16], who approximated the angular distributions by adjusting the ratio of the semi-long axis (b) to the semi-short axis (a) of a geometrical ellipsoid. The ellipsoidal distribution of leaf angles is represented as follows:
g ( α ) = 2 χ 3 sin α Λ ( cos 2 α + χ 2 sin 2 α ) 2
where α is the leaf inclination angle ( 0 α π / 2 ) ,   χ   = b / a , and
Λ = { χ + ( sin 1 ε ) / ε χ < 1 χ + ln [ ( 1 + ε ) ( 1 ε ) ] 2 ε χ χ > 1 2 χ = 1
where   ε = ( 1 χ 2 ) 1 2 . Zheng and Moskal [19] found that the χ value could be calculated from the plant leaf angle distribution function inferred from discrete terrestrial LiDAR points as defined in Equation (3)
χ 2 = 1 3 sin 2 α + 1
where α is the leaf inclination angle with the highest frequency. Using the TLS point cloud from each of the four scanning dates, we computed the specific χ value for each date by inferring the specific α computed from the normal operator in PCL. Keeping the stalk points will have an effect on the computation of leaf angle distribution functions. So we remove the stalk points from leaf points using DoN method for χ value calculation. The DoN operator is the difference of leaf points normal and stalk points normal. The corn stalk is erect from ground commonly, so the stalk points normal are horizontal. In comparison, the corn leaved will distributed in all directions and their normal are irregular. Therefore, there are distinctly different normal between the normal of corn leaves and stalk. The DoN operator is the difference of leaf points normal and stalk points normal [28], which is used to remove the stalk points from leaf points in this study. Then the computation of corn leaf angle distribution function using corn leaf points is as followed.
Step 1: Computation of centroid for points set from neighboring pi.
For a given TLS point pi, there are k neighboring points which can organized as a point set Pi (Pi = pi1, pi2, pi3, pi4,…, pik). A surface Ti can be constructed using the points set Pi. n i is the normal vector of surface Ti, and p ¯ i is the centroid of points set Pi. And p ¯ i is computed as:
p i = 1 k i = 1 k p i k
Step 2: Computation of normal vector for point pi using PCA method.
Firstly, the covariance matrix M of point set Pi (Pi = pi1, pi2, pi3, pi4,…, pik) is computed using PCA algorithm. Secondly, the eigenvalues (λ1, λ2, λ3) and eigenvectors ( e 1 , e 2 , e 3 ) of covariance matrix M are calculated. Lastly, exploring the eigenvectors of the critical eigenvalues, which is the normal vector for point pi. The covariance matrix M is computed as:
M = 1 k i = 1 k ( p i p ) ( p i p ) T
where, ( p i p ¯ ) is a column vector; ( p i p ¯ ) T is the transpose of column vector; p ¯ is the centroid of neighboring points.
Step 3: Eliminating the directional ambiguity of normal vectors.
There is directional ambiguity randomly for all points: there are two equal but opposite (negative and positive) normals for any object surface, both of them are mathematically valid. Disambiguation of the normals is done by constructing the Riemann surface, and the normal with the minimum spanning tree in the Riemann surface is the direction of the disambiguated normal for the given point.
Step 4: Computation of leaf angle distribution function.
Normal vectors for each point cloud are calculated by step 1 to step 3. Therefore, the corn leaf angle distribution function can be calculated by counting all leaf points normal, and the leaf angle is ranging from 0° to 90°.

2.3.3. LUT-Based LAI Retrieval Strategy Based on PROSAIL Model

1.
LUT generation
LAI retrieval using PROSAIL requires adjusting the values of the input canopy biophysical variables to best match the bidirectional reflectance factors measured by the sensor for a range of directions and wavelength bands [35]. However, retrieving biophysical parameters from PROSAIL is nontrivial. First, the model uses a nonlinear function of the canopy optical and structural parameters. Second, the retrieval solution is not always unique because two or more combinations of canopy parameters may produce similar canopy reflectance spectra [36]. Furthermore, there are uncertainties in the sensor measurements and the input biophysical variables of the canopy model [27]. Therefore, LAI retrieval using PROSAIL is “an ill-posed problem” [35,37,38]. To solve this problem, several algorithms based on minimization have been developed [39], including lookup tables, quasi-Newton algorithms, neural networks, and tangent linear models. Because of it’s ease of use and general robustness, we used the LUT approach in this study.
The determination of LUT dimension is an important issue for LUT generation for LAI retrieval [40]. In order to take balance between the flexibility in determining the optimal band combination for LAI retrieval and mitigation of ill-posed problem, the LUT dimension is designed ranging from blue band to NIR band. Consequently, the spectral wavelength ranging from 450 nm–900 nm is simulated by the PROSAIL model and used to generate LUT, which produces the size of the LUT in column direction. Subsequently, the size of the LUT in row direction should be specified, which is determined by the range/value and interval of input variables [40]. The sun-view geometry (tts, tto, psi) and skyl are coming from the head files of remote images, which are unique values. So, there are only three kinds of parameters in our lookup table: leaf optical properties, canopy structure features and psoil. Weiss et al. [41] found that an LUT based on 100,000 modeled spectra provides an optimal compromise between model accuracy and required computer-resources. Therefore, we combine the input variables using the following range/value and interval of input variables (see Table 3). The version of PROSPECT model used in this study is PROSPECT-5B [42], so the brown pigments (Cbp) is fitted to 0 [43], which is not listed in Table 3.
2.
Cost function
Based on the generated LUT, the program identifies the remote sensing image pixel by pixel. According to the range/value of input variables (see Table 3), there are about 100,000 spectra loaded from the LUT library for LAI retrieval. We compare the pixel values from the blue band to the near infrared band of the Landsat images with the spectrum ranging from 450 nm–900 nm in the lookup table, aiming at finding the best fit(s) for LAI retrieval. The cost function is used to measure the discrepancies between the observed pixel values and the simulated reflectance values in the lookup table. The pair with minimum discrepancy is the best fit/match, which will be stored as the pixel value of LAI output image [40]. Referring to previous studies [35,40], the cost functions based on the root-mean-square error (RMSE) between simulated and observed reflectance values at multiple bands were used for LAI retrieving in this study, which is defined as:
RMSE = 1 n λ = 1 n ( R sim ( λ ) R L ( λ ) ) 2
where RL(λ) is the reflectance in band λ for the Landsat image, Rsim (λ) is the simulated reflectance in band λ, and n is the number of wavelength bands. The retrieved LAI is found when the RMSE approaches 0.

3. Results and Analysis

3.1. Sensitivity Analysis of PROSAIL for Simulating Corn Canopy Reflectance

The sensitivity analysis is done to determine which inputs are sensitive parameters of PROSAIL model and if the leaf angle distribution function (LIDFa) is a sensitive variable [11]. And the sensitive parameters would be set as a ranging value for LUT construction. Figure 5 is the sensitivity analysis result for the 12 input parameters of the PROSAIL model in the blue, green, red and NIR bands. Figure 5a–c,k,l show that the simulated reflectance differs obviously when LAI, LIDFa, N, tts, and tto are changed and other inputs are fixed. This shows that LAI, LIDFa, tts, and tto are sensitive input parameters within blue, green, red and NIR bands. Figure 5f–h show that the simulated reflectance differs within NIR band when Cw, Cm and hspot are changed and other inputs are fixed. This shows that Cw, Cm and hspot are sensitive within NIR band. Figure 5e,i,j show that the simulated reflectance is few different within all bands. This shows that Car, psoil and skyl are insensitive inputs for the PROSAIL model.

3.2. Inferred Leaf Angle Distribution Function from the TLS Scanner Data

We obtained leaf angle distribution functions for the four phenological stages from the TLS scanner data. The χ values simulated by using the Campbell ellipsoid were 1.223 on 10 July, 1.206 on 26 July, 1.214 on 19 August, and 1.195 on 4 September (Table 4). Figure 6 shows the leaf angle distribution functions for the corn canopy modeled as an ellipsoid distribution ( b   >   a ) .
Figure 6 indicates that the leaf angle at the maximum probability density decreased when χ increased through time: the leaf angle intervals around the maximum probability density are from 40° to 50° on 10 July (the stem elongation stage), from 50° to 60° on 26 July (the heading stage) and on 19 August (the flowering stage), and from 60° to 70° on 4 September (the grain-filling stage). Thus, the dominant leaf angle increased over the course of the growing season.

3.3. Retrieved Corn Canopy LAIs

LAI values at four phenological stages were retrieved through using the PROSAIL model based on the lookup table and the leaf angle function inferred from the terrestrial LiDAR data. Figure 7 shows maps of the retrieved LAI derived for the stem elongation stage (10 July), the heading stage (26 July), the flowering stage (19 August), and the grain-ripening stage (4 September). Figure 7a shows low LAI values during the early growing season (10 July, at the seven-leaf stage for corn in study area), with values ranging from 0.5 m2/m2 to 2.7 m2/m2. By 26 July (heading stage; Figure 7b), LAI had increased to 3.7. By late August (flowering stage; Figure 7c), LAI reached values of up to 7.2. By early September (grain filling stage; Figure 7d), LAI reached its seasonal maximum, with values reaching up to 7.8. The retrieved LAI maps therefore show a realistic progression of LAI throughout the growing season.
In order to validate these LAI retrieval results, a further assessment of the accuracy is completed. The extensive in-situ measured LAIs on 4 September are used to determine whether the introduction of inferred leaf angle distributions from TLS data points can improve the LAI retrieval accuracy using the PROSAIL model. In addition, the LAI retrieval results on 10 July, 26 July, 19 August and 4 September with default and inferred corn leaf angle distribution functions are all compared with MODIS LAI products in the whole study area. This comparison is done to evaluate the LAI retrieved results in a wide spatial range and temporal range. Figure 8 shows the relationships between the in-situ measured LAI and the retrieved LAIs using inferred and default leaf angle function for the retrieved results of 4 September. The results show a strong and significant fit for the retrieved LAI using leaf angle function inferred from TLS scanner data (R2 = 0.82, p < 0.001) and the retrieved LAI using default Campbell leaf angle function (R2 = 0.76, p < 0.001). Both retrieved values did not differ significantly from in-situ values (t-test, p < 0.001). Moreover, we found a low RMSE value of 0.31 m2/m2 for the retrieval using leaf angle function inferred from TLS scanner data, compared with the RMSE value of 0.56 m2/m2 for the retrieval using the default Campbell leaf angle function. In addition, the difference between the retrieved LAI using inferred leaf angle function and in-situ value ranged from 0 to 0.6 m2/m2, and the difference between the retrieved using default Campbell leaf angle function and in-situ value ranged from 0 to 1.0 m2/m2. This comparation revealed that the use of leaf angle distribution inferred from TLS data worked better than the default value for this retrieval.
To further evaluate the retrieved LAI results, we compared the LAI retrieving values with inferred and default leaf angle distribution function, Figure 9a is the difference between them for LAI retrieving results on 4 September 2014, i.e., LAIwith inferred LAD-LAIwith default LAD. Figure 9b is the zoomed difference LAI map in the middle of the study area. We can see that the LAI result with inferred leaf angle distribution using TLS data is higher than the LAI result with default leaf angle function. The difference between LAIwith inferred LAD and LAIwith default LAD is ranging from 0 m2/m2 to 1.5 m2/m2. This analysis shows that the use of leaf angle distribution inferred from TLS data increases the LAI retrieved result.
Finally, all LAI retrieved results using default leaf angle and referred leaf angle in the whole study area are compared with homologous MODIS LAI products. Figure 10 shows the compared results of LAI retrieved on 10 July, 26 July, 19 August, 4 September 2014, respectively. Blue boxes represent the retrieved LAI using default corn leaf angle, red boxes represent the retrieved LAI using inferred leaf angle distribution functions from TLS points, and green boxes represent the LAI value from MODIS LAI products. In both panels, solid lines within the boxes represent the mean, box top and bottom represent the 75th and 25th percentile, respectively. The comparison revealed that the LAI variation trend from 10 July to 4 September using inferred corn leaf angle distribution functions is more consistent with the MODIS LAI variation trend than the LAI variation trend using default leaf angle function. Both plots raised quickly from 10 July to 26 July when the corn grew quickly, and both plots raised slowly from 26 July to 4 September when the corn growth was slow. In addition, the curve of retrieved LAI using inferred corn leaf angle distribution functions is higher than the MODIS LAI curve and the curve of retrieved LAI using default leaf angle. Previous comparisons proposed that the MODIS LAI products underestimate the LAIs of corn crops [44] and other vegetation species [10,45,46,47,48]. Fortunately, our experiment shows that this underestimate can be alleviated for corn canopy LAI retrieval by considering the leaf angle difference in different phenological stages.

4. Discussion

Previous research has documented the importance of leaf angle distribution function for LAI retrieval using the PROSAIL model [11,16,18,34]. However, only the Campbell ellipsoidal leaf angle distribution function [18] is commonly used to characterize vegetation leaf angle by defining the ratio of ellipse a and b of a fitted ellipsoid without considering the difference of leaf inclination between different vegetation species and different growing stage, as it is known that most vegetation leaf angle distributions are not precisely ellipsoid. In addition, the leaf inclination angles change during the whole growth period for some crops such as corn. This study proposes the method of estimating the leaf angle distribution functions from terrestrial LiDAR-scanned points at four phenological stages of corn for LAI retrieval using PROSAIL model. Millimeter-scale measurement precision of terrestrial LiDAR points ensures the accuracy of leaf angle distribution function.
Sensitivity analysis results reveal that the LIDFa is a sensitive input from blue band to NIR band for the PROSAIL model. These four bands are the commonly used bands for LAI retrieval. This indicates that the description of the leaf angle distribution function will strongly affect the LAI retrieving accuracy. Therefore, the attempt of inferring leaf angle distribution function from the terrestrial laser scanner points data in this study is meaningful for LAI retrieval based on the PROSAIL model. At the same time, the specification of leaf inclination angle does good to mitigate the ill-posed problem of LAI retrieval [49].
This study has demonstrated that the leaf angle distribution functions are different for the corn plants at different phenological growing stages. Our results show that the maximum corn leaf angle is ranging from 40° to 50° on 10 July (stem elongation stage). The uppermost leaf is almost up-right and the other 6–8 leaves are oblique upward. During the growing period, the corn leaves will flat gradually for intercepting more sunlight. Thus, the maximum corn leaf angle changes from 50° to 60° on 26 July (heading stage) and 19 August (flowering stage) in this study. In the coming growing season, the corn leaf will be curled and the leaf tip will hang down gradually. And the maximum corn leaf angle changes from 60° to 70° on 4 September (grain-filling stage) for our studied corn plants. This leaf angle distribution function changing is in line with the corn leaf growing rhythm. This result confirms the assumption that the leaf angle distribution of corn canopy is clearly changing during different phenological phases over the growing season.
Our LAI retrieved results in four corn phenological stages show that the incorporation of derived specific corn leaf angle distribution functions distributes the improvement of LAI retrieval using the PROSAIL model. Two accuracy assessments are completed. One is the comparison of the retrieved LAIs using inferred and default leaf angle distribution functions with the measured LAIs on the in-situ measuring quadrats on 4 September, 2014. The comparison revealed that the accuracy of LAI retrieval using the inferred leaf angle function from TLS scanner data is higher than using the default Campbell leaf angle function. This improvement is a result of the fact our method takes into account the phenological change of corn leaf angle distributions using TLS points echoing on corn leaves. Of course, there are relevant input variables are dependent on the phenological stage and influence LAI retrieval next to leaf angle, such as leaf carotenoid content (Cab), dry matter (Cm) etc. We are focusing on depicting the leaf angle change in this study. And we will expand the LAI retrieval improvement by bringing in other input variables in future work. The other validation is done in a wide spatial range (the whole study area) and multi-temporal (10 July, 26 July, 19 August and 4 September, 2014). Comparison results show that the spatial pattern and temporal trajectory of LAI retrieved results using inferred corn leaf angle distributions are closer to MODIS LAI products. This conclusion reveals that the LAI retrieval improvements can exist in a time series LAI retrieval and generalize to wider corn planted areas. Much research has indicated that the retrieved LAI using a remote sensing technique is underestimated [50,51,52,53,54], especially for MODIS LAI [10,44,45,46,47,48]. Therefore, the inferred leaf angle distribution improves the underestimation of LAI retrieval using remote sensing images.
This study demonstrates the ability of derived corn leaf angle distribution function to improve LAI retrieval accuracy using PROSAL model. This method leads to the improvement of LAI retrieval using a physical model. However, some limitations are worth noting. One is that we only depict the leaf angle distribution function for corn canopy, and there is no other crop studied currently. We will extend this method to other crops, such as sugarcane, wheat, rice, potato etc. The other limitation is that the terrestrial LiDAR scanner is costly. Namely, the LiDAR observations are valuable and often limited in availability. Therefore, we will explore the possibility of deriving the leaf angle distribution functions of different phenological dates by taking stereoscopic photographs using a photographic measurement method in future work.

5. Conclusions

In this study, we propose a method for retrieving multitemporal corn canopy LAI values by using a TLS scanner data to infer an empirical leaf angle probability density function suitable for use in the PROSAIL model. The novelty of this study is that different empirical leaf angle distribution functions for the corn canopy were estimated from terrestrial LiDAR data at four phenological stages. The results show that the inferred leaf angle function improved the accuracy of LAI retrieval compared with the default Campbell leaf angle function. In addition, we localized the PROSAIL model to retrieve corn canopy LAI by accounting for its sensitivity to changes in parameters such as the leaf angle distribution function, LAI, and leaf chlorophyll (a+b) content.
The retrieved LAI results demonstrate the potential of retrieving multitemporal corn canopy LAI by inferring specific leaf angle distribution functions for use in the PROSAIL model and the ability of this method to account for phenological changes during the growing season. By adopting this method, it should be possible to improve the monitoring of crop growth, as well as the detection of stresses such as drought and yield prediction. Future studies can be focus on the following directions. First, exploring how to infer other sensitive inputs to improve vegetation canopy parameter retrieval, and ultimately resolve the ill-posed problem. Second, harmonizing the prior knowledge of crop growth to the canopy parameter retrieval, for example, LAI change, chlorophyll content change etc. in the whole growing season. (3) Developing different cost functions to account for model effects, including 4D-VAR method, which considers observation errors and model errors to determine the best fits between observed and simulated reflectance signatures.

Author Contributions

This work is cooperated by our research team, and the contributions are as followed. Conceptualization, W.S. and J.H.; methodology, W.S. and M.Z.; software, M.Z.; writing—original draft preparation, W.S.; writing—review and editing, D.L.

Funding

This study was funded by the National Natural Science Foundation of China under the project ‘Growth process monitoring of corn by combining time series spectral remote sensing images and terrestrial laser scanning data’ (41671433), ‘Estimating the leaf area index of corn in whole growth period using terrestrial LiDAR data’ (41371327), Science and Technology Facilities Council of UK- Newton Agritech Programme (Project No. ST/N006798/1) and the National Keypoint Research and Invention Program of the Thirteenth 5-year Plan (2017YFD0300903).

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. Huang, J.X.; Tian, L.Y.; Liang, S.L.; Ma, H.Y.; Becker-Reshef, I.; Su, W.; Huang, Y.B.; Zhang, X.D.; Zhu, D.H.; Wu, W.B. 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]
  3. Huang, J.; Ma, H.; Su, W.; Zhang, X.; Huang, Y.; Fan, J.; Wu, W.B. Jointly assimilating MODIS LAI and ET products into the SWAP model for winter wheat yield estimation. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 4060–4071. [Google Scholar] [CrossRef]
  4. Huang, J.X.; Sedano, F.; Huang, Y.B.; Ma, H.Y.; Li, X.L.; Liang, S.L.; Tian, L.Y.; Zhang, X.D.; Fan, J.L.; Wu, W.B. Assimilating a synthetic Kalman filter leaf area index series into the WOFOST model to estimate regional winter wheat yield. Agric. For. Meteorol. 2016, 216, 188–202. [Google Scholar] [CrossRef]
  5. Rudolph, S.; Kruk, J.V.D.; Hebel, C.V.; Ali, M.; Herbst, M.; Montzka, C.; Pätzold, S.; Robinson, D.A.; Vereecken, H.; Weihermüller, L. Linking satellite derived LAI patterns with subsoil heterogeneity using large-scale ground-based electromagnetic induction measurements. Geoderma 2015, s241–s242, 262–271. [Google Scholar] [CrossRef]
  6. Zheng, G.; Moskal, L.M.; Kim, S.H. Retrieval of effective leaf area index in heterogeneous forests with terrestrial laser scanning. IEEE Trans. Geosci. Remote Sens. 2013, 51, 777–786. [Google Scholar] [CrossRef]
  7. Campos-Taberner, M.; García-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]
  8. Fang, H.L.; Ye, Y.C.; Liu, W.L.; Wei, S.S.; Ma, L. Continuous estimation of canopy leaf area index (LAI) and clumping index over broadleaf crop fields: An investigation of the PASTIS-57 instrument and smartphone applications. Agric. For. Meteorol. 2018, 253, 48–61. [Google Scholar] [CrossRef]
  9. Bréda, N.J.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] [PubMed]
  10. 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]
  11. Jacquemoud, S.; Verhoef, W.; Baret, F.; Bacour, C.; Zarcotejada, P.J.; Asner, G.P.; François, C.; Ustin, S.L. PROSPECT+SAIL models: A review of use for vegetation characterization. Remote Sens. Environ. 2009, 113, S56–S66. [Google Scholar] [CrossRef]
  12. 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]
  13. Jacquemoud, S.; Baret, F. PROSPECT: A model of leaf optical properties spectra. Remote Sens. Environ. 1990, 34, 75–91. [Google Scholar] [CrossRef]
  14. Baret, F.; Jacquemoud, S.; Guyot, G.; Leprieur, C. Modeled analysis of the biophysical nature of spectral shifts and comparison with information content of broad bands. Remote Sens. Environ. 1992, 41, 133–142. [Google Scholar] [CrossRef]
  15. Chen, J.M.; Cihlar, J. Quantifying the effect of canopy architecture on optical measurements of leaf area index using two gap size analysis methods. IEEE Trans. Geosci. Remote Sens. 1995, 33, 777–787. [Google Scholar] [CrossRef]
  16. Campbell, G.S. Derivation of an angle density function for canopies with ellipsoidal leaf angle distributions. Agric. For. Meteorol. 1990, 49, 173–176. [Google Scholar] [CrossRef]
  17. Zou, X.; Mõttus, M. Sensitivity of Common Vegetation Indices to the Canopy Structure of Field Crops. Remote Sens. 2017, 9, 994. [Google Scholar] [CrossRef]
  18. Campbell, G.S. Extinction coefficients for radiation in plant canopies calculated using an ellipsoidal inclination angle distribution. Agric. For. Meteorol. 1986, 36, 317–321. [Google Scholar] [CrossRef]
  19. Zheng, G.; Moskal, L.M. Leaf orientation retrieval from terrestrial laser scanning (TLS) data. IEEE Trans. Geosci. Remote Sens. 2012, 50, 3970–3979. [Google Scholar] [CrossRef]
  20. Greaves, H.E.; Vierling, L.A.; Eitel, J.U.H.; Boelman, N.T.; Magney, T.S.; Prager, C.M.; Griffin, K.L. Estimating aboveground biomass and leaf area of low-stature arctic shrubs with terrestrial lidar. Remote Sens. Environ. 2015, 164, 26–35. [Google Scholar] [CrossRef]
  21. Bailey, B.N.; Mahaffee, W.F. Rapid measurement of the three-dimensional distribution of leaf orientation and the leaf angle probability density function using terrestrial LiDAR scanning. Remote Sens. Environ. 2017, 194, 63–76. [Google Scholar] [CrossRef]
  22. Vicari, M.B.; Pisek, J.; Disney, M. New estimates of leaf angle distribution from terrestrial LiDAR Comparison with measured and modelled estimates from nine broadleaf tree species. Agric. For. Meteorol. 2019, 264, 322–333. [Google Scholar] [CrossRef]
  23. Hosoi, F.; Nakai, Y.; Omasa, K. Estimating the leaf inclination angle distribution of the wheat canopy using a portable scanning lidar. J. Agric. Meteorol. 2009, 65, 297–302. [Google Scholar] [CrossRef] [Green Version]
  24. Li, Y.M.; Su, Y.J.; Hu, T.Y.; Xu, G.C.; Guo, Q.H. Retrieving 2-D leaf angle distributions for deciduous trees from terrestrial laser scanner data. IEEE Trans. Geosci. Remote Sens. 2018, 56, 4945–4955. [Google Scholar] [CrossRef]
  25. Xu, Q.F.; Cao, L.; Xue, L.F.; Chen, B.C.; An, F.; Yun, T. Extraction of leaf biophysical attributes based on a computer graphic-based algorithm using terrestrial laser scanning data. Remote Sens. 2019, 11, 15. [Google Scholar] [CrossRef]
  26. Li, D.; Pang, Y.; Yue, C.R.; Zhao, D.; Xu, G.C. Extraction of individual tree DBH and height based on terrestrial laser scanner data. J. Beijing For. Univ. 2012, 34, 79–86. [Google Scholar]
  27. Su, W.; Zhu, D.H.; Huang, J.X.; Guo, H. Estimation of the vertical leaf area profile of corn (Zea mays) plants using terrestrial laser scanning (TLS). Comput. Electron. Agric. 2018, 150, 5–13. [Google Scholar] [CrossRef]
  28. Su, W.; Zhang, M.Z.; Liu, J.M.; Sun, Z.P. Automated extraction of corn leaf points from unorganized terrestrial lidar point clouds. Int. J. Agric. Biol. Eng. 2018, 11, 166–170. [Google Scholar] [CrossRef]
  29. Rusu, R.B.; Cousins, S. 3D is here: Point Cloud Library (PCL). In Proceedings of the IEEE International Conference on Robotics and Automation (ICRA), Shanghai, China, 9–13 May 2011; pp. 1–4. [Google Scholar]
  30. Su, W.; Jiang, F.F.; Zhu, D.H.; Zhan, J.G.; Ma, H.Y.; Zhang, X.D. Extraction of maize planting area based on decision tree and mixed-pixel unmixing methods. Trans. Chin. Soc. Agric. Mach. 2015, 46, 289–295. [Google Scholar]
  31. Markwell, J.; Osterman, J.C.; Mitchell, J.L. Calibration of the minolta spad-502 leaf chlorophyll meter. Photosynth. Res. 1995, 46, 467–472. [Google Scholar] [CrossRef]
  32. 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]
  33. Quan, X.; He, B.; Li, X.; Liao, Z. Retrieval of grassland live fuel moisture content by parameterizing radiative transfer model with interval estimated LAI. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 9, 910–920. [Google Scholar] [CrossRef]
  34. Propastin, P.; Panferov, O. Retrieval of remotely sensed LAI using landsat ETM+ data and ground measurements of solar radiation and vegetation structure: Implication of leaf inclination angle. Int. J. Appl. Earth Obs. Geoinf. 2013, 25, 38–46. [Google Scholar] [CrossRef]
  35. Combal, B.; Baret, F.; Weiss, M.; Trubuil, A.; Macé, D.; Pragnère, 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]
  36. Ding, Y.L.; Zhang, H.Y.; Li, Z.W.; Xin, X.P.; Zheng, X.M.; Zhao, K. Comparison of fractional vegetation cover estimations using dimidiate pixel models and look-up table inversions of the PROSAIL model from Landsat 8 OLI data. J. Appl. Remote Sens. 2016, 10, 036022. [Google Scholar] [CrossRef]
  37. Duan, S.B.; Li, Z.L.; Wu, H.; Tang, B.H.; Ma, L.; Zhao, E.Y.; Li, C.R. 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]
  38. Zurita-Milla, R.; Laurent, V.C.E.; Gijsel, J.A.E.V. Visualizing the ill-posedness of the inversion of a canopy radiative transfer model: A case study for sentinel-2. Int. J. Appl. Earth Obs. Geoinf. 2015, 43, 7–18. [Google Scholar] [CrossRef]
  39. Chaabouni, S.; Kallel, A. Inversion of vegetation canopy reflectance based on variationnel multiscale approach. In Proceedings of the 2nd International Conference on Advanced Technologies for Signal and Image Processing—ATSIP, Monastir, Tunisia, 21–24 March 2016; pp. 494–498. [Google Scholar]
  40. Locherer, M.; Hank, T.; Danner, M.; Mauser, W. Retrieval of seasonal leaf area index from simulated EnMAP data through optimized LUT-based inversion of the PROSAIL model. Remote Sens. 2015, 7, 10321–10346. [Google Scholar] [CrossRef]
  41. 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]
  42. Féret, J.B.; Gitelson, A.A.; Noble, S.D.S. Jacquemoud. PROSPECT-D: Towards modeling leaf optical properties through a complete lifecycle. Remote Sens. Environ. 2017, 193, 204–215. [Google Scholar] [CrossRef]
  43. García-Haro, F.J.; Campos-Taberner, M.; Muñoz-Marí, J.; Laparra, V.; de Coca, F.C.; Sánchez-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]
  44. Yang, F.; Sun, J.L.; Zhang, B.; Yao, Z.F.; Wang, Z.M.; Wang, J.L.; Yue, X.F. Assessment of MODIS LAI product accuracy based on the PROSAIL model, TM and field measurements. Trans. Chin. Soc. Agric. Eng. 2010, 26, 192–197. [Google Scholar]
  45. Yuan, H.; Dai, Y.; Xiao, Z.; Ji, D.; Shangguan, W. Reprocessing the MODIS leaf area index products for land surface and climate modelling. Remote Sens. Environ. 2011, 115, 1171–1187. [Google Scholar] [CrossRef]
  46. Fang, H.; Jiang, C.; Li, W.; Wei, S.; Baret, F.; Chen, J.M.; Garcia-Haro, J.; Liang, S.; Liu, R.; Myneni, R.B.; et al. Characterization and intercomparison of global moderate resolution leaf area index (LAI) products: Analysis of climatologies and theoretical uncertainties. J. Geophys. Res. Biogeosci. 2013, 118, 529–548. [Google Scholar] [CrossRef] [Green Version]
  47. Li, Z.W.; Tang, H.; Zhang, B.H.; Yang, G.X.; Xin, X.P. Evaluation and intercomparison of MODIS and GEOV1 global leaf area index products over four sites in north China. Sensors 2015, 15, 6196–6216. [Google Scholar] [CrossRef]
  48. Xu, X.J.; Du, H.Q.; Zhou, G.M.; Li, P.H. Method for improvement of MODIS leaf area index products based on pixel-to-pixel correlations. Eur. J. Remote Sens. 2016, 49, 57–72. [Google Scholar] [CrossRef] [Green Version]
  49. Justice, C.; Belward, A.; Morisette, J.; Lewis, P.; Privette, J.; Baret, F. Developments in the ‘validation’ of satellite sensor products for the study of the land surface. Int. J. Remote Sens. 2000, 21, 3383–3390. [Google Scholar] [CrossRef]
  50. Morisette, J.T.; Baret, F.; Privette, J.L.; Myneni, R.B.; Nickeson, J.E.; Garrigues, S.; Shabanov, N.V.; Weiss, M.; Fernandes, R.A.; Leblanc, S.G.; et al. Validation of global moderate-resolution LAI products: A framework proposed within the CEOS land product validation subgroup. IEEE Trans. Geosci. Remote Sens. 2006, 44, 1804–1817. [Google Scholar] [CrossRef]
  51. Weiss, M.; Baret, F.; Garrigues, S.; Lacaze, R. LAI and fAPAR cyclopes global products derived from vegetation. Part 2: Validation and comparison with MODIS collection 4 products. Remote Sens. Environ. 2007, 110, 317–331. [Google Scholar] [CrossRef]
  52. Fang, H.; Wei, S.; Jiang, C.; Scipal, K. Theoretical uncertainty analysis of global MODIS, CYCLOPES and GLOBCARBON LAI products using a triple collocation method. Remote Sens. Environ. 2012, 124, 610–621. [Google Scholar] [CrossRef]
  53. Zou, X.; Haikarainen, I.; Haikarainen, I.P.; Mäkelä, P.; Mõttus, M.; Pellikka, P. Effects of crop leaf angle on LAI-sensitive narrow-band vegetation indices derived from imaging spectroscopy. Appl. Sci.-Basel 2018, 8, 1435. [Google Scholar] [CrossRef]
  54. Huang, J.X.; Ma, H.Y.; Sedano, F.; Lewis, P.; Liang, S.L.; Wu, Q.L.; Su, W.; Zhang, X.D.; Zhu, D.H. Evaluation of regional estimates of winter wheat yield by assimilating three remotely sensed reflectance datasets into the coupled WOFOST–PROSAIL model. Eur. J. Agron. 2019, 102, 1–13. [Google Scholar] [CrossRef]
Figure 1. Location of the study area (a), the distribution of sampling quadrats in field work (b), and sampling pattern used in the fieldwork (c).
Figure 1. Location of the study area (a), the distribution of sampling quadrats in field work (b), and sampling pattern used in the fieldwork (c).
Remotesensing 11 00572 g001
Figure 2. In-situ data collection activities conducted in each sampling plot: (a) obtaining a 3D image of the plot using the TLS scanner, (b) photographing the leaves for manual calculation of leaf angles, and (c,d) measuring LAI using the LI-COR LAI-2200C plant canopy analyzer.
Figure 2. In-situ data collection activities conducted in each sampling plot: (a) obtaining a 3D image of the plot using the TLS scanner, (b) photographing the leaves for manual calculation of leaf angles, and (c,d) measuring LAI using the LI-COR LAI-2200C plant canopy analyzer.
Remotesensing 11 00572 g002
Figure 3. Top view (a) and side view (b) of the laser scanning scheme for the corn canopy, and (c) sample image of the scanning data for a corn canopy.
Figure 3. Top view (a) and side view (b) of the laser scanning scheme for the corn canopy, and (c) sample image of the scanning data for a corn canopy.
Remotesensing 11 00572 g003
Figure 4. Overview of the LAI retrieval based on improving the leaf angle distribution function of the PROSAIL model.
Figure 4. Overview of the LAI retrieval based on improving the leaf angle distribution function of the PROSAIL model.
Remotesensing 11 00572 g004
Figure 5. Results of the sensitivity analysis for the 12 input parameters to which PROSAIL is sensitive.
Figure 5. Results of the sensitivity analysis for the 12 input parameters to which PROSAIL is sensitive.
Remotesensing 11 00572 g005
Figure 6. The leaf angle probability density functions for the corn canopy on 10 July, 26 July, 19 August, and 4 September, 2014, respectively.
Figure 6. The leaf angle probability density functions for the corn canopy on 10 July, 26 July, 19 August, and 4 September, 2014, respectively.
Remotesensing 11 00572 g006
Figure 7. The LAI maps derived using the inferred leaf angle function from TLS scanner data on (a) 10 July (stem elongation stage), (b) 26 July (heading stage), (c) 19 August (flowering state), and (d) 4 September (grain-filling stage).
Figure 7. The LAI maps derived using the inferred leaf angle function from TLS scanner data on (a) 10 July (stem elongation stage), (b) 26 July (heading stage), (c) 19 August (flowering state), and (d) 4 September (grain-filling stage).
Remotesensing 11 00572 g007
Figure 8. The relationships between the in-situ measured LAI and the retrieved LAIs using inferred and default leaf angle function. Values represent means ± SD based on the in-situ LAI measurements at each sample plot.
Figure 8. The relationships between the in-situ measured LAI and the retrieved LAIs using inferred and default leaf angle function. Values represent means ± SD based on the in-situ LAI measurements at each sample plot.
Remotesensing 11 00572 g008
Figure 9. LAI bias (m2/m2) between the retrieved LAI with and without inferred leaf angle distribution function on 4 September 2014.
Figure 9. LAI bias (m2/m2) between the retrieved LAI with and without inferred leaf angle distribution function on 4 September 2014.
Remotesensing 11 00572 g009
Figure 10. Comparation between the retrieved LAIs with default and inferred leaf angle distribution function on 10 July, 26 July, 19 August, 4 September, 2014, respectively.
Figure 10. Comparation between the retrieved LAIs with default and inferred leaf angle distribution function on 10 July, 26 July, 19 August, 4 September, 2014, respectively.
Remotesensing 11 00572 g010
Table 1. The parameters of the FARO Focus3D X330 LiDAR surveying instrument.
Table 1. The parameters of the FARO Focus3D X330 LiDAR surveying instrument.
ParameterRange of Values
Scanning distance (m)0.6 to 330
Scanning speed (points/s)122,000 to 976,000
Ranging error (mm)±2
Resolution (pixels)7 × 107
Vertical field of view (°)300
Horizontal field of view (°)360
Laser classClass 1
Wavelength (nm)1550
GPSIntegrated GPS receiver
Table 2. The parameters of the Landsat 7 and 8 images.
Table 2. The parameters of the Landsat 7 and 8 images.
DateSensorUTM TimeSun Elevation Angle (°)Sun Azimuth Angle (°)Viewing Zenith Angle (°)Viewing Azimuth Angle (°)
10 JulyETM+02:51:2764.77124.82090
26 JulyETM+02:51:3162.52128.55090
19 AugustOLI02:53:5958.03138.33090
4 SeptemberOLI02:54:0253.72144.98090
Table 3. Range and distribution of the input variables used to establish the synthetic corn canopy reflectance database in the lookup table.
Table 3. Range and distribution of the input variables used to establish the synthetic corn canopy reflectance database in the lookup table.
Model VariablesRange or ValueDistribution
CanopyLAILeaf area index (m2 m−2)0.1 to 7.0Uniform
LIDFaLeaf angle distribution (º)0 to 90Gaussian
hspotHotspot parameter (m m−1)0.1-
LeafNLeaf structural parameter in PROSPECT1.518-
CabChlorophyll a+b content in PROSPECT (μg cm−2)0.1 to 60.0Uniform
CarCarotenoid content in PROSPECT (μg cm−2)8-
CwEquivalent water thickness in PROSPECT (cm)0.05 to 0.3Gaussian
CmDry matter content in PROSPECT (g cm−2)0.002 to 0.012Gaussian
Soil and skypsoilSoil reflectance assumed to be Lambertian (1) or not (0)0–1Gaussian
skylRatio of diffuse to total incident radiationCalculated by tts-
Sun-sensorttsSolar zenith angle (°)//
ttoViewing zenith angle (°)//
psiRelative azimuth angle (v)//
The symbol “/” represents values obtained from the Landsat data header, “-” represents the one set value, and skyl is calculated by tts using skyl = 0.847 − 1.61 × sin((90 − tts) × rd) + 1.04 × sin((90 − tts) × rd) × sin((90 − tts) × rd). Notes: The boundary conditions for the corn canopy, leaves, and soil were selected to describe the characteristics of all growth conditions for the corn canopy in our study area. Sun and sensor viewing conditions corresponded to the measurement conditions when the satellite passed overhead. All possible combinations of all variables were calculated from the corn leaf and canopy input ranges.
Table 4. Proportion of the leaf angles (in 10° intervals) on the four phenological dates.
Table 4. Proportion of the leaf angles (in 10° intervals) on the four phenological dates.
DateχProportion of Leaf Angles (% of Total)
0°–10°10°–20°20°–30°30°–40°40°–50°50°–60°60°–70°70°–80°80°–90°
10 July1.2237.8311.3314.216.116.6815.0510.635.682.5
26 July1.2066.498.9810.7312.4214.1815.4114.6711.16.02
19 August1.2147.110.1111.311.8913.414.7914.0210.856.54
4 September1.1956.178.7610.0210.7312.1914.5415.7913.538.27
Notes: Boldfaced values represent the interval with the highest proportion of the total.

Share and Cite

MDPI and ACS Style

Su, W.; Huang, J.; Liu, D.; Zhang, M. Retrieving Corn Canopy Leaf Area Index from Multitemporal Landsat Imagery and Terrestrial LiDAR Data. Remote Sens. 2019, 11, 572. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11050572

AMA Style

Su W, Huang J, Liu D, Zhang M. Retrieving Corn Canopy Leaf Area Index from Multitemporal Landsat Imagery and Terrestrial LiDAR Data. Remote Sensing. 2019; 11(5):572. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11050572

Chicago/Turabian Style

Su, Wei, Jianxi Huang, Desheng Liu, and Mingzheng Zhang. 2019. "Retrieving Corn Canopy Leaf Area Index from Multitemporal Landsat Imagery and Terrestrial LiDAR Data" Remote Sensing 11, no. 5: 572. https://0-doi-org.brum.beds.ac.uk/10.3390/rs11050572

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