Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
SSCNet: A Spectrum-Space Collaborative Network for Semantic Segmentation of Remote Sensing Images
Previous Article in Journal
Dual-Channel Semi-Supervised Adversarial Network for Building Segmentation from UAV-Captured Images
Previous Article in Special Issue
AM–GM Algorithm for Evaluating, Analyzing, and Correcting the Spatial Scaling Bias of the Leaf Area Index
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Enhancing Leaf Area Index Estimation with MODIS BRDF Data by Optimizing Directional Observations and Integrating PROSAIL and Ross–Li Models

1
School of Geographic and Environmental Sciences, Tianjin Normal University, Tianjin 300387, China
2
School of Mechatronical Engineering, Beijing Institute of Technology, Beijing 100081, China
3
Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2023, 15(23), 5609; https://doi.org/10.3390/rs15235609
Submission received: 6 October 2023 / Revised: 27 November 2023 / Accepted: 29 November 2023 / Published: 2 December 2023

Abstract

:
The Leaf Area Index (LAI) is a crucial vegetation parameter for climate and ecological models. Reflectance anisotropy contains valuable supplementary information for the retrieval of properties of an observed target surface. Previous studies have utilized multi-angular reflectance data and physically based Bidirectional Reflectance Distribution Function (BRDF) models with detailed vegetation structure descriptions for LAI estimation. However, the optimal selection of viewing angles for improved inversion results has received limited attention. By optimizing directional observations and integrating the PROSAIL and Ross–Li models, this study aims to enhance LAI estimation from MODIS BRDF data. A dataset of 20,000 vegetation parameter combinations was utilized to identify the directions in which the PROSAIL model exhibits higher sensitivity to LAI changes and better consistency with the Ross–Li BRDF models. The results reveal significant variations in the sensitivity of the PROSAIL model to LAI changes and its consistency with the Ross–Li model over the viewing hemisphere. In the red band, directions with high sensitivity to LAI changes and strong model consistency are mainly found at smaller solar and viewing zenith angles. In the near-infrared band, these directions are distributed at positions with larger solar and viewing zenith angles. Validation using field measurements and LAI maps demonstrates that the proposed method achieves comparable accuracy to an algorithm utilizing 397 viewing angles by utilizing reflectance data from only 30 directions. Moreover, there is a significant improvement in computational efficiency. The accuracy of LAI estimation obtained from simulated multi-angle data is relatively high for LAI values below 3.5 when compared with the MODIS LAI product from two tiles. Additionally, there is also a slight improvement in the results when the LAI exceeds 4.5. Overall, our results highlight the potential of utilizing multi-angular reflectance in specific directions for vegetation parameter inversion, showcasing the promise of this method for large-scale LAI estimation.

1. Introduction

The leaf area index (LAI), occupying a crucial role in the vegetation canopy [1], is defined as the hemi-surface area of all leaves or needles in the vegetation canopy divided by the horizontal ground surface area [2,3]. LAI serves as a key indicator of vegetation health, productivity, and overall ecosystem functioning [4]. It provides essential information about the size and distribution of the leaf canopy, relating to various processes such as photosynthesis [5], evaporation and transpiration [6,7], and carbon assimilation [8]. Accurately estimating the LAI is of utmost importance as it has been employed in estimating biophysical parameters of vegetation [9,10] and in monitoring vegetation growth [11], as well as in the study of global climate change [12], ecosystem productivity, biogeochemistry, hydrology, and ecology [13]. LAI also plays a crucial role in modeling the exchanges of mass, energy, and momentum between the biosphere and the atmosphere [14,15]. Therefore, ensuring the accuracy of LAI estimation is pivotal for effective ecological research, environmental monitoring, and sustainable land management strategies.
The measurement of the LAI can be achieved through various approaches [16]. Direct methods entail partially or completely defoliating the canopy to assess the total leaf area of plants or trees, which involves destructiveness and a time-consuming process [17,18]. Indirect approaches employ mathematical algorithms to characterize the passage of light through the canopy, utilizing Beer’s Law for estimating the total leaf area [16,19]. In recent years, technological advancements have revolutionized the way we measure leaf area. Surface reflectance data captured by moderate-resolution sensors have been employed in the generation of several global LAI products [20,21,22,23,24], and these products have been extensively validated using field measurements and upscaled high-resolution LAI reference maps [25,26,27,28,29]. These technologies offer the potential to enhance our understanding of vegetation dynamics over vast geographical areas, providing valuable information for global climate change studies. However, accurate measurements face challenges posed by factors like cloud cover, satellite orbit constraints, and the requirement for precise multi-angle data [30].
Ground-based experiments have demonstrated that utilizing abundant multi-angle information can significantly improve LAI estimation accuracy [31]. This improvement is attributed to the fact that multi-angle data provide a wealth of additional information about vegetation parameters. Satellite-observed reflectance at multiple viewing angles, along with physical BRDF models that provide detailed descriptions of vegetation structure, are often directly utilized for the retrieval of vegetation structure parameters [32]. Physical BRDF models, based on physical laws, describe how radiation interacts with vegetation canopies. These models enable the calculation of reflectance under different viewing and illumination conditions, taking into account a wide range of leaf and canopy structure parameters. By employing different inversion techniques, physical BRDF models can estimate these vegetation parameters using reflectance data as input [33]. For example, using the PROSAIL model can simulate the interaction of light with vegetation at different scales, from individual leaves to entire canopies, and provides outputs such as leaf reflectance, transmittance, and absorptance, as well as canopy reflectance [34,35]. By employing the PROSAIL model in inverse mode, it is possible to obtain biophysical parameters such as LAI, chlorophyll content, and vegetation water content from canopy reflectance [36].
Although multi-angular measurements provide enhanced information regarding the structural properties of vegetation [37,38,39,40], it is important to acknowledge that the capability of satellites to acquire multi-angle data is constrained by the characteristics of their orbits as well as by the presence of clouds and shadows [41], leading to a loss of anisotropic information. The number of angle samples in remote sensing observations often fails to meet the sample quantity requirements of the physical BRDF model inversion. Additionally, reflectance measurements with similar geometric configurations tend to exhibit significant autocorrelation, while random noise in these measurements is closely associated with vegetation structure [42]. Moreover, during the inverse process, iterative optimization techniques may encounter the challenge of converging to a local minimum instead of the global minimum [31], which can result in potentially inaccurate outcomes. All the above-mentioned factors may have an impact on the accuracy of LAI inversion.
By accumulating observations over time from satellite sensors with a wide field of view, such as those onboard the Moderate Resolution Imaging Spectroradiometer (MODIS) [43] and the Multi-angle Imaging SpectroRadiometer (MISR) [44], along with Ross–Li models, BRDF products have been successfully generated at both global and regional scales. The surface BRDF product provides a means to characterize the directional properties of the underlying surface reflection across the entire observation hemisphere [45,46]. Coupling different models is a commonly employed strategy to utilize their individual strengths and enhance the accuracy of modeling surface reflection and the estimation of vegetation parameters. A model-to-model approach has been proposed to avoid most of the challenges associated with inverting LAI, which by using multi-angular reflectance can be resolved [47]. In this kind of approach, the multi-angular measurements are used to fit a simplified BRDF model that requires only a few parameters to describe the surface reflectance anisotropy. For example, by fitting the Ross–Li BRDF model to the multi-angular reflectance data, the model provides estimates of the BRDF parameters, which can quantitatively describe how surface reflectance varies with different viewing and illumination angles [42,48]. Parameters of the model are then estimated to facilitate the inversion process. The BRDF information reconstructed from these satellite observations has been successfully applied in numerous parameter estimations, such as albedo, clumping index and canopy height [43,49,50,51,52,53], laying the foundation for improving LAI satellite inversion algorithms. Although the model-to-model approach provides a solution for using multi-angle observation to inverse LAI, it may introduce additional sources of error or uncertainty due to the additional modeling step.
Zhang et al. [30] coupled the PROSAIL model with the hotspot-corrected RossThick–LiSparseReciprocal (RTLSR_C) BRDF model [54,55] to inverse LAI. They used the PROSAIL model to simulate reflectance for 397 viewing directions and established a lookup table based on a dataset of 20,000 average distributed vegetation parameter combinations, and then utilized MODIS BRDF products to simulate directional reflectance, and subsequently inverted the data to obtain estimates of LAI. There is a disadvantage to this kind of method; that is, they did not search for the optimal number and positions of viewing angles that are needed to obtain a better inversion result. There are also studies suggesting that incorporating more observation directions introduces redundant information, thereby increasing the noise in the input data and leading to greater uncertainty in the retrieval of parameters [56]. Due to the varying ability of physical models to accurately reproduce anisotropy across different viewing directions, increasing the number of observations can also lead to a decrease in inversion accuracy [40]. Hence, it is crucial to enhance LAI estimation methods by incorporating reconstructed BRDF information and determining the optimal number and positions of viewing angles. This will improve the accuracy and efficiency of large-scale LAI estimation.
This study focuses on the collaborative use of the RTLSR_C and PROSAIL models, specifically examining the PROSAIL model’s sensitivity to varying LAI at different observation positions. Our emphasis lies in ensuring consistency between the PROSAIL and RTLSR_C models across diverse observation points, with a primary goal of identifying optimal viewing angles for improved inversion results. The paper starts with an introduction highlighting the importance of LAI estimation. The methodology section details the coupling of models and the PROSAIL sensitivity analysis. Data utilization involves high-quality MODIS BRDF products and an improved look-up table for LAI inversion. Results from optimal viewing angles are presented, followed by validation against ground-based LAI measurements, high-resolution LAI maps, and MODIS LAI products.

2. Materials and Methods

The flowchart of this research is illustrated in Figure 1, which includes four main components. The first component focuses on analyzing the PROSAIL model’s sensitivity to variations. The second component investigates the coherence between PROSAIL and RTLSR_C BRDF models under different observation geometries. The third component involves using the results obtained from the previous two steps to determine the optimal direction based on the analysis conducted with the 397 viewing directions. A new lookup table that captures the relationship between vegetation parameters and reflectance in the optimal observation direction using the PROSAIL model is created. The fourth component involves simulating the reflectance in the optimal observation direction based on MODIS BRDF products. By comparing the simulated reflectance and the corresponding reflectance in the new lookup table, the minimum cost function is calculated to determine the LAI value. The results will be validated using ground-based LAI measurements and MODIS LAI products.

2.1. PROSAIL Model for Multi-Angular Reflectance Simulations

The PROSAIL model, a reliable canopy radiative transfer model, was chosen for vegetation parameter inversion due to its simple input parameters and its ability to maintain a consistent BRDF with the Ross–Li model [57]. This study references the latest version of the 4SAIL canopy BRDF model, which incorporates the hotspot effect, along with the PROSPECT-5 model that provides a comprehensive description of leaf optical reflectance and transmittance (http://teledetection.ipgp.jussieu.fr/prosail/, accessed on 10 August 2021) [58,59,60]. The PROSAIL model has the ability to simulate canopy reflectance with a resolution of 1 nm for any sun-viewing geometry within the spectra range of 400–2500 nm.
Each parameter in the PROSAIL model is associated with a common value. This study utilizes a comprehensive simulation dataset of 20,000 vegetation parameter combinations, which was used in a previous study [30]. These combinations were generated using the Satellite periodic function through the uniform sampling of seven leaf and canopy parameters, i.e., leaf structure parameter (Ns), chlorophyll a and b content (Cab), equivalent water thickness (Cw), leaf mass per unit leaf area (Cm), leaf area index (LAI), average leaf angle (ALA), and soil coefficient (Psoil). Reasonable ranges have been set for all seven parameters, and the rest of the parameters have been set to the constant values. The specific information about these leaf and canopy parameters can be found in Zhang et al., 2021 [30].
By inputting the vegetation parameters into the PROSAIL model, the canopy reflectance spectra for the entire view hemisphere at different solar zenith angles can be simulated. The red and near-infrared (NIR) bands are preferred for retrieving most LAI products due to their heightened sensitivity [61]. Consequently, PROSAIL simulations were executed at the central wavelengths within the two MODIS bands, specifically at 645 nm and 858 nm.

2.2. Kernel-Driven Ross–Li BRDF Model and MODIS BRDF

The kernel-driven Ross–Li model [42,48] has been widely adopted to reconstruct BRDF data in the sun-viewing hemisphere from limited observation. The RossThick–LiSparseReciprocal (RTLSR) model was selected as the operational algorithm for generating the MODIS BRDF parameter product [43] The general expression of the RTLSR model is provided in Equation (1). R is the surface bidirectional reflectance in the wavelength λ; fiso, fvol, and fgeo are the spectrally dependent BRDF model parameters; Kvol and Kgeo are kernel functions of the volumetric [62,63] and geometric optical [48,64] scattering, respectively; Ω is a direction vector; and i and r represent the meanings of incident and reflected radiance.
R Ω i , Ω r , λ = f i s o λ + f v o l λ K v o l Ω i , Ω r , λ + f g e o λ K g e o Ω i , Ω r , λ
Through the utilization of multi-angular reflectance data in the model, the three BRDF model parameters (fiso, fvol, and fgeo) can be derived by employing least-squares regression. Finally, by utilizing the extrapolation capability of the model, it is possible to simulate directional reflectance at arbitrary orientations [43,48].
The global 500 m BRDF parameter products provided by MODIS consist of spectrally dependent BRDF model parameters as well as quality data. In subsequent studies, only high-quality (full inversion, quality flag = 0/1) MODIS BRDF products [41,43] were used. In this study, a modified version of the Ross–Li model called RTLSR_C was employed [54,55]. The RTLSR_C model incorporates the corrected exponential hotspot function developed by Chen and Cihlar [65] into the volumetric and geometric-optical scattering kernels, resulting in an enhanced capability for simulating hotspots compared to the original RTLSR model. The optimal hotspot parameters (C1 for height and C2 for width) of the RTLSR_C model in various typical bands were determined through an extensive search using abundant hotspot data [54]. Due to the localized modification of reflectance near the hotspot, without affecting other directions, the MODIS BRDF parameters were directly applied to the hotspot-corrected kernels in order to simulate multi-angular reflectance. The hotspot parameters for MODIS at 645 nm (C1 = 0.5, C2 = 3.4°) and 858 nm (C1 = 0.5, C2 = 3.0°) were employed to simulate the directional reflectance using the MODIS BRDF parameters.

2.3. Determination of the Optimal Direction

When establishing a lookup table, it is important to choose an adequate number of directions to capture the anisotropic reflection characteristics of the Earth’s surface. However, selecting too many directions can impact computational efficiency. It is necessary to evaluate the representativeness of reflectance in each selected direction. Based on the 397 sun-viewing geometry used in the previous study [30], this study focuses on identifying the directions in which the PROSAIL model exhibits a higher sensitivity to LAI changes and demonstrates better consistency with the RTLSR_C model. The range of solar zenith angles for the 397 observations is 0°–60° with a spacing of 15°. The range of viewing zenith angles is 0°–80° with a spacing of 10°. The range of relative azimuth angles is 0°–330° with a spacing of 30°. The reason for filtering the reflectance based on these directions is that the lookup table based on these directions yields a high accuracy of LAI, and this approach could avoid reflectance correlation in neighboring directions.

2.3.1. Sensitivity Analysis of the PROSAIL Model to Changes in LAI

To determine the most sensitive directions of the PROSAIL model to changes in LAI, the dataset of 20,000 sets of vegetation parameter combinations was expanded. For each vegetation parameter combination, the LAI value was varied from 0.5 to 10 with a step size of 0.5 while keeping the other parameters unchanged. Subsequently, the expanded vegetation parameters for each sample were sequentially inputted into the PROSAIL model, and the standard deviation (σ) (Equation (2)) of the simulated reflectance ρi, Ωr) for each combination was calculated in any direction within the observed hemisphere. n represents the sample size, and in this study its value is 20. ρ ¯ is the mean value of the 20 simulated sets of reflectance ρ. A larger value of σ indicates a higher sensitivity. To ensure a more representative outcome, the study finally employed the average value of σ obtained from 20,000 datasets.
σ ( Ω i , Ω r ) = j = 1 n ρ j ( Ω i , Ω r ) ρ ( Ω i , Ω r ) ¯ 2 n 1

2.3.2. The Consistency between the Models

First, the 20,000 sets of vegetation parameters are sequentially input into the PROSAIL model to simulate the ρi, Ωr) in various directions. Then, the ρi, Ωr) is inputted into the RTLSR_C BRDF model to calculate the model parameters, and the reflectance ρ0i, Ωr) in any direction simulated by the kernel-driven model can be obtained. The root mean square error (RMSE) between the reflectance obtained by the two models (Equation (3)) is used to evaluate the consistency of the two models. The RMSE is computed for each direction using ρ and ρ₀, with a value of k set at 20,000. A lower RMSE value signifies greater consistency between the two models.
R M S E ( Ω i , Ω r ) = j = 1 k ρ j ( Ω i , Ω r ) ρ 0 , j ( Ω i , Ω r ) 2 k 1

2.4. LAI Estimation from MODIS BRDF Data

Once the optimal directions were determined, in Section 2.3, the PROSAIL model was employed to simulate the reflectance in the red and NIR bands at these directions using 20,000 different combinations of vegetation parameters. As a result, a new lookup table was created to establish the relationship between vegetation parameters and multi-angle reflectance.
When calculating LAI using MODIS BRDF, the process begins by simulating the reflectance in the red and NIR bands at the specified directions using the MODIS BRDF product and the RTLSR_C model. Subsequently, the simulated reflectance values are compared with the reflectance values produced by the PROSAIL model in the new lookup table, and the RMSEc is computed over all available viewing angles and wavelengths. The calculation formula for RMSEc is the same as Equation (3), with the difference that ρ represents the simulated reflectance in the red or NIR band at the selected direction when a specific set of vegetation parameter combinations is input into the PROSAIL model, while ρ0 represents the corresponding simulated data from MODIS BRDF. Finally, the 50 records with the lowest RMSEc values are selected, and the corresponding LAI values are averaged to derive the retrieved LAI result [30].
To improve the computational efficiency, this study also introduced the linear relationships between the Area Leaf Angle Distribution (ALA) and the fraction of vegetation volume (fvol) in the NIR band [30]. The formula is as follows:
A L A = 186.54 × f N I R v o l + 13.88   ( A L A     [ 10 ° ,   85 ° ] ,   f N I R v o l     [ 0 ,   0.3813 ] )
During the actual computation process, the size of the lookup table is constrained by incorporating a variation of ±3° around the empirically estimated Angular Leaf Alignment (ALA). To ensure that all data can participate in the calculation, when the fvol in the NIR band is not within the range of 0–0.3813, all 20,000 sets of data will be involved in the calculation.

2.5. Validation with LAI Measurements and MODIS LAI Product

We utilized the field-measured LAI values at the 500 m plot level and 30 m LAI maps to validate the proposed method. The validation process was conducted at two sites, Honghe (47°39′N, 133°31′E) and Hailun (47°24′~47°26′N, 126°47′~126°51′E), located in Heilongjiang Province, northeastern China. The crop types include rice, corn, soybean, and sorghum. A total of 180 sets of field-measured LAI data were collected by Fang et al. in 2012, 2013, and 2016 (https://doi.pangaea.de/10.1594/PANGAEA.900090, accessed on 10 August 2021) [61,66,67]. The LAI data were collected using the LAI-2200 instrument and underwent comprehensive validation during the majority of the growing seasons. Furthermore, they were utilized to validate the precision of various crucial satellite LAI products [61]. In addition to validating the accuracy of the results in this study, the field measurements were also utilized in determining the optimal number of observation directions.
LAI maps with a spatial resolution of 30 m, obtained from a previous study [61], were used alongside the corresponding LAI measurements. These high-resolution reference LAI maps were generated using HJ-1, Landsat 7, and Sentinel-2A images and demonstrated strong agreement with the field-measured LAI values. For the purpose of cross-validation, our proposed method was valid using 284 sets of aggregated LAI values at a 1.5 km scale, derived from the 30 m LAI maps in this study. The MODIS BRDF data with spatiotemporal consistency with these measurements were collected by Zhang et al. in 2021 [30]. To ensure data integrity, high-quality BRDF data from the nearest available time were incorporated as alternative data when the BRDF had a fill value or exhibited poor quality.
The algorithm’s accuracy was also validated using high-quality MODIS LAI products in two tiles. The MCD15A3H Version 6.1 MODIS Level 4 LAI product is a 4-day composite dataset with 500 m pixel size. Tile h26v04 is located in northeastern China, while tile h12v04 is situated in northeastern America. According to the statistics derived from the Annual International Geosphere-Biosphere Programme (IGBP) classification at the MODIS pixel scale, tile h26v04 is primarily characterized by grassland, comprising 53% of the tile’s coverage, along with 20% cropland. On the other hand, tile h12v04 is predominantly covered by forests, with 24% Mixed Forests, 18% Deciduous Broadleaf Forests, 8% Woody Savannas, and 5% Cropland. The selected time for h26v04 is days 181–193 of the year 2020, while the selected time for h12v04 is days 245–257. The vegetation during these selected time periods corresponds to the growing season, where the vegetation exhibits thriving and robust growth. To ensure high-quality standards, it is imperative to mandate that all MODIS LAI product data within a 16-day period exhibit high quality, with the condition that the standard deviation of the four LAI values remains below 0.4. Furthermore, the reflectance values in the red and NIR bands for 397 directions were simulated using the high quality BRDF product (MCD43A1) and the kernel-driven BRDF model at the corresponding 4 days. It is required that the RMSE between the multi-angular reflectance values of different days remains below 0.01.

3. Results

3.1. Sensitivity of the PROSAIL Model to LAI and Consistency with the Kernel-Driven BRDF Model

3.1.1. Sensitive Directions to LAI Variations

The dataset of 20,000 sets of vegetation parameter combinations was expanded by setting the LAI to range from 0.5 to 10 with a step size of 0.5. Figure 2 illustrates a three-dimensional shape of simulated reflectance under different LAI parameter conditions based on the common value and PROSAIL models with a solar zenith angle of 45°. The results demonstrate significant variations in the simulated bidirectional reflectance with changing LAI values. Figure 3 presents examples of directional reflectance variations with changing LAI values, simulated using the common value and PROSAIL model in three specific directions. It shows that as LAI is low, there are significant changes in surface reflectance with the increase in LAI. However, as LAI further increases (>3.5), the reflectance tends to remain relatively constant. It also demonstrates that the intensity of LAI variations differs across different directions and spectral bands, and the σ of reflectance in a specific direction can be used to measure the magnitude of reflectance variation. In other words, directions with higher σ values indicate that the PROSAIL model is more sensitive to LAI variations.
Figure 3 only shows an example of the distribution of σ for the common value over the viewing hemisphere. We employed a total of 20,000 sets of vegetation parameter combinations, and by calculating the average of these 20,000 sets of σ, we can determine the directions that are more sensitive to variations in LAI. The distribution of mean σ at different solar zenith angles over the viewing hemisphere is shown in Figure 4. The results indicate that in the red band, when the solar zenith angle is small, high sensitivity is observed near the small view zenith angles area. However, in the NIR band, under large solar zenith angles area, high sensitivity is observed near hotspots. The sensitivity of the PROSAIL model to changes in LAI varies across different directions in different spectral bands, highlighting the importance of considering each band individually.

3.1.2. Consistent between Two Models

The 20,000 sets of vegetation parameter combinations were sequentially inputted into the PROSAIL model to simulate the observation of hemispherical reflectance (ρ) at different solar zenith angles. These reflectance values were then used as inputs in the RTLSR_C BRDF model to calculate model parameters, allowing for the simulation of reflectance (ρ0) in any given direction. Figure 5 shows an example of the three-dimensional BRDF shapes of ρ0 based on ρ simulated using the common value and PROSAIL model. Compared to Figure 2b,e, the results indicate a high level of consistency between the BRDF simulated by the RTLSR_C BRDF model and the results obtained from the PROSAIL model for the common value dataset. Comparative analysis of these two BRDF models enables a quantitative evaluation of the consistency between them.
To assess the consistency between the two models in any direction, 20,000 sets of ρ and corresponding ρ0 from both models are compared. Figure 6 shows the comparison between the reflectance of the two models in two directions in the NIR band when the solar zenith angle is 30°. The consistency of the two models varies significantly across the two directions, which can be quantified using RMSE. Typically, the consistency in the nadir direction, which is commonly studied, is not significant, with an RMSE of 0.042. In comparison, another direction with a view zenith angle of 40° in the principal plane exhibits a lower RMSE of 0.012, indicating a stronger consistency between the two models in that particular direction.
By calculating the RMSE between the reflectance from the two models across various viewing directions, we can quantitatively evaluate the level of agreement between them. Figure 7 shows the distribution of the RMSE between the reflectance of the two models over the whole viewing hemisphere under different solar zenith angles. A lower RMSE indicates a higher level of consistency, suggesting that the two models yield similar results in terms of reflectance predictions. The consistent direction varies with solar zenith angle changes. Although efforts have been made to correct for the hot-spot effects, the consistency between the two models remains relatively low around the hot-spot direction. Additionally, in the red band, when the solar zenith angle is between 15° and 45°, the consistency is poor in some directions within the forward hemisphere. On the other hand, in the NIR band, when the solar zenith angle is greater than 30°, there is a weaker correlation near the nadir direction.

3.2. Optimal Observation Geometry for LAI Retrieval

To jointly utilize these two models for LAI inversion, it is necessary to select directions that are sensitive to LAI variations and exhibit high consistency between the two models. Based on the previous analysis using 397 directions, and considering the sensitivity of the PROSAIL model to LAI changes and model consistency, the optimized observation directions can be selected by adjusting the threshold. To determine the optimal number of observations for LAI retrieval, different numbers of observations were selected and studied, and ultimately the optimal number of observations was determined. The σ and RMSE thresholds set for different bands and observation numbers are shown in Table 1.
Figure 8 illustrates the distribution of observations in the red and NIR bands for 30, 60, and 90 selected observations. Considering both models have a symmetrical distribution about the principal plane, only half of the observation hemisphere’s viewing directions were taken into account. Similar to the previous findings, the data in the red band are concentrated in observation directions with smaller solar zenith angles. However, in the NIR band, there is a larger proportion of observations with larger solar zenith angles.
Once the observation directions were determined, the reflectance for various bands in those directions was simulated using the PROSAIL model. By incorporating the input vegetation parameters, a new lookup table was generated to establish the connection between the reflectance at different observation directions and the vegetation parameters. This lookup table will be utilized in the subsequent LAI inversion, which is based on a MODIS BRDF product.

3.3. Validation of the LAI Estimations Using Field Measurements and LAI Maps

A total of 180 field LAI measurements at the 500 m plot and the high-resolution (30 m) LAI maps were utilized for algorithm validation. High-quality MODIS BRDF products corresponding to the spatiotemporal coverage of these surface observations were employed for LAI inversion. In cases where the MODIS BRDF product quality was low, it was replaced with the nearest neighboring high-quality BRDF product to ensure data integrity.
During the LAI inversion process, the reflectance for selected observation directions was first simulated based on the MODIS BRDF product. Then, the simulated MODIS reflectance was compared with the reflectance in the new lookup table. To limit the size of the lookup table and improve the efficiency of the inversion process, the empirical ALA was calculated using the fvol-ALA relationship when the fvol in the NIR band was less than 0.3813 [30]. Finally, the RMSEc was calculated between the simulated MODIS reflectance and the reflectance values in the lookup table. The average LAI value corresponding to the 50 sets of data with the lowest RMSEc was considered to be the inversion LAI.
Figure 9 shows the comparison of 180 LAI measurements in the 500 m plot with LAI estimated from MODIS BRDF. Among them, Figure 9a is the result of previous research provided by Zhang et al., 2021 [30], using 397 observations. After the sensitivity analysis and consistency analysis, the inversion results based on 30 directions had a similar accuracy (RMSE = 1.37, Bias = −0.16, R2 = 0.43) to the result of Figure 9a. In comparison, the coefficient of determination (R2) slightly decreased, the RMSE slightly increased, but the absolute magnitude of bias decreased. Despite the significant reduction in the number of observations, the inverted LAI still exhibited a high level of consistency with the surface observation data, and also exhibited higher accuracy compared to the MODIS LAI product generated with the main algorithm, as evidenced by the evaluation metrics: RMSE of 1.50, bias of −0.23, and an R2 value of 0.24 [68]. Figure 9c shows the R2 between the results obtained based on different numbers of observations and surface observation data. The results indicate that when reducing the number of observations, the highest consistency with ground observations is achieved through inversion using 30 directions. In the following research, the retrieval of LAI will focus on utilizing the new lookup table specifically designed for the case with 30 observations.
The 30 m high-quality LAI maps were initially aggregated to a 500 m scale, and the inverted LAI values were also obtained at the 500 m resolution. The comparison of 30 m high-quality LAI maps was conducted at a scale of 1.5 km (3 × 3 MODIS pixels), resulting in a total of 284 data sets. Figure 10a is the result of previous research provided by Zhang et al., 2021 [30], using 397 observations per band; Figure 10b is the result based on 30 observations. Similar to the previous findings, despite the reduction in the number of observations used for inversion, the accuracy of the inversion results was not affected and there was a better agreement between the retrieved LAIs and the LAI maps.

3.4. Validation of LAI Estimation Based on MODIS LAI Product

To filter the selected MODIS BRDF data, it is necessary to ensure that both datasets are of high quality. Additionally, the LAI product within the 16-day period should exhibit stability, and the variations in the 397 directional reflectance data simulated by the MODIS BRDF model should be minimal. The choice of 397 directions instead of 30 directions is primarily aimed at facilitating a comparison between the proposed method and the previous approach.
The comparison between LAIs retrieval using MODIS BRDF and the MODIS LAI product is shown in Figure 11. For tile h26v04, the LAI values are relatively small, with the majority of pixels having an LAI less than 3.5. In the case of tile h12v04, the MODIS LAI values are divided into two parts: one part is below 3.5, while the other part is concentrated around 5. Compared to the inversion results based on 397 angles, the results based on 30 observation angles exhibit higher consistency. The RMSE of the two scenes decreases from 0.354 to 0.340 and from 2.424 to 2.038, respectively. The results also demonstrate that the LAI accuracy obtained from simulated multi-angle data inversion is relatively high when the LAI is less than 3.5. However, when the LAI exceeds 4.5, the LAI estimated from multi-angle data inversion tends to be significantly underestimated. This can be attributed to the use of the look-up table method, the uniform distribution of vegetation parameters employed in the PROSAIL model, and the reflectance simulated by the PROSAIL model, which shows relatively small variations when LAI is greater than 3.5. The mean value of LAI corresponding to the 50 sets of minimum RMSEc data can cause the inversion results to underestimate LAI in higher LAI areas.

4. Discussion

The accurate inversion of LAI not only plays a crucial role in advancing ecological research but also holds practical applications in improving agricultural productivity and maintaining the health of vegetation. In this paper, we examined the improvement in LAI estimation by utilizing multi-angle reflectance data and the PROSAIL model, with the application of MODIS BRDF products and an enhanced lookup table method. The results indicate that by selecting highly sensitive and consistent directions, even though the number of observation directions is reduced, the improved lookup table based on surface observations can still achieve high accuracy in LAI retrieval. Additionally, reducing the number of observation directions significantly improves the efficiency of inversion. This study also provides new insights for the coupled inversion of surface parameters using both physical and semi-empirical kernel-driven models.
The anisotropic reflectance characteristics of the Earth’s surface play a crucial role in the inversion of surface parameters, as they provide valuable information about its structural features [46]. The physical BRDF model accurately describes the surface reflection behavior. By incorporating detailed vegetation structure descriptions, this model enables the simulation of directional reflectance for various parameter combinations [60]. Through an inversion strategy, the vegetation structure can be estimated using reflectance data. Numerous studies have successfully inverted vegetation with high accuracy by leveraging observed multi-angle reflectance [31,53]. Coupling the kernel-driven and PROSAIL models by analyzing the consistency between the models as well as their sensitivity to land surface structural parameters enables the creation of a lookup table that establishes the correlation between specific direction reflectance and land surface parameters. In order to fully utilize the surface’s anisotropic reflection characteristics and provide a simplified LAI retrieval method, this study couples the PROSAIL BRDF model with the RTLSR_C model to invert LAI. Although there are discrepancies between the models, a previous study demonstrated a strong agreement in BRDF when integrating the two models, even when the PROSAIL input parameters underwent substantial variations [57]. This serves as a fundamental basis for the current study.
The most significant improvement compared to previous studies is that this study searched for the optimal number and positions of viewing angles that are needed to obtain the best inversion results. By examining the changes in directional reflectance simulated by the PROSAIL model as a function of LAI, the directions that are sensitive to changes were identified. By inputting the output of the PROSAIL model into the kernel-driven model, the consistency between the two models was explored. This process helped to determine the optimal directions and provided a basis for selecting the observation directions. In general, the optimal directions tend to exhibit a clustered distribution, meaning that under different solar zenith angles, the more sensitive and highly correlated directions are often concentrated within a smaller range. Considering the strong correlation of the model-simulated reflectance in adjacent directions, the study evaluated the existing discrete distribution of 397 observation directions [30], fully considering the surface’s anisotropic information. It is worth noting that although the improved RTLSR_C model captured the bidirectional reflectance characteristics of the hotspot, the 30 selected directions did not include the reflectance near the hotspot. This is because there were still significant differences between the two models in the surroundings of the hotspot area [54], and further research is needed to improve the characterization of the hotspot.
This study is based on high-quality MODIS 500 m BRDF products for LAI inversion. The improved algorithm significantly enhances computational efficiency. In comparison to LAI derived from surface measurements, it achieves similar levels of accuracy with an algorithm based on 397 observations [30]. Compared to the MODIS LAI product, the retrieved LAI results demonstrate better consistency with both the LAI measurement and the LAI maps. This also indicates that multi-angle observation data can provide advantageous assistance in improving the accuracy of LAI retrieval [40,69]. For high-resolution remote sensing data that have multi-angle observations, these simulated BRDF data from the RTLSR_C model hold promise as compensatory means. Certainly, further validation of this method is needed with additional surface measurements.
Based on the validation of MODIS LAI products, the improved lookup table method shows relatively better performance when LAI is less than 3.5. There are certain improvements in terms of correlation coefficient and RMSE. However, for larger LAI values, although there is some improvement, both methods yield relatively poor results. There are several reasons that contribute to this phenomenon. Firstly, the radiative transfer models like PROSAIL consist of a set of input variables that always appear in combination, making the inversion of PROSAIL an ill-posed problem [70]. This could be a common issue with lookup table methods. The results are based on the average of the 50 data sets with the lowest RMSE, which can lead to a significant underestimation of LAI values, particularly when the LAI is large. Another possible reason is that this issue is also influenced by the configuration of model parameters. In the study, the parameter combinations used were uniformly distributed [30], which did not accurately represent the actual surface characteristics. Since each parameter has an impact on the final BRDF and can introduce deviations in the results, it is crucial to optimize the parameter combinations based on the specific vegetation features of the study area [31]. This optimization process would further enhance the accuracy of the inversion results. In addition, it is important to consider that the coarse resolution of MODIS limits its ability to fully capture the detailed surface parameter characteristics of a pixel with an average parameter combination [71]. The simulated dataset based on MODIS BRDF does not fully capture the complexity and diversity of the actual land surface. Further research can explore higher-resolution BRDF data and more complex canopy structures to better understand the potential advantages of multi-angle reflectance data in LAI estimation.
In summary, the use of the improved lookup table approach, compared to the original algorithm, enables the maintenance of high LAI accuracy while improving the efficiency of the inversion process. The estimation process described in this study can be widely expanded to incorporate a range of semi-empirical kernel driven and physical BRDF models for extracting various surface parameters in large, homogeneous regions. The next step will involve optimizing the distribution of vegetation parameters based on surface characteristics and further validating the method by incorporating high-resolution BRDF data obtained from surface measurements. This will provide a more comprehensive validation of the approach and enhance the accuracy of LAI estimation.

5. Conclusions

Coupling the kernel-driven BRDF model with the PROSAIL model, and utilizing the reflectance in directions that are sensitive to LAI changes, while maintaining high consistency between two models, can improve the accuracy of LAI retrieval. In this study, we focused on determining the optimal number and positions of viewing angles required to enhance the inversion results for LAI extraction using the high-quality MODIS BRDF product and PROSAIL model. By utilizing 20,000 sets of surface vegetation parameter combinations, we explored the sensitivity of the PROSAIL model to variations in LAI and its consistency with an RTLSR_C model. Through the application of a threshold, we identified the optimal viewing directions for LAI retrieval. According to the validation using field-measured LAI data, it was found that the best results for LAI retrieval were obtained when 30 viewing angles were selected in the inversion using the MODIS BRDF model. The proposed method was further validated using field measurements and high-quality, high-resolution LAI maps, which indicated that the inversion results were comparable in accuracy to previous algorithms based on 397 viewing angles. Validation using two tiles of high-quality MODIS LAI products revealed that the proposed method exhibited a high level of consistency with the MODIS LAI product when the LAI values were less than 3.5. However, in forest-dominated regions where LAI values were higher, although there was some improvement compared to previous methods, the extent of improvement was relatively weak. This may be attributed to the PROSAIL model’s ability to capture the anisotropic characteristics of the Earth’s surface, the algorithm’s utilization of uniformly distributed input parameters, and the minimal variations in PROSAIL-simulated reflectance at high LAI.
In the future, unmanned aerial vehicles will be employed to acquire high-resolution surface BRDF datasets. The focus will be on redesigning the distribution of vegetation parameter combinations based on the actual characteristics of surface vegetation. Furthermore, the algorithm will be enhanced by integrating the kernel-driven BRDF model with additional physical models, aiming to improve the inversion accuracy in regions with high LAI values.

Author Contributions

Conceptualization, H.Z. and Y.L. (Yan Liu); methodology, H.Z. and X.Z.; validation, Y.L. (Yi Lian) and H.C.; formal analysis, L.C. (Lei Cui); investigation, Y.D.; data curation, L.C. (Lei Chen), X.Z. and Q.X.; writing—original draft preparation, H.Z.; writing—review and editing, Y.L. (Yan Liu), X.Z. and Y.D.; funding acquisition, H.Z. and Y.L. (Yan Liu). All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (41971306).

Data Availability Statement

All satellite remote sensing and field measured data used in this study are openly and freely available. LAI measurements at the 500 m plot level and 30 m LAI maps can be accessed via https://doi.pangaea.de/10.1594/PANGAEA.900090 (accessed on 10 August 2021). The Collection 6 MODIS BRDF parameter product (MCD43A1 and MCD43A2) is available at https://search.earthdata.nasa.gov/search (accessed on 10 March 2023).

Acknowledgments

We express our sincere gratitude to the research team of Hongliang Fang at the Chinese Academy of Sciences for generously sharing their LAI measurements as public data. Our gratitude also goes to the PROSAIL team for making their model code freely accessible online. Additionally, we appreciate the research team of the MODIS BRDF products for freely providing global data. Lastly, we are thankful for the thorough review and valuable feedback provided by the anonymous reviewers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yan, G.; Hu, R.; Luo, J.; Weiss, M.; Jiang, H.; Mu, X.; Xie, D.; Zhang, W. Review of indirect optical measurements of leaf area index: Recent advances, challenges, and perspectives. Agric. For. Meteorol. 2019, 265, 390–411. [Google Scholar] [CrossRef]
  2. Watson, D.J. Comparative Physiological Studies on the Growth of Field Crops: I. Variation in Net Assimilation Rate and Leaf Area between Species and Varieties, and within and between Years. Ann. Bot. 1947, 11, 41–76. [Google Scholar] [CrossRef]
  3. Chen, J.M.; Black, T.A. Measuring leaf area index of plant canopies with branch architecture. Agric. For. Meteorol. 1991, 57, 1–12. [Google Scholar] [CrossRef]
  4. Wilson, T.B.; Meyers, T.P. Determining vegetation indices from solar and photosynthetically active radiation fluxes. Agric. For. Meteorol. 2007, 144, 160–179. [Google Scholar] [CrossRef]
  5. Duchemin, B.; Hadria, R.; Erraki, S.; Boulet, G.; Maisongrande, P.; Chehbouni, A.; Escadafal, R.; Ezzahar, J.; Hoedjes, J.; Kharrou, M. Monitoring wheat phenology and irrigation in Central Morocco: On the use of relationships between evapotranspiration, crops coefficients, leaf area index and remotely-sensed vegetation indices. Agric. Water Manag. 2006, 79, 1–27. [Google Scholar] [CrossRef]
  6. Cleugh, H.A.; Leuning, R.; Mu, Q.; Running, S.W. Regional evaporation estimates from flux tower and MODIS satellite data. Remote Sens. Environ. 2007, 106, 285–304. [Google Scholar] [CrossRef]
  7. Chen, J.M.; Chen, X.; Ju, W.; Geng, X. Distributed hydrological model for mapping evapotranspiration using remote sensing inputs. J. Hydrol. 2005, 305, 15–39. [Google Scholar] [CrossRef]
  8. Chen, B.; Chen, J.M.; Ju, W. Remote sensing-based ecosystem–atmosphere simulation scheme (EASS)—Model formulation and test with multiple-year data. Ecol. Model. 2007, 209, 277–300. [Google Scholar] [CrossRef]
  9. Zhao, J.; Li, J.; Liu, Q.; Xu, B.; Hu, Z. Estimating fractional vegetation cover from leaf area index and clumping index based on the gap probability theory. Int. J. Appl. Earth Obs. Geoinf. 2020, 90, 102112. [Google Scholar] [CrossRef]
  10. Clark, D.B.; Olivas, P.C.; Oberbauer, S.F.; Clark, D.A.; Ryan, M.G. First direct landscape-scale measurement of tropical rain forest Leaf Area Index, a key driver of global primary productivity. Ecol. Lett. 2008, 11, 163–172. [Google Scholar] [CrossRef]
  11. Jin, H.; Li, A.; Wang, J.; Bo, Y. Improvement of spatially and temporally continuous crop leaf area index by integration of CERES-Maize model and MODIS data. Eur. J. Agron. 2016, 78, 1–12. [Google Scholar] [CrossRef]
  12. Belen, F.; Eric, V.; Jean-Claude, R.; Emilie, M.; Inbal, B.R.; Chris, J.; Martin, C.; Jyoteshwar, N.; Ivan, C.; Dave, M. A 30+Year AVHRR Land Surface Reflectance Climate Data Record and Its Application to Wheat Yield Monitoring. Remote Sens. 2017, 9, 296. [Google Scholar]
  13. 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]
  14. Asner, G.; Scurlock, J.; Hicke, J. Global synthesis of leaf area index observations: Implications for ecological and remote sensing studies: Global leaf area index. Glob. Ecol. Biogeogr. 2003, 12, 191–205. [Google Scholar] [CrossRef]
  15. Tian, Y.; Zheng, Y.; Zheng, C.; Xiao, H.; Fan, W.; Zou, S.; Wu, B.; Yao, Y.; Zhang, A.; Liu, J. Exploring scale-dependent ecohydrological responses in a large endorheic river basin through integrated surface water-groundwater modeling. Water Resour. Res. 2015, 51, 4065–4085. [Google Scholar] [CrossRef]
  16. 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]
  17. Mora, M.; Avila, F.; Carrasco-Benavides, M.; Maldonado, G.; Olguín-Cáceres, J.; Fuentes, S. Automated computation of leaf area index from fruit trees using improved image processing algorithms applied to canopy cover digital photograpies. Comput. Electron. Agric. 2016, 123, 195–202. [Google Scholar] [CrossRef]
  18. Rahimikhoob, H.; Delshad, M.; Habibi, R. Leaf area estimation in lettuce: Comparison of artificial intelligence-based methods with image analysis technique. Measurement 2023, 222, 113636. [Google Scholar] [CrossRef]
  19. Jonckheere, I.; Fleck, S.; Nackaerts, K.; Muys, B.; Coppin, P.; Weiss, M.; Baret, F. Review of methods for in situ leaf area index determination: Part I. Theories, sensors and hemispherical photography. Agric. For. Meteorol. 2004, 121, 19–35. [Google Scholar] [CrossRef]
  20. 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]
  21. Xiao, Z.; Liang, S.; Wang, J.; Xiang, Y.; Zhao, X.; Song, J. Long-Time-Series Global Land Surface Satellite Leaf Area Index Product Derived from MODIS and AVHRR Surface Reflectance. IEEE Trans. Geosci. Remote Sens. 2016, 54, 5301–5318. [Google Scholar] [CrossRef]
  22. Baret, F.; Hagolle, O.; Geiger, B.; Bicheron, P.; Miras, B.; Huc, M.; Berthelot, B.; Niño, 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]
  23. Tum, M.; Günther, K.P.; Böttcher, M.; Baret, F.; Bittner, M.; Brockmann, C.; Weiss, M. Global Gap-Free MERIS LAI Time Series (2002–2012). Remote Sens. 2016, 8, 69. [Google Scholar] [CrossRef]
  24. Yan, K.; Park, T.; Chen, C.; Xu, B.; Song, W.; Yang, B.; Zeng, Y.; Liu, Z.; Yan, G.; Knyazikhin, Y.; et al. Generating Global Products of LAI and FPAR From SNPP-VIIRS Data: Theoretical Background and Implementation. IEEE Trans. Geosci. Remote Sens. 2018, 56, 2119–2137. [Google Scholar] [CrossRef]
  25. Garrigues, S.; Lacaze, R.; Baret, F.; Morisette, J.T.; Weiss, M.; Nickeson, J.E.; Fernandes, R.; Plummer, S.; Shabanov, N.V.; Myneni, R.B.; et al. Validation and intercomparison of global Leaf Area Index products derived from remote sensing data. J. Geophys. Res. Biogeosci. 2008, 113, G02028. [Google Scholar] [CrossRef]
  26. Fang, H.; Wei, S.; Liang, S. Validation of MODIS and CYCLOPES LAI products using global field measurement data. Remote Sens. Environ. 2012, 119, 43–54. [Google Scholar] [CrossRef]
  27. Claverie, M.; Vermote, E.F.; Weiss, M.; Baret, F.; Hagolle, O.; Demarez, V. Validation of coarse spatial resolution LAI and FAPAR time series over cropland in southwest France. Remote Sens. Environ. 2013, 139, 216–230. [Google Scholar] [CrossRef]
  28. 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]
  29. Xiao, Z.; Liang, S.; Jiang, B. Evaluation of four long time-series global leaf area index products. Agric. For. Meteorol. 2017, 246, 218–230. [Google Scholar] [CrossRef]
  30. Zhang, X.; Jiao, Z.; Zhao, C.; Yin, S.; Cui, L.; Dong, Y.; Zhang, H.; Guo, J.; Xie, R.; Li, S.; et al. Retrieval of Leaf Area Index by Linking the PROSAIL and Ross-Li BRDF Models Using MODIS BRDF Data. Remote Sens. 2021, 13, 4911. [Google Scholar] [CrossRef]
  31. Roosjen, P.P.J.; Brede, B.; Suomalainen, J.M.; Bartholomeus, H.M.; Kooistra, L.; Clevers, J.G.P.W. Improved estimation of leaf area index and leaf chlorophyll content of a potato crop using multi-angle spectral data—Potential of unmanned aerial vehicle imagery. Int. J. Appl. Earth Obs. Geoinf. 2018, 66, 14–26. [Google Scholar] [CrossRef]
  32. Xiao, Z.; Liang, S.; Wang, J.; Xie, D.; Song, J.; Fensholt, R. A Framework for Consistent Estimation of Leaf Area Index, Fraction of Absorbed Photosynthetically Active Radiation, and Surface Albedo from MODIS Time-Series Data. IEEE Trans. Geosci. Remote Sens. 2015, 53, 3178–3197. [Google Scholar] [CrossRef]
  33. Jacquemoud, S.; Bacour, C.; Poilvé, H.; Frangi, J.P. Comparison of Four Radiative Transfer Models to Simulate Plant Canopies Reflectance: Direct and Inverse Mode. Remote Sens. Environ. 2000, 74, 471–481. [Google Scholar] [CrossRef]
  34. 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]
  35. Jacquemoud, S.; Baret, F. PROSPECT: A model of leaf optical properties spectra. Remote Sens. Environ. 1990, 34, 75–91. [Google Scholar] [CrossRef]
  36. Darvishzadeh, R.; Skidmore, A.; Schlerf, M.; Atzberger, C. Inversion of a radiative transfer model for estimating vegetation LAI and chlorophyll in a heterogeneous grassland. Remote Sens. Environ. 2008, 112, 2592–2604. [Google Scholar] [CrossRef]
  37. Chen, J.M.; Liu, J.; Leblanc, S.G.; Lacaze, R.; Roujean, J.-L. Multi-angular optical remote sensing for assessing vegetation structure and carbon absorption. Remote Sens. Environ. 2003, 84, 516–525. [Google Scholar] [CrossRef]
  38. Widlowski, J.L.; Pinty, B.; Gobron, N.; Verstraete, M.M.; Diner, D.J.; Davis, A.B. Canopy Structure Parameters Derived from Multi-Angular Remote Sensing Data for Terrestrial Carbon Studies. Clim. Change 2004, 67, 403–415. [Google Scholar] [CrossRef]
  39. Duan, S.; Li, Z.; Wu, H.; Tang, B.; Ma, L.; Zhao, E.; Li, C. Inversion of the PROSAIL model to estimate leaf area index of maize, potato, and sunflower fields from unmanned aerial vehicle hyperspectral data. Int. J. Appl. Earth Obs. Geoinf. 2014, 26, 12–20. [Google Scholar] [CrossRef]
  40. Dorigo, W.A. Improving the Robustness of Cotton Status Characterisation by Radiative Transfer Model Inversion of Multi-Angular CHRIS/PROBA Data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2012, 5, 18–29. [Google Scholar] [CrossRef]
  41. Jin, Y.; Schaaf, C.B.; Gao, F.; Li, X.; Strahler, A.H.; Lucht, W.; Liang, S. Consistency of MODIS surface bidirectional reflectance distribution function and albedo retrievals: 1. Algorithm performance. J. Geophys. Res. Atmos. 2003, 108, 4158. [Google Scholar] [CrossRef]
  42. Roujean, J.L.; Leroy, M.; Deschanps, P.Y. A bidirectional reflectance model of the Earth's surface for the correction of remote sensing data. J. Geophys. Res. Atmos. 1992, 972, 20455–20468. [Google Scholar] [CrossRef]
  43. Schaaf, C.B.; Gao, F.; Strahler, A.H.; Lucht, W.; Li, X.; Tsang, T.; Strugnell, N.C.; Zhang, X.; Jin, Y.; Muller, J.P. First operational BRDF, albedo nadir reflectance products from MODIS. Remote Sens. Environ. 2002, 83, 135–148. [Google Scholar] [CrossRef]
  44. Lucht, W.; Lewis, P. Theoretical noise sensitivity of BRDF and albedo retrieval from the EOS-MODIS and MISR sensors with respect to angular sampling. Int. J. Remote Sens. 2000, 21, 81–98. [Google Scholar] [CrossRef]
  45. Chen, J.M.; Li, X.; Nilson, T.; Strahler, A. Recent advances in geometrical optical modelling and its applications. Remote Sens. Rev. 2000, 18, 227–262. [Google Scholar] [CrossRef]
  46. Strahler, A.H. Vegetation canopy reflectance modeling—Recent developments and remote sensing perspectives. Remote Sens. Rev. 1997, 15, 179–194. [Google Scholar] [CrossRef]
  47. Qi, J.; Cabot, F.; Moran, M.S.; Dedieu, G. Biophysical parameter estimations using multidirectional spectral measurements. Remote Sens. Environ. 1995, 54, 71–83. [Google Scholar] [CrossRef]
  48. Wanner, W.; Li, X.; Strahler, A.H. On the derivation of kernels for kernel-driven models of bidirectional reflectance. J. Geophys. Res. Atmos. 1995, 100, 21077–21089. [Google Scholar] [CrossRef]
  49. He, T.; Liang, S.; Wang, D.; Cao, Y.; Gao, F.; Yu, Y.; Feng, M. Evaluating land surface albedo estimation from Landsat MSS, TM, ETM+, and OLI data based on the unified direct estimation approach. Remote Sens. Environ. 2018, 204, 181–196. [Google Scholar] [CrossRef]
  50. Qu, Y.; Liu, Q.; Liang, S.; Wang, L.; Liu, N.; Liu, S. Direct-Estimation Algorithm for Mapping Daily Land-Surface Broadband Albedo from MODIS Data. IEEE Trans. Geosci. Remote Sens. 2014, 52, 907–919. [Google Scholar] [CrossRef]
  51. Jiao, Z.; Zhang, H.; Dong, Y.; Liu, Q.; Xiao, Q.; Li, X. An Algorithm for Retrieval of Surface Albedo from Small View-Angle Airborne Observations through the Use of BRDF Archetypes as Prior Knowledge. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2015, 8, 3279–3293. [Google Scholar] [CrossRef]
  52. Jiao, Z.; Dong, Y.; Schaaf, C.B.; Chen, J.M.; Román, M.; Wang, Z.; Zhang, H.; Ding, A.; Erb, A.; Hill, M.J. An algorithm for the retrieval of the clumping index (CI) from the MODIS BRDF product using an adjusted version of the kernel-driven BRDF model. Remote Sens. Environ. 2018, 209, 594–611. [Google Scholar] [CrossRef]
  53. Wang, Z.; Schaaf, C.B.; Lewis, P.; Knyazikhin, Y.; Schull, M.A.; Strahler, A.H.; Yao, T.; Myneni, R.B.; Chopping, M.J.; Blair, B.J. Retrieval of canopy height using moderate-resolution imaging spectroradiometer (MODIS) data. Remote Sens. Environ. 2011, 115, 1595–1601. [Google Scholar] [CrossRef]
  54. Dong, Y.; Jiao, Z.; Cui, L.; Zhang, H.; Guo, J. Assessment of the Hotspot Effect for the PROSAIL Model with POLDER Hotspot Observations Based on the Hotspot-Enhanced Kernel-Driven BRDF Model. IEEE Trans. Geosci. Remote Sens. 2019, 57, 8048–8064. [Google Scholar] [CrossRef]
  55. Jiao, Z.; Schaaf, C.B.; Dong, Y.; Román, M.; Hill, M.J.; Chen, J.M.; Wang, Z.; Zhang, H.; Saenz, E.; Poudyal, R. A method for improving hotspot directional signatures in BRDF models used for MODIS. Remote Sens. Environ. 2016, 186, 135–151. [Google Scholar] [CrossRef]
  56. 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]
  57. Zhang, X.; Jiao, Z.; Dong, Y.; Zhang, H.; Li, Y.; He, D.; Ding, A.; Yin, S.; Cui, L.; Chang, Y. Potential Investigation of Linking PROSAIL with the Ross-Li BRDF Model for Vegetation Characterization. Remote Sens. 2018, 10, 437. [Google Scholar] [CrossRef]
  58. Wang, L.; Dong, T.; Zhang, G.; Niu, Z. LAI Retrieval Using PROSAIL Model and Optimal Angle Combination of Multi-Angular Data in Wheat. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2013, 6, 1730–1736. [Google Scholar] [CrossRef]
  59. 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]
  60. Jacquemoud, S.; Verhoef, W.; Baret, F.; Bacour, C.; Zarco-Tejada, 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]
  61. Fang, H.; Zhang, Y.; Wei, S.; Li, W.; Ye, Y.; Sun, T.; Liu, W. Validation of global moderate resolution leaf area index (LAI) products over croplands in northeastern China. Remote Sens. Environ. 2019, 233, 111377. [Google Scholar] [CrossRef]
  62. Myneni, R.B.; Asrar, G.; Hall, F.G. A three-dimensional radiative transfer method for optical remote sensing of vegetated land surfaces. Remote Sens. Environ. 1992, 41, 105–121. [Google Scholar] [CrossRef]
  63. Ross, J.K. The Radiation Regime and Architecture of Plant Stands; Dr. W. Junk: Norwell, MA, USA, 1981; p. 392. [Google Scholar]
  64. Li, X.; Strahler, A.H. Geometric-optical bidirectional reflectance modeling of the discrete crown vegetation canopy: Effect of crown shape and mutual shadowing. IEEE Trans. Geosci. Remote Sens. 1992, 30, 276–292. [Google Scholar] [CrossRef]
  65. Chen, J.M.; Cihlar, J. A hotspot function in a simple bidirectional reflectance model for satellite applications. J. Geophys. Res. Atmos. 1997, 102, 25907–25913. [Google Scholar] [CrossRef]
  66. Fang, H.; Li, W.; Wei, S.; Jiang, C. Seasonal variation of leaf area index (LAI) over paddy rice fields in NE China: Intercomparison of destructive sampling, LAI-2200, digital hemispherical photography (DHP), and AccuPAR methods. Agric. For. Meteorol. 2014, 198–199, 126–141. [Google Scholar] [CrossRef]
  67. Fang, H.; Ye, Y.; Liu, W.; Wei, 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–254, 48–61. [Google Scholar] [CrossRef]
  68. Yan, K.; Park, T.; Yan, G.; Liu, Z.; Yang, B.; Chen, C.; Nemani, R.R.; Knyazikhin, Y.; Myneni, R.B. Evaluation of MODIS LAI/FPAR product collection 6. Part 2: Validation and intercomparison. Remote Sens. 2016, 8, 460. [Google Scholar] [CrossRef]
  69. Schaepman, M.E.; Koetz, B.; Schaepman-Strub, G.; Itten, K.I. Spectrodirectional remote sensing for the improved estimation of biophysical and -chemical variables: Two case studies. Int. J. Appl. Earth Obs. Geoinf. 2005, 6, 271–282. [Google Scholar] [CrossRef]
  70. 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]
  71. Pan, F.; Wu, X.; Tang, R.; Zeng, Q.; Li, Z.; Wang, J.; You, D.; Wen, J.; Xiao, Q. A coarse pixel scale ground “truth” dataset based on the global in situ site measurements to support validation and bias correction of satellite surface albedo products. Earth Syst. Sci. Data Discuss. 2023, 2023, 1–25. [Google Scholar]
Figure 1. Flowchart of LAI estimation by linking the PROSAIL and Ross–Li BRDF models using the MODIS BRDF product. Part A is the PROSAIL model’s sensitivity to LAI, part B is the coherence between PROSAIL and RTLSR_C BRDF models, part C is to determine the optimal direction, and part D is the inversion and validation of LAI based on MODIS BRDF.
Figure 1. Flowchart of LAI estimation by linking the PROSAIL and Ross–Li BRDF models using the MODIS BRDF product. Part A is the PROSAIL model’s sensitivity to LAI, part B is the coherence between PROSAIL and RTLSR_C BRDF models, part C is to determine the optimal direction, and part D is the inversion and validation of LAI based on MODIS BRDF.
Remotesensing 15 05609 g001
Figure 2. Three-dimensional BRDF shape in the red (ac) and NIR (df) bands simulated based on the common value and PROSAIL model under different LAI parameter conditions. Different colors represent the magnitude of the reflectance. In the bottom coordinate plane, the radius represents the zenith angle, while the polar angle represents the azimuth angle. The vertical axis is used to plot the BRDF values.
Figure 2. Three-dimensional BRDF shape in the red (ac) and NIR (df) bands simulated based on the common value and PROSAIL model under different LAI parameter conditions. Different colors represent the magnitude of the reflectance. In the bottom coordinate plane, the radius represents the zenith angle, while the polar angle represents the azimuth angle. The vertical axis is used to plot the BRDF values.
Remotesensing 15 05609 g002
Figure 3. Here, (a,c) refer to the directional reflectance variations with changing LAI values, simulated using the common value and PROSAIL model in the nadir (green line), 45° forward (45F, black line) and 45° backward (45B, red line) in the red and NIR bands under a solar zenith angle of 45°. The σ of directional reflectance with changing LAI is also demonstrated. Then, (b,d) refer to the distribution of σ in the viewing hemisphere. The radius represents the zenith angle, while the polar angle represents the azimuth angle. Different colors represent the magnitude of the σ.
Figure 3. Here, (a,c) refer to the directional reflectance variations with changing LAI values, simulated using the common value and PROSAIL model in the nadir (green line), 45° forward (45F, black line) and 45° backward (45B, red line) in the red and NIR bands under a solar zenith angle of 45°. The σ of directional reflectance with changing LAI is also demonstrated. Then, (b,d) refer to the distribution of σ in the viewing hemisphere. The radius represents the zenith angle, while the polar angle represents the azimuth angle. Different colors represent the magnitude of the σ.
Remotesensing 15 05609 g003
Figure 4. The distribution of the average σ of 20,000 sets of data over the viewing hemisphere at solar zenith angles (SZA) of 0° (a,f), 15° (b,g), 30° (c,h), 45° (d,i), and 60° (e,j) in the red (ae) and NIR (fj) bands.
Figure 4. The distribution of the average σ of 20,000 sets of data over the viewing hemisphere at solar zenith angles (SZA) of 0° (a,f), 15° (b,g), 30° (c,h), 45° (d,i), and 60° (e,j) in the red (ae) and NIR (fj) bands.
Remotesensing 15 05609 g004
Figure 5. The three-dimensional BRDF shapes based on the RTLSR_C model and the multi-angular reflectance simulated using the Common Value and PROSAIL models in the red (a) and NIR (b) bands.
Figure 5. The three-dimensional BRDF shapes based on the RTLSR_C model and the multi-angular reflectance simulated using the Common Value and PROSAIL models in the red (a) and NIR (b) bands.
Remotesensing 15 05609 g005
Figure 6. Comparison of reflectance in the NIR band between the PROSAIL and RTLSR_C models in the nadir direction (black points) and at a backscattering angle of 40° (red points) when the solar zenith angle is 30°.
Figure 6. Comparison of reflectance in the NIR band between the PROSAIL and RTLSR_C models in the nadir direction (black points) and at a backscattering angle of 40° (red points) when the solar zenith angle is 30°.
Remotesensing 15 05609 g006
Figure 7. The distribution of the RMSE between the reflectance from the PROSAIL and RTLSR_C models over the whole viewing hemisphere at solar zenith angles of 0° (a,f), 15° (b,g), 30° (c,h), 15° (d,i), and 60° (e,j) in the red (ae) and NIR (fj) bands.
Figure 7. The distribution of the RMSE between the reflectance from the PROSAIL and RTLSR_C models over the whole viewing hemisphere at solar zenith angles of 0° (a,f), 15° (b,g), 30° (c,h), 15° (d,i), and 60° (e,j) in the red (ae) and NIR (fj) bands.
Remotesensing 15 05609 g007
Figure 8. The distribution of observations in the red (ac) and NIR (df) bands for 30 (a,d), 60 (b,e), and 90 (c,f) selected directions. Solid circles represent the position of the sun, while hollow rectangles represent the observation locations. The colors red, green, blue, gold, and purple correspond to observations taken at solar zenith angles of 0°, 15°, 30°, 45°, and 60°, respectively. Different-sized rectangles are used to indicate the positions of the observation data in order to avoid overlap.
Figure 8. The distribution of observations in the red (ac) and NIR (df) bands for 30 (a,d), 60 (b,e), and 90 (c,f) selected directions. Solid circles represent the position of the sun, while hollow rectangles represent the observation locations. The colors red, green, blue, gold, and purple correspond to observations taken at solar zenith angles of 0°, 15°, 30°, 45°, and 60°, respectively. Different-sized rectangles are used to indicate the positions of the observation data in order to avoid overlap.
Remotesensing 15 05609 g008
Figure 9. The comparison of 180 LAI measurements at the 500 m plot with LAI estimated from MODIS BRDF based on different numbers of observations. Among them, (a) is the result based on 397 observations; (b) is the result based on 30 directions; and (c) is the R2 between the results obtained based on different numbers of observations and surface observation data. The blue line in (a,b) represents the fitted relationship between these two sets of data.
Figure 9. The comparison of 180 LAI measurements at the 500 m plot with LAI estimated from MODIS BRDF based on different numbers of observations. Among them, (a) is the result based on 397 observations; (b) is the result based on 30 directions; and (c) is the R2 between the results obtained based on different numbers of observations and surface observation data. The blue line in (a,b) represents the fitted relationship between these two sets of data.
Remotesensing 15 05609 g009
Figure 10. The comparison of 30 m high-quality LAI maps with LAIs retrieval from MODIS BRDF at a scale of 1.5 km: (a) is the result based on 397 observations provided by Zhang et al., 2021, [10] and (b) is the result based on 30 selected observations. The blue line represents the fitted relationship between these two sets of data.
Figure 10. The comparison of 30 m high-quality LAI maps with LAIs retrieval from MODIS BRDF at a scale of 1.5 km: (a) is the result based on 397 observations provided by Zhang et al., 2021, [10] and (b) is the result based on 30 selected observations. The blue line represents the fitted relationship between these two sets of data.
Remotesensing 15 05609 g010
Figure 11. The comparison between LAI retrieval using MODIS BRDF and the MODIS LAI product. Specifically, (a,b) show the results for tile h26v04 during days 181–193 of the year 2020, while (c,d) display the results for tile h12v04 during days 245–257 of the year 2020. Here, (a,c) refer to LAI retrieval from 397 directional reflectance per band, while (b,d) refer to LAI retrieval from 30 directional reflectance per band. Different colors represent the number of pixels in each category, with gray indicating a pixel number of less than 5.
Figure 11. The comparison between LAI retrieval using MODIS BRDF and the MODIS LAI product. Specifically, (a,b) show the results for tile h26v04 during days 181–193 of the year 2020, while (c,d) display the results for tile h12v04 during days 245–257 of the year 2020. Here, (a,c) refer to LAI retrieval from 397 directional reflectance per band, while (b,d) refer to LAI retrieval from 30 directional reflectance per band. Different colors represent the number of pixels in each category, with gray indicating a pixel number of less than 5.
Remotesensing 15 05609 g011
Table 1. The σ and RMSE thresholds set for the red and NIR bands with respect to different observation numbers.
Table 1. The σ and RMSE thresholds set for the red and NIR bands with respect to different observation numbers.
No.RedNIR
σRMSEσRMSE
150.02350.00380.08600.0148
300.02200.00400.08000.0150
600.02010.00400.07500.0150
900.01810.00410.06760.0152
1500.01450.00450.05000.0325
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zhang, H.; Zhang, X.; Cui, L.; Dong, Y.; Liu, Y.; Xi, Q.; Cao, H.; Chen, L.; Lian, Y. Enhancing Leaf Area Index Estimation with MODIS BRDF Data by Optimizing Directional Observations and Integrating PROSAIL and Ross–Li Models. Remote Sens. 2023, 15, 5609. https://doi.org/10.3390/rs15235609

AMA Style

Zhang H, Zhang X, Cui L, Dong Y, Liu Y, Xi Q, Cao H, Chen L, Lian Y. Enhancing Leaf Area Index Estimation with MODIS BRDF Data by Optimizing Directional Observations and Integrating PROSAIL and Ross–Li Models. Remote Sensing. 2023; 15(23):5609. https://doi.org/10.3390/rs15235609

Chicago/Turabian Style

Zhang, Hu, Xiaoning Zhang, Lei Cui, Yadong Dong, Yan Liu, Qianrui Xi, Hongtao Cao, Lei Chen, and Yi Lian. 2023. "Enhancing Leaf Area Index Estimation with MODIS BRDF Data by Optimizing Directional Observations and Integrating PROSAIL and Ross–Li Models" Remote Sensing 15, no. 23: 5609. https://doi.org/10.3390/rs15235609

APA Style

Zhang, H., Zhang, X., Cui, L., Dong, Y., Liu, Y., Xi, Q., Cao, H., Chen, L., & Lian, Y. (2023). Enhancing Leaf Area Index Estimation with MODIS BRDF Data by Optimizing Directional Observations and Integrating PROSAIL and Ross–Li Models. Remote Sensing, 15(23), 5609. https://doi.org/10.3390/rs15235609

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