Research Article  Open Access
Guangjian Yan, Hailan Jiang, Jinghui Luo, Xihan Mu, Fan Li, Jianbo Qi, Ronghai Hu, Donghui Xie, Guoqing Zhou, "Quantitative Evaluation of Leaf Inclination Angle Distribution on Leaf Area Index Retrieval of Coniferous Canopies", Journal of Remote Sensing, vol. 2021, Article ID 2708904, 15 pages, 2021. https://doi.org/10.34133/2021/2708904
Quantitative Evaluation of Leaf Inclination Angle Distribution on Leaf Area Index Retrieval of Coniferous Canopies
Abstract
Both leaf inclination angle distribution (LAD) and leaf area index (LAI) dominate optical remote sensing signals. The Gfunction, which is a function of LAD and remote sensing geometry, is often set to 0.5 in the LAI retrieval of coniferous canopies even though this assumption is only valid for spherical LAD. Large uncertainties are thus introduced. However, because numerous tiny leaves grow on conifers, it is nearly impossible to quantitatively evaluate such uncertainties in LAI retrieval. In this study, we proposed a method to characterize the possible change of Gfunction of coniferous canopies as well as its effect on LAI retrieval. Specifically, a MultiDirectional Imager (MDI) was developed to capture stereo images of the branches, and the needles were reconstructed. The accuracy of the inclination angles calculated from the reconstructed needles was high. Moreover, we analyzed whether a spherical distribution is a valid assumption for coniferous canopies by calculating the possible range of the Gfunction from the measured LADs of branches of Larch and Spruce and the true Gfunctions of other species from some existing inventory data and threedimensional (3D) tree models. Results show that the constant G assumption introduces large errors in LAI retrieval, which could be as large as 53% in the zenithal viewing direction used by spaceborne LiDAR. As a result, accurate LAD estimation is recommended. In the absence of such data, our results show that a viewing zenith angle between 45 and 65 degrees is a good choice, at which the errors of LAI retrieval caused by the spherical assumption will be less than 10% for coniferous canopies.
1. Introduction
Leaf area index (LAI) is an important parameter of vegetation canopy structure and is critical in the quantitative calculation of global mass and energy exchange between the biosphere and the atmosphere [1, 2]. Beer’s law [3] is widely used to retrieve LAI from gap fraction. An important factor that highly affects the retrieval accuracy of LAI based on Beer’s law is the leaf projection coefficient (Gfunction). The Gfunction is the effective projection proportion of leaves along the observing directions [4] and is dependent on the leaf inclination angle distribution (LAD) (assuming a uniform azimuth angle). Previous studies showed that LAD could be highly variable depending on a number of factors such as tree species, light exposure, growing season, and canopy height [5–9]. However, due to the inconvenience of in situ measurements, the spherical distribution assumption () of LAD is often used in LAI retrieval [10–13]. Some groundbased measurements utilize the gap fraction data close to a view zenith angle of 57.5°, where the values of the Gfunction with various LADs are close to 0.5 [14–17]. However, assuming G to be 0.5 may introduce large errors to the retrieval of LAI in other view directions.
The LAD measurement methods can be classified into two categories: direct and indirect. Direct methods measure leaf inclination angles by direct contact. A protractor with a lead hammer, an inclinometer, or a digitizer is used to measure the inclination angles of a certain number of leaves [18, 19]. These methods are timeconsuming, laborintensive, and fraught with potential errors [20, 21]. Indirect methods can be divided into several types: (i) a digital photography measurement was proposed by Ryu et al. (2010) [22]. It uses a horizontal digital camera to view canopies and then manually measure the inclination angles of leaves which are approximately perpendicular to the viewing direction of the camera. This method was used in many studies and was then applied to Unmanned Aerial Vehicles (UAV) platforms [9, 23–25]. However, a large amount of humancomputer interactions are required, and it is difficult to measure the leaf inclination angles on tall forest canopies. (ii) Another method was proposed to derive leaf inclination angles based on the relationship between light transmittance and LADs by using Beer’s law under the assumption of randomly distributed leaves [26, 27]. However, this method only provides the statistical average of leaf inclination angles rather than detailed LAD. In addition, the accuracy is affected by the nonrandom distribution of leaves and light conditions. (iii) With the development of Light Detection And Ranging (LiDAR) techniques, the use of Terrestrial Laser Scanning (TLS) has become more prevalent [21, 28–32]. However, it is difficult to extract individual needles from the LiDAR point clouds. As a result, current TLSbased methods mainly focus on the broadleaved trees. (iv) Stereo imaging, which can be used to reconstruct vegetation canopies with photos taken from two or more positions based on the principle of geometric perspective, is a suitable approach to extract LAD. Wang et al. (2009) [33] reconstructed a threedimensional (3D) corn model with photos and measured its LAD. MüllerLinow et al. (2015) [34] derived 3D surface models of plants from stereo images taken by two cameras aligned in a fixed geometry and computed geometric properties from them. Mobile phone photography is more portable and efficient compared with using the professional cameras. Qi et al. (2019) [35] proposed a new method to reconstruct 3D models of leaves using multiangle mobile photos that can be used to calculate accurate LADs.
The above methods are suitable mainly for broadleaf canopies. Convenient and efficient methods for measuring LADs for coniferous canopies have been less studied. Inspired by the method proposed by MüllerLinow et al. (2015) [34], we developed a portable and inexpensive measurement equipment called MultiDirectional Imager (MDI) that can capture stereo images of canopies, from which the leaf inclination angles can be computed. 3D line segment models were generated from stereo images taken from multiple directions, based on line segment detector and linebased 3D reconstruction method. Line segments can be used as image features to reconstruct the 3D scene [36]. Some methods have been proposed using the pure linebased StructurefromMotion (SfM) approach, which is challenging for pose estimation [37–39]. Camera pose information was used to overcome this problem by some subsequent Multiview Stereo (MVS) approaches (e.g., SURE [40]). However, MVS methods are typically accompanied by high computational costs. In this paper, we extracted 3D models with a low computational costs algorithm called Line3D++, developed by Hofer et al. (2015) [41].
In addition to the LAD, another factor that greatly affects the Beer’s lawbased LAI retrieval is the clumping effect, which is caused by the nonrandomly spatial distribution of leaves in the canopy [1, 42]. Few studies both correct the clumping effect and calculate the true Gfunction in LAI retrieval [43]. The clumping effect usually occurs at different spatial scales [44], including the large gapsinduced betweencrown clumping, crownshape clumping, and withincrown clumping due to the nonrandomness of foliage in the crown; these are difficult to be quantified. Thus, the clumping effect is not in the scope of this paper. Directly applying Beer’s law under the assumption of a randomly spatial distribution of leaves, we can obtain the effect of LAI. Given that measuring the angles of all the needles of a tree is challenging, we explored how to obtain the possible LADs and range of Gfunction of our fieldmeasured species. To expand the calculation of LAD and Gfunction on more trees with various species, we further used existing inventory data and 3D tree models. Our objective was to explore the possible effect of LAD on effective LAI retrieval for coniferous canopies, which we achieved by the following: (1)Developing a method for retrieving LAD for coniferous branches using MDI(2)Obtaining the possible range of the Gfunction for the canopies of the studied Spruce and Larch species(3)Calculating LADs and Gfunctions of other species based on existing inventory data and 3D tree models(4)Analyzing the possible theoretical impact of LAD on Beer’s lawbased LAI retrieval
2. Materials
Normally, the leaf area of a single needle is far smaller than that of a broadleaved tree. It is difficult to distinguish needles in field measurement, even with the use of highresolution digital cameras. Hence, we selected typical branch samples and measured the LADs of their needles in the laboratory. Numerous measurements were carried out to calculate the possible LADs and the range of Gfunction.
2.1. Field Measurement Using MultiDirectional Imager
2.1.1. MultiDirectional Imager
MDI (Figure 1) was developed based on our Multispectral Canopy Imager (MCI) [45]. The improvement included upgrading three highresolution digital cameras and adding a synchronous exposure device, which is useful to overcome the possible movements of the samples. MDI is a lowcost and portable instrument and can be used to obtain canopy images from multiple perspectives simultaneously without the help of any direction marker.
MDI consists of the following components: (i) three digital cameras (Canon 750D) with a suite of zoom lenses (EFS 18135 mm f/3.55.6 IS STM) per set. The CMOS size of each camera is , and the maximum resolution is per camera; (ii) a platform and supporting metal arms. The platform, which can be rotated from 0° to 360° in the azimuth direction and 0° to 90° in the zenith direction, is designed for controlling the camera’s orientation arbitrarily; (iii) a set of exposure controllers, which is used to control the simultaneous exposure of the three cameras; (iv) a tripod, which is used to support the platform and the cameras equipped on it.
Before the experiment, calibration is needed to obtain the interior orientation parameters of the cameras for the accurate measurement of the leaf inclination angles. The calibration of the MDI measuring system mainly consists of two parts: the calibration of the platform and that of the cameras. The platform was manufactured by a professional company and had been calibrated before use. The camera calibrations were carried out in the laboratory by using a checkerboard. Table 1, which can be used as a reference for selecting the camera focal length in the experiments, summarizes the theoretical maximum pixel size in an image when the actual observation distance is 10 m for two different camera resolutions ( and , respectively).

2.1.2. Study Sites and Samples
Field measurement campaigns of LADs were carried out from July 21 to July 26, 2016, in Chengde Saihanba National Forest Park (42.3°N, 117.2°E), including the local sampling of several conifer branches (Table 2): four Larch samples and three Spruce samples (see Section 4.2 for more details). There are some differences in measurement between Larch and Spruce samples: (i) for Larch sample I and II, the samples were fixed but the direction of the MDI observation changed; (ii) for Larch sample III and IV, the branches were rotated in different orientations but the MDI was fixed; (iii) for Spruce samples, both the samples and the direction of the MDI observation were fixed. This was done to test whether the direction of the MDI observation and the rotation of the branches have some influence on LAD and Gfunction (Details in section 4.3). In Table 2, the underlined sample ID denotes different observing directions of MDI around a specific sample; the corresponding number of photos taken by MDI are listed. In the experiment, the branch samples were scissored at different crown heights from a few representative trees. The branch samples were measured in the laboratory in one day after cutting to prevent the atrophy of the needles. The lengths of branches were measured, and the number of needle clusters and the number of needles in each cluster were counted. A white cloth was used as the background when taking photos in the laboratory to reduce the impact of background on object detection in photography.

2.2. 3D Tree Models and Inventory Data
What we did in section 2.1 is to obtain the possible LADs and range of Gfunction of Larch and Spruce through rotation of the branches. To investigate more species, we further used other kinds of data: (i) 3D tree models: 71 coniferous tree models of 13 species from OnyxTREE CONIFER library (http://www.onyxtree.com/conifer.html) with known true LADs (see Table 3); (ii) inventory data: 23 trees of Järvselja Scots Pine (https://ramibenchmark.jrc.ec.europa.eu/HTML/RAMIIV/EXPERIMENTS4/ACTUAL_CANOPIES/JARVSELJA_SUMMER_PINESTAND/JARVSELJA_SUMMER_PINESTAND.php) and Ofenpass Mountain Pine (https://ramibenchmark.jrc.ec.europa.eu/HTML/RAMIIV/EXPERIMENTS4/ACTUAL_CANOPIES/OFENPASS_WINTER_PINESTAND/OFENPASS_WINTER_PINESTAND.php) from RAMI IV. The measured LADs of these trees, one species measured in summer and the other in winter, are provided. More details about the inventory data can be found on the RAMI IV website. The number of trees, the average, and the range of tree heights for each species are summarized in Table 3. Front views of some of the studied 3D tree models from OnyxTREE are shown in Figure 2 as examples.

3. Method
3.1. Estimation of Leaf Inclination Angle Distribution
3.1.1. Detection and Reconstruction of Needles
The exact leaf inclination angles can be calculated from 3D leaf models. LSD [46] and Line3D++ [41] were used to generate 3D line segment models from stereo images of coniferous branches. To be specific, LSD was used to extract 2D lines from each image, and Line3D++ was used to develop 3D line models.
LSD is an efficient line segment detector that can be used to acquire accurate subpixel results [47]. It follows the method proposed by Burns et al. (1986) [48] and uses a contrarian validation approach proposed by Desolneux et al. (2000) [49] to control false detections. LSD algorithms can be roughly divided into three steps. First, the image is segmented into a series of small continuous regions with the same gradient angle. Then, the most approximate line segment is determined in each successive small region. Finally, a validation criterion is used to control the number of false detections.
Line3D++ is a novel algorithm using line segments as fundamental features to abstract 3D scenes, which consists of several steps. First, a large number of potential line correspondences are established using epipolar geometry constraints. The epipolar geometry can be computed through known camera poses, which can be obtained by running any StructurefromMotion (SfM) pipelines. Then, the most suitable correspondence is selected for each line segment. Hofer et al. (2017) [50] used a scoring formulation composed of angular similarity and positional similarity to distinguish between correct and incorrect matches. The final 3D line model is obtained by fusing 2D line correspondences through an arbitrary graphclustering algorithm.
Figure 3 shows an example of line segment detection and reconstruction using LSD and Line3D++ algorithms. Apart from some highly occluded regions, most of the needles are detected from the original images. The result (Figure 3(c)) shows that the majority of needles are reconstructed effectively.
3.1.2. Calculation of Leaf Inclination Angle Distribution
The data processing flow chart of the LAD retrieval is shown in Figure 3. First, MDI was used to obtain the multiangle stereo images of the branch samples (Figure 3(a)). Then, dense matching and linear feature extraction techniques were used to get the orientation of each needle automatically from the coaxial stereo image pairs obtained by MDI. Finally, the view zenith angle of MDI was used to calculate the final LAD (Figure 3(d)) based on the extracted inclination angles of the needles.
Through the above steps, we can obtain the distribution of the needles in the camera coordinate system. It has to be converted to the world coordinate system to obtain the real inclination angles of the needles. From the schematic diagram of the MDI coordinate system shown in Figure 1, represents a realworld local coordinate system. is the direction of the metal arms, and is the zenith direction. represents the camera space coordinate system of the intermediate camera. The angles of the platform in the vertical direction ( in Figure 1(b)) can be read directly from the graduated collar on the rotatable platform.
represents the real vertical direction, which is transformed from the world space to the camera space. The inclination of the needle in world space is defined as where indicates the normal vector of a needle in camera space, .
The final LAD is obtained by weighting the relative areas with different leaf inclination angles.
3.1.3. Validation of the Retrieved Leaf Inclination Angles
To validate the accuracy of our method, the inclination angles of needles on a coniferous tree branch were measured manually with a protractor and were assumed to be groundtruth. It is unrealistic to measure the inclination angle of each needle directly because the needles on a branch are dense and are sheltered from each other; hence, most of the needles were cut off from the branch, and only 16 remaining needles were measured manually (Figure 4(a)). Meanwhile, 3D needle models were reconstructed from the stereo images obtained from MDI, and the inclination angles were calculated automatically. Then, the coefficient of determination () and the root mean square error (RMSE) were calculated. The results are shown in Section 4.1.
3.2. GFunction Calculation
The definition of the Gfunction proposed by Stenberg (2006) [51] was used in this study: where is the view zenith angle, is the needle inclination angle, and is the azimuth divergence angle between the needle and the observation directions.
The measured leaf inclination angles are binned to the leaf inclination angle intervals ranging from 0 to with a step of 5°. Thus, the relationship between the discrete leaf inclination angle and the Gfunction is where is the total number of leaf inclination angle intervals and , is the center needle inclination angle of each interval, is the leaf area fraction of the needle inclination angle interval centered at , and is the value of the Gfunction when the needle inclination angle equals to .
3.3. Calculating the Possible Range of GFunction of the Canopy
For the field measured species including Larch and Spruce, it is impossible to measure every needle on a tree due to the huge number of needles and highly complex canopy background. Thus, what we can obtain is just the LAD of a certain branch, not that of a tree. However, we believe that the possible range of LAD and Gfunction of a tree can be obtained by rotating the branches to different orientations, assuming that a tree is made up of branches in different directions. Thus, the reconstructed 3D needle leaf models were rotated in many different orientations in the 3D space to represent possible orientations of needles and branches in the real canopy, and then, the corresponding changes in LAD and Gfunction were studied. To obtain the 3D needle leaf models with different orientations, we used CloudCompare (http://cloudcompare.org) to rotate each reconstructed 3D needle leaf models in fifteen different orientations. The LAD and Gfunction were subsequently computed using the rotated 3D needle leaf models. At last, all of the possible Gfunctions of the Larch and Spruce samples were calculated, which can be treated as possible situations in real canopies of these two species. The real Gfunctions of the canopies should be in the range of these simulated Gfunctions due to the fact that more branches in different orientations in reality tend to smooth the change of the Gfunction.
For other species from inventory data and 3D tree models, the range of Gfunction of the canopy is calculated directly from all the trees we collected. This is due to the fact that the LAD is a statistic of leaf inclination angles and is not related to spatial distribution of trees in the canopy.
3.4. Calculation of the Impact of Leaf Inclination Angle Distribution on LAI Retrieval
The traditional Beer’s law of light transmission through a turbid medium is widely used in indirect LAI measurements [3, 52], assuming that the medium of the vegetation is homogeneous: where denotes the gap fraction in the viewing zenith angle ; is the Gfunction, assuming of spherical distribution, . LAI is the effective LAI assuming a random distribution of leaves without considering the clumping effect.
We assumed the true Gfunction value in view zenith angle as , and the relative error of LAI caused by the assumption of spherical distribution is
From this equation, we infer that if we use the spherical distribution assumption, a 10% error in the Gfunction results in a 10% error in Beer’s lawbased LAI retrieval.
4. Results and Discussion
4.1. The Accuracy of Leaf Inclination Angles Retrieved from Stereo Photos
The leaf inclination angles of needles were measured manually as the groundtruth to validate the accuracy of our proposed method. Needles on a branch (Figure 4(a)) were manually measured, and the inclination angles were calculated automatically after they were detected. The enlarged details of the automatically extracted needles and the reconstructed point clouds are shown in Figure 4(b). As what is shown in Table 4, 14 of a total of 16 needles were extracted (except needles 11 and 15).

The comparison of manual measurements and algorithm extraction is shown in Figure 5. The algorithm performs well, with an of 0.98 and a RMSE of 1.63°. It shows that the proposed method can be used for accurate and rapid measurement of needle leaf inclination angle.
4.2. Leaf Inclination Angle Distribution and GFunction
Both the LADs and Gfunctions were estimated from different conifer samples: 4 Larch samples and 3 Spruce samples. Specifically, Larch sample I and II and the Spruce samples were fixed on the cantilever brackets at the time of data acquisition, and the position of MDI was moved. Sample III and IV on the cantilever brackets were manually rotated in different orientations. Photos around the samples were taken in different directions using MDI (Figure 6(a)), and at least 8 photos (details are listed in Table 2) with good quality were selected for line segment detection and 3D needle line segment reconstruction. The final LADs were calculated from the reconstructed needles (Figure 6(b)). Different trends were found in the results. For sample I and II, there were small differences in both the LADs and the mean leaf inclination angles (MLAs) measured in different directions, indicating that the observing direction of MDI has little impact on LAD. For sample III, the difference of MLAs between a31 and a32 was small, which was probably due to the relatively close orientations of them. However, there were relatively large differences in both the LADs and the MLAs before and after rotation for sample IV, and the maximum difference of MLA was about 10°.
Spruce samples were fixed on the cantilever brackets when data were acquired. Photos for different samples and the extracted LADs were shown in Figures 7(a1–a3) and 7(b1–b3), respectively. Results showed that the LADs might be different even with different samples of the same tree species, and the corresponding MLAs might have a large difference.
Figure 8 illustrates the Gfunction computed from abovementioned measurement data in this section. The Gfunction values calculated using the data obtained in different directions of the same conifer sample were close to each other, and the overall trend was consistent except for sample IV (black lines in Figure 8(a)). However, the Gfunction values of the three Spruce samples (Figure 8(b)) differ greatly.
4.3. Possible Range of the GFunction by Rotating the Reconstructed Needles
The main purpose of this section is to study the changes in LADs and Gfunctions with different needle orientations in the 3D space. We assume that the possible range of the Gfunction of the trees can be obtained by arbitrary rotations of the branches.
Larch sample IVa43 was used as an example to illustrate the experimental processes and results. First, each reconstructed 3D needle leaf model was rotated in 15 different orientations in the 3D space using CloudCompare (Figure 9). We saved the rotated 3D needle leaf models and computed LADs (Figure 10) and Gfunctions (Figure 11).
The other remaining samples were also processed according to the above steps. Several Gfunctions to represent the Larch and Spruce samples in the 3D space were acquired (Figure 12). Figure 12 contains the Gfunction calculated from both the original reconstructed 3D needle leaf models and those after rotation.
Figures 10–12 show that large differences exist in the corresponding LADs and Gfunctions when the 3D needle leaf models are rotated in different orientations in the 3D space, which also exists in MLA. The range of Gfunctions is largest when the view zenith angle is 0°, which is 0.13 for Larch samples and is up to 0.21 for Spruce samples. However, the maximum range of Gfunction is only 0.05 and 0.09 for Larch and Spruce samples, respectively, when the view zenith angles are between 40° and 70°, which is commonly used in field measurement campaigns. It can be inferred that changes in the values of the Gfunction are not significant in this range across the whole forest, because a real canopy has a greater number of randomly orientated branches and needles.
In addition, the overall variation range of the Gfunctions of Larch samples is smaller than that of Spruce samples. The reason may lie in the different growth structure of the two different conifer species.
4.4. GFunctions of Inventory Data and 3D Tree Models
The Gfunctions of trees from RAMI IV inventory data and 3D tree models with known LADs were also calculated. The range of Gfunction is the largest when the view zenith angle is 0°, which is 0.09 for the inventory data (Figure 13(a)) and is up to 0.31 for the studied 3D tree models (Figure 13(b)). However, the maximum range of the Gfunction is only 0.02 and 0.09 when the view zenith angles are in 40°~70° and is much smaller when the view zenith angles are in 45°~65°. In addition, the overall variation range of the Gfunction of the inventory data is smaller than that of the 3D tree models, which is because both the number of trees and tree species of the 3D tree models are significantly larger than those of the inventory data, i.e., the 3D tree models reflect more conditions of realistic coniferous trees.
4.5. The Effects of Leaf Inclination Angle Distribution on LAI Retrieval
Due to the inconvenience of in situ measurements of LAD of coniferous trees, the spherical distribution assumption is usually used in LAI retrieval. Previous studies showed that a spherical distribution is not a valid assumption for most broadleaf tree species [53]. However, relevant research in coniferous species is still lacking. In this section, whether the spherical distribution () is a valid assumption for coniferous canopies and the possible impact on LAI retrieval will be discussed. Specifically, we are interested in three kinds of view zenith angles: (i) 40°~70°, which is commonly used in field measurement campaigns; (ii) 0°, considering that the view zenith angle for spaceborne LiDAR instruments including Geoscience Laser Altimeter System (GLAS) and the Global Ecosystem Dynamics Investigation (GEDI) are all close to 0° (up to 1° for Geoscience Laser Altimeter System (GLAS) [54], and up to 6° for Global Ecosystem Dynamics Investigation (GEDI) [55]); (iii) 0°~15°, considering that the view zenith angle for airborne LiDAR is usually less than 15° and is widely used in LAI retrieval [56].
The relative errors of LAI retrieval under a constant Gfunction assumption of 0.5 were compared with those retrieved using the true Gfunction value for coniferous canopies (Figure 14). Results show that the relative errors of retrieved LAI depend on the view zenith angles and increase gradually with the deviation from 57.5°. When the view zenith angles are in 40°~70°, the relative errors of LAI are less than 15%. When the view zenith angle is at 0°, LAI could be underestimated by as much as 53%; when the view zenith angles range of 0°~15°, the relative errors of LAI are in the range of 53%~9%. However, the relative errors are less than 10% when the view zenith angles are in the range of 45°~65°.
Based on the above analysis, we recommend calculating LAD accurately in LAI retrieval, especially at nadir viewing. If the exact LAD measurement data is absent, we recommend using the remotely sensed data captured in the zenithal angle between 45° and 65°. In such cases, the possible LAI retrieval errors caused by fixing the Gfunction value as 0.5 is less than 10%.
4.6. Limitations and Suggestions for Future Research
It is difficult to obtain a general knowledge of the effect of the LAD on LAI retrieval for coniferous canopies just based on the field measurement. This is due to the fact that the field measurement is limited to a few branches of two species, since it is challenging to measure all the branches of a tree. Thus, the LAD of a tree is hard to be obtained. Two strategies were used in this study to overcome these problems: (i) rotating branches to obtain the possible range of the Gfunction and (ii) collecting as much data as possible to find a general impact of leaf inclination angles on LAI retrieval. Inevitably, there are some limitations of this work which should be further studied: (1)More branches can be measured in the future, considering that the selected branch samples may have some influence on the LAD and Gfunction estimation(2)More tree models can be collected to make the finding more general(3)More data of the same species can be collected in different seasons since the angles of leaves also change with the season. Even though both LAD data measured in summer and winter are included in the inventory data of RAMI IV, and the 71 tree models from a conifer database already contain diverse branches and leaves, which we believe can indirectly take the impact of season into consideration, a careful study on this issue is still an interesting topic(4)Both the LAD and clumping effect can be analyzed on the influence of true LAI retrieval. This study only focused on the impact of LAD on effective LAI. Previous studies showed that the clumping effect might cause a 30%~70% underestimation of true LAI [42, 57, 58]. Computer simulation, which can be used to generate diverse forest scenes with known gap fraction, LAD, and true LAI, might be a pathway to study the coupling effects of the clumping and the LAD to make sense of the total LAI retrieval error
5. Conclusions
Leaf inclination angle distribution (LAD) is an important factor affecting Beer’s lawbased LAI retrieval. Despite the high variability of LAD in reality, a spherical distribution assumption is often used for coniferous canopies due to the difficulties in measuring it. However, the error caused by this assumption has been seldom studied. We estimated the possible range of the Gfunction of two species by rotating the reconstructed needles of the branches measured by our developed equipment—MDI, which was an alternative to the manual measurement of the leaf inclination angle. To study more species, both 3D tree models and existing inventory data with known LADs were used. Results show that the range of the Gfunction is highly variable at different view zenith angles, indicating the necessity to measure LAD accurately to improve LAI retrieval. We found that the relative error of LAI with spherical distribution assumptions could be as large as 53% in the zenithal view direction, which is nonnegligible for LAI retrieval from airborne and spaceborne LiDAR. Besides, our results show that 45°~65° is a good range of the view zenith angle to reduce the impact of LAD on LAI retrieval (relative error less than 10%) for coniferous canopies under the assumption of spherical distribution when the fieldmeasured data is absent.
Data Availability
The data is freely available upon request.
Additional Points
General. Thank you for anonymous reviewers who gave constructive suggestions.
Conflicts of Interest
The authors declare no competing interests.
Authors’ Contributions
G.Y. proposed the method and wrote the paper with J.L., H.J., and X.M. conceptualized the method and revised the paper. F.L., J.Q., and R.H. contributed to the field measurement, and D.X. and G.Z. revised the paper.
Acknowledgments
This work was supported by the key program of the National Natural Science Foundation of China (NSFC) (Grant No. 42090013); Guangxi Innovative Development Grand Grant under the grant number: Guike AA18118038; and the China Scholarship Council, Grant No. 201906040055.
References
 G. Yan, R. Hu, J. Luo et al., “Review of indirect optical measurements of leaf area index: recent advances, challenges, and perspectives,” Agricultural and Forest Meteorology, vol. 265, pp. 390–411, 2019. View at: Publisher Site  Google Scholar
 G. Yan, R. Hu, Y. Wang et al., “Scale effect in indirect measurement of leaf area index,” IEEE Transactions on Geoscience and Remote Sensing, vol. 54, no. 6, pp. 3475–3484, 2016. View at: Publisher Site  Google Scholar
 C. T. De Wit, “Photosynthesis of leaf canopies. Agricultural Research Report No 663,” in Center for AgriCultural Publication and Documentation, pp. 1–57, Wageningen, The Netherlands, 1965. View at: Google Scholar
 J. Ross, The Radiation Regime and Architecture of Plant Stands, Springer, The Hague, Netherlands, 1981. View at: Publisher Site
 H. Utsugi, M. Araki, T. Kawasaki, and M. Ishizuka, “Vertical distributions of leaf area and inclination angle, and their relationship in a 46yearold Chamaecyparis obtusa stand,” Forest Ecology and Management, vol. 225, no. 13, pp. 104–112, 2006. View at: Publisher Site  Google Scholar
 K. Omasa, F. Hosoi, and A. Konishi, “3D Lidar imaging for detecting and understanding plant responses and canopy structure,” Journal of Experimental Botany, vol. 58, no. 4, pp. 881–898, 2006. View at: Publisher Site  Google Scholar
 F. Hosoi and K. Omasa, “Detecting seasonal change of broadleaved Woody canopy leaf area density profile using 3D portable Lidar imaging,” Functional Plant Biology, vol. 36, no. 11, pp. 998–1005, 2009. View at: Publisher Site  Google Scholar
 K. F. Huemmrich, “Simulations of seasonal and latitudinal variations in leaf inclination angle distribution: implications for remote sensing,” Advances in Remote Sensing, vol. 2, no. 2, pp. 93–101, 2013. View at: Publisher Site  Google Scholar
 K. Raabe, J. Pisek, O. Sonnentag, and K. Annuk, “Variations of leaf inclination angle distribution with height over the growing season and light exposure for eight broadleaf tree species,” Agricultural and Forest Meteorology, vol. 214215, pp. 2–11, 2015. View at: Publisher Site  Google Scholar
 A. R. G. Lang, R. E. McMurtrie, and M. L. Benson, “Validity of surface area indices of Pinus radiata estimated from transmittance of the sun’s beam,” Agricultural and Forest Meteorology, vol. 57, no. 13, pp. 157–170, 1991. View at: Publisher Site  Google Scholar
 J. M. Chen and T. A. Black, “Foliage area and architecture of plant canopies from sunfleck size distributions,” Agricultural and Forest Meteorology, vol. 60, no. 34, pp. 249–266, 1992. View at: Publisher Site  Google Scholar
 Y. Ryu, J. Verfaillie, C. Macfarlane et al., “Continuous observation of tree leaf area index at ecosystem scale using upwardpointing digital cameras,” Remote Sensing of Environment, vol. 126, pp. 116–125, 2012. View at: Publisher Site  Google Scholar
 R. Hu, G. Yan, X. Mu, and J. Luo, “Indirect measurement of leaf area index on the basis of path length distribution,” Remote Sensing of Environment, vol. 155, pp. 239–247, 2014. View at: Publisher Site  Google Scholar
 H. H. Neumann, G. Den Hartog, and R. H. Shaw, “Leaf area measurements based on hemispheric photographs and leaflitter collection in a deciduous forest during autumn leaffall,” Agricultural and Forest Meteorology, vol. 45, no. 34, pp. 325–345, 1989. View at: Publisher Site  Google Scholar
 S. G. Leblanc and J. M. Chen, “A practical scheme for correcting multiple scattering effects on optical LAI measurements,” Agricultural and Forest Meteorology, vol. 110, no. 2, pp. 125–139, 2001. View at: Publisher Site  Google Scholar
 D. L. B. Jupp, D. S. Culvenor, J. L. Lovell, G. J. Newnham, A. H. Strahler, and C. E. Woodcock, “Estimating forest LAI profiles and structural parameters using a groundbased laser called ‘Echidna®,” Tree Physiology, vol. 29, no. 2, pp. 171–181, 2008. View at: Publisher Site  Google Scholar
 D. Culvenor, G. Newnham, A. Mellor, N. Sims, and A. Haywood, “Automated insitu laser scanner for monitoring forest leaf area index,” Sensors, vol. 14, no. 8, pp. 14994–15008, 2014. View at: Publisher Site  Google Scholar
 A. R. G. Lang, “Leaf orientation of a cotton plant,” Agricultural Meteorology, vol. 11, pp. 37–51, 1973. View at: Publisher Site  Google Scholar
 H. Sinoquet and P. Rivet, “Measurement and visualization of the architecture of an adult tree based on a threedimensional digitising device,” Trees, vol. 11, no. 5, pp. 265–270, 1997. View at: Publisher Site  Google Scholar
 C. S. T. Daughtry, “Direct measurements of canopy structure,” Remote Sensing Reviews, vol. 5, no. 1, pp. 45–60, 1990. View at: Publisher Site  Google Scholar
 B. N. Bailey and W. F. Mahaffee, “Rapid measurement of the threedimensional distribution of leaf orientation and the leaf angle probability density function using terrestrial lidar scanning,” Remote Sensing of Environment, vol. 194, pp. 63–76, 2017. View at: Publisher Site  Google Scholar
 Y. Ryu, O. Sonnentag, T. Nilson et al., “How to quantify tree leaf area index in an open savanna ecosystem: a multiinstrument and multimodel approach,” Agricultural and Forest Meteorology, vol. 150, no. 1, pp. 63–76, 2010. View at: Publisher Site  Google Scholar
 J. Pisek, Y. Ryu, and K. Alikas, “Estimating leaf inclination and Gfunction from leveled digital camera photography in broadleaf canopies,” Trees, vol. 25, no. 5, pp. 919–924, 2011. View at: Publisher Site  Google Scholar
 X. Zou, M. Mõttus, P. Tammeorg et al., “Photographic measurement of leaf angles in field crops,” Agricultural and Forest Meteorology, vol. 184, pp. 137–146, 2014. View at: Publisher Site  Google Scholar
 B. E. McNeil, J. Pisek, H. Lepisk, and E. A. Flamenco, “Measuring leaf angle distribution in broadleaf canopies using UAVs,” Agricultural and Forest Meteorology, vol. 218219, pp. 204–208, 2016. View at: Publisher Site  Google Scholar
 X. Li, J. Wang, and Y. Xiang, “Measuring leaf angle distribution and leaf area index of conifer canopies simultaneously on wide angle image,” Journal of Remote Sensing, vol. 1, pp. 54–61, 1997. View at: Google Scholar
 C. J. Kucharik, J. M. Norman, and S. T. Gower, “Measurements of leaf orientation, light distribution and sunlit leaf area in a boreal aspen forest,” Agricultural and Forest Meteorology, vol. 91, no. 12, pp. 127–148, 1998. View at: Publisher Site  Google Scholar
 F. Hosoi, K. Nakabayashi, and K. Omasa, “3D modeling of tomato canopies using a highresolution portable scanning lidar for extracting structural information,” Sensors, vol. 11, no. 2, pp. 2166–2174, 2011. View at: Publisher Site  Google Scholar
 K. Zhao, M. García, S. Liu et al., “Terrestrial lidar remote sensing of forests: maximum likelihood estimates of canopy profile, leaf area index, and leaf angle distribution,” Agricultural and Forest Meteorology, vol. 209210, pp. 100–113, 2015. View at: Publisher Site  Google Scholar
 F. Hosoi and K. Omasa, “Estimating leaf inclination angle distribution of broadleaved trees in each part of the canopies by a highresolution portable scanning lidar,” Journal of Agricultural Meteorology, vol. 71, no. 2, pp. 136–141, 2015. View at: Publisher Site  Google Scholar
 G. Zheng and L. M. Moskal, “Leaf orientation retrieval from terrestrial laser scanning (TLS) data,” IEEE Transactions on Geoscience and Remote Sensing, vol. 50, no. 10, pp. 3970–3979, 2012. View at: Publisher Site  Google Scholar
 L. Ma, G. Zheng, J. U. H. Eitel, T. S. Magney, and L. M. Moskal, “Retrieving forest canopy extinction coefficient from terrestrial and airborne lidar,” Agricultural and Forest Meteorology, vol. 236, pp. 1–21, 2017. View at: Publisher Site  Google Scholar
 H. Wang, W. Zhang, G. Zhou, G. Yan, and N. Clinton, “Imagebased 3D corn reconstruction for retrieval of geometrical structural parameters,” International Journal of Remote Sensing, vol. 30, no. 20, pp. 5505–5513, 2009. View at: Publisher Site  Google Scholar
 M. MüllerLinow, F. PintoEspinosa, H. Scharr, and U. Rascher, “The leaf angle distribution of natural plant populations: assessing the canopy with a novel software tool,” Plant Methods, vol. 11, no. 1, p. 11, 2015. View at: Publisher Site  Google Scholar
 J. Qi, D. Xie, L. Li, W. Zhang, X. Mu, and G. Yan, “Estimating leaf angle distribution from smartphone photographs,” IEEE Geoscience and Remote Sensing Letters, vol. 16, no. 8, pp. 1190–1194, 2019. View at: Publisher Site  Google Scholar
 C. Baillard, C. Schmid, A. Zisserman, and A. Fitzgibbon, “Automatic line matching and 3D reconstruction of buildings from multiple views,” in Paper presented at the ISPRS Conference on Automatic Extraction of GIS Objects from Digital Imagery, Munich, Germany, 1999. View at: Google Scholar
 A. Bartoli and P. Sturm, “Structurefrommotion using lines: representation, triangulation, and bundle adjustment,” Computer Vision and Image Understanding, vol. 100, no. 3, pp. 416–441, 2005. View at: Publisher Site  Google Scholar
 H. Bay, A. Ess, A. Neubeck, and L. Van Gool, “3D from line segments in two poorlytextured, uncalibrated images,” in Third International Symposium on 3D Data Processing, Visualization, and Transmission (3DPVT'06), pp. 496–503, Chapel Hill, NC, USA, 2006. View at: Publisher Site  Google Scholar
 L. Zhang and R. Koch, “Structure and motion from line correspondences: representation, projection, initialization and sparse bundle adjustment,” Journal of Visual Communication and Image Representation, vol. 25, no. 5, pp. 904–915, 2014. View at: Publisher Site  Google Scholar
 M. Rothermel, K. Wenzel, D. Fritsch, and N. Haala, “SURE: photogrammetric surface reconstruction from imagery,” in Paper presented at the Proceedings LC3D Workshop, pp. 1–9, Berlin, Germany, 2012. View at: Google Scholar
 M. Hofer, M. Maurer, and H. Bischof, “Line3D: Efficient 3D Scene Abstraction for the Built Environment,” in Paper presented at the Pattern Recognition, pp. 237–248, Springer, 2015. View at: Publisher Site  Google Scholar
 M. Weiss, F. Baret, G. J. Smith, I. Jonckheere, and P. Coppin, “Review of methods for in situ leaf area index (LAI) determination: Part II. Estimation of LAI, errors and sampling,” Agricultural and Forest Meteorology, vol. 121, no. 12, pp. 37–53, 2004. View at: Publisher Site  Google Scholar
 X. Zhu, A. K. Skidmore, T. Wang et al., “Improving leaf area index (LAI) estimation by correcting for clumping and woody effects using terrestrial laser scanning,” Agricultural and Forest Meteorology, vol. 263, pp. 276–286, 2018. View at: Publisher Site  Google Scholar
 L. Ma, G. Zheng, X. Wang, S. Li, Y. Lin, and W. Ju, “Retrieving forest canopy clumping index using terrestrial laser scanning data,” Remote Sensing of Environment, vol. 210, pp. 452–472, 2018. View at: Publisher Site  Google Scholar
 J. Zou, G. Yan, L. Zhu, and W. Zhang, “Woodytototal area ratio determination with a multispectral canopy imager,” Tree Physiology, vol. 29, no. 8, pp. 1069–1080, 2009. View at: Publisher Site  Google Scholar
 R. G. von Gioi, J. Jakubowicz, J. M. Morel, and G. Randall, “LSD: a fast line segment detector with a false detection control,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 32, no. 4, pp. 722–732, 2010. View at: Publisher Site  Google Scholar
 R. G. von Gioi, J. Jakubowicz, J.M. Morel, and G. Randall, “LSD: a line segment detector,” Image Processing On Line, vol. 2, pp. 35–55, 2012. View at: Publisher Site  Google Scholar
 J. B. Burns, A. R. Hanson, and E. M. Riseman, “Extracting straight lines,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. PAMI8, no. 4, pp. 425–455, 1986. View at: Publisher Site  Google Scholar
 A. Desolneux, L. Moisan, and J.M. Morel, “Meaningful alignments,” International Journal of Computer Vision, vol. 40, no. 1, pp. 7–23, 2000. View at: Publisher Site  Google Scholar
 M. Hofer, M. Maurer, and H. Bischof, “Efficient 3D scene abstraction using line segments,” Computer Vision and Image Understanding, vol. 157, pp. 167–178, 2017. View at: Publisher Site  Google Scholar
 P. Stenberg, “A note on the _G_ function for needle leaf canopies,” Agricultural and Forest Meteorology, vol. 136, no. 12, pp. 76–79, 2006. View at: Publisher Site  Google Scholar
 J. L. Monteith, “Light distribution and photosynthesis in field crops,” Annals of Botany, vol. 29, no. 1, pp. 17–37, 1965. View at: Publisher Site  Google Scholar
 J. Pisek, O. Sonnentag, A. D. Richardson, and M. Mõttus, “Is the spherical leaf inclination angle distribution a valid assumption for temperate and boreal broadleaf tree species?” Agricultural and Forest Meteorology, vol. 169, pp. 186–194, 2013. View at: Publisher Site  Google Scholar
 A. C. Brenner, H. J. Zwally, C. R. Bentley et al., “Derivation of range and range distributions from laser pulse waveform analysis for surface elevations,” Roughness, Slope, and Vegetation Heights, Geoscience Laser Altimeter System Algorithm Theoretical Basis Document, Version 4, 2003, http://www.csr.utexas.edu/glas/pdf/Atbd_20031224.pdf. View at: Google Scholar
 H. Tang and J. Armston, Algorithm Theoretical Basis Document (ATBD) for GEDI L2B Footprint Canopy Cover and Vertical Profile Metrics, Goddard Space Flight Center, Greenbelt, MD, USA, 2020, https://lpdaac.usgs.gov/documents/588/GEDI_FCCVPM_ATBD_v1.0.pdf.
 R. Hu, G. Yan, F. Nerry et al., “Using airborne laser scanner and path length distribution model to quantify clumping effect and estimate leaf area index,” IEEE Transactions on Geoscience and Remote Sensing, vol. 56, no. 6, pp. 3196–3209, 2018. View at: Publisher Site  Google Scholar
 J. M. Chen and J. Cihlar, “Plant canopy gapsize analysis theory for improving optical measurements of leafarea index,” Applied Optics, vol. 34, no. 27, pp. 6211–6222, 1995. View at: Publisher Site  Google Scholar
 P. Stenberg, “Correcting LAI2000 estimates for the clumping of needles in shoots of conifers,” Agricultural and Forest Meteorology, vol. 79, no. 12, pp. 1–8, 1996. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2021 Guangjian Yan et al. Exclusive Licensee Aerospace Information Research Institute, Chinese Academy of Sciences. Distributed under a Creative Commons Attribution License (CC BY 4.0).