Stand Canopy Closure Estimation in Planted Forests Using a Geometric-Optical Model Based on Remote Sensing

: Canopy closure, which is the ratio of the vertical projection area of the crowns to the area of forest land, can indicate the growth and tending situation of a forest and is of great signiﬁcance for forest management planning. In this study, a geometric-optical model (GOST model) was used to simulate the canopy gap fraction of a forest. Then, a canopy closure estimation method using the gap fraction was discussed. In this study, three typical planted forest farms (the Mengjiagang (MJG), Gaofeng (GF), and Wangyedian (WYD) forest farms) containing the most commonly planted tree species in the north and south regions of China were selected, and ﬁeld measurements were executed. The results show that the gap fraction ( P vg-c ) had a higher correlation with the average projected area of the tree crowns, and the relationship was an exponential function, with R 2 and RMSE values of 0.5619 and 0.0723, respectively. Finally, the applicability and accuracy of this method were evaluated using line transects, and a ﬁsheye camera measured the canopy closure. The accuracy of the canopy closure estimated by the P vg-c was 86.69%. This research can provide a reference for canopy closure estimation using a geometric-optical model.


Introduction
Globally, planted forests are an important type of forest. According to a number of studies, planted forest areas have continued to increase due to industrial demands for wood shifting from natural forests to planted forests [1,2]. Planted forests reduce harvesting from natural forests by 26%, and they have significant ecological benefits [3]. In order to improve the management level of China's planted forests, it is essential to plan the afforestation process and to scientifically arrange the forest management strategy. At the same time, it is also important to monitor and evaluate the resources in planted forests precisely and meticulously [4].
Canopy closure, which is the ratio of the vertical projection area of tree crowns to the ground area [5][6][7][8], plays a very important role in the state of the forest ecosystem as well as in environmental evaluation, and it is widely used in forestry evaluations [9]. Canopy closure is an important index that reflects the spatial structure of an ecosystem and the tree stand density, and it is an important investigation indicator for planted forests [10,11]. Accurate canopy closure measurements and estimations provide an important reference that can be used to evaluate the quality of a plantation.
Canopy closure can be measured using a wide variety of ground-based techniques. These ground-based methods mainly include ocular estimates, hemispherical photography, transects, sample points, the line intercept method, canopy projections, the visual observation method, and canopy instrument analysis, among others [12]. Among these methods, transects and sample points have a higher measurement accuracy than other ground-based methods [13,14]. Hemispherical photography is also a commonly used method that is used to estimate the canopy closure [13]. The development of digital cameras and their increasing availability has made hemispherical photography more widely available for forest inventory evaluations [15]. However, the cost and required resources for hemispherical photography still preclude many forest managers from using it as a monitoring tool [11]. In addition, canopy closure is underestimated by using the instruments, such as fisheye cameras, due to gaps in the canopy [13]. Macfarlane et al. concluded that future investigations of fisheye camera methods should concentrate on obtaining an accurate gap fraction distribution method that can separate the effects of foliage angle distribution from those of foliage clumping [16]. The transect method is considered to be the most reliable method for investigating canopy closure, and it can directly verify the canopy closure estimated by remote sensing images [17]. However, even though ground-based methods always have better accuracy, human and material resource inputs need to be considered. Additionally, it is difficult to obtain the canopy closure distribution at a regional or larger scale [18].
As an efficient and low-cost data resource, remote sensing is regarded as one of the most effective ways to estimate canopy closure in regional or large areas [19][20][21][22]. Aerial platforms are used to obtain photos of the forest region during the early stages, and originally, these aerial photos were obtained to estimate forest information such as the canopy closure, LAI, and forest volume [23][24][25]. With the development of laser ranging technology, light detection and ranging (LiDAR) data have been used to estimate forest structural parameters, such as tree height as well as tree crown and single tree properties [26,27]. Placing LiDAR sensors on the aerial platforms is a good option for tree canopy closure estimation [28]. With the emergence of various types of LiDAR sensors on the market, more and more research on canopy closure estimation using LIDAR data has been published [29][30][31]. Additionally, the high-precision canopy height model (CHM) is extracted using LIDAR data and image segmentation, and canopy width extraction technology is used to estimate canopy closure [32,33]. Bode et al. proposed a canopy closure estimation approach using a light penetration index (LPI) based on airborne LiDAR data at watershed scales [34]. However, huge economic costs must be considered, especially when aerial platforms are being used to estimate the canopy closure in large areas [35]. Additionally, the costs make it very difficult to determine the forest inventory at a high frequency and in a large-scale region [36]. Compared to aerial platforms, satellite remote sensing platforms have advantages, such as good real-time, large-scale, low-cost, accurate, and multiple types of sensors. Satellite remote sensing images have been widely used to determine the forest inventory at both local and large scales [37][38][39].
To estimate canopy closure in a large area, statistical or physical models are used to estimate the forest canopy closure. Statistical models represent economical and efficient canopy closure estimation methods that establish the relationship between remote sensing variables and measured field data [40][41][42].
Among the various statistical models, multiple linear regression models are usually used to establish the relationship between the canopy closure and the remote sensing variables obtained from remote sensing images, such as spectral information, the vegetation index, and texture information [43][44][45]. However, the method does not fully consider the influence of the spatial resolution of remote sensing images on the accuracy of canopy closure estimations [46,47]. Additionally, the effect of mixed pixels on the accuracy of canopy closure estimations is ignored in this method. The pixel binary model considers mixed pixels, it cannot separate arbors and shrubs, and it is also affected by the background soil type, meaning that it is difficult for this method to provide accurate stand canopy closure information [48][49][50][51].
Physical models are other types of canopy closure estimation models that are based on the energy transmission procedures of remote sensing. Typical physical models include the geometric-optical model and the radiative transfer model [52]. The Li-Strahler geometric-optical model is suitable for canopy closure estimations in flat study areas [53,54]. Additionally, the core of this method is its ability to accurately extract proportions of four scene components-sunlit foliage, sunlit ground, shaded foliage, and shaded ground-in a pixel [52,55]. It is suitable for simulating canopy reflectance by considering forest distribution and structure. The PROSAIL model is one of the most famous radiative transfer models used for canopy reflection simulation and combines the PROSPECT leaf model and SAIL canopy model (arbitrary oblique leaf-scattering) [56][57][58]. The PROSAIL model is also a commonly used radiation transmission model for canopy closure estimation [59,60]. The geometric-optical model and the radiation transmission model are more complex and require more parameters, but parameter sensitivity analysis can reduce the complexity of these models. Additionally, canopy closure estimation models that are established with fewer parameters are less affected by the region and measured dataset, resulting in them having better robustness [61].
In summary, there is a lot of uncertainty and unknowns when estimating canopy closure using remote sensing techniques. In this study, we considered the canopy from two perspectives-an opaque canopy body or with gaps in the tree crown. Then, the Geometric-Optical Model for Sloping Terrains (GOST) model, which is a geometric-optical model that considers terrain effects on the input parameters, was used to simulate the canopy gaps from these two perspectives. Finally, the stand-scale canopy closure was estimated in planted forests. The estimated results were tested using the measured canopy closure results obtained from transects and fisheye camera photos to answer the following questions: (1) For remote sensing data with high spatial resolution, which parameter is appropriate to estimate the stand-scale canopy closure with the highest accuracy? (2) For remote sensing data with low or medium spatial resolution, how can a geometric-optical model be used to establish a better robustness canopy closure model while considering the mixed pixel problem? (3) How can we evaluate the feasibility and accuracy of the canopy closure method presented in this study in a planted forest?

Study Area
In China, the main forest types include evergreen broadleaf planted forests, deciduous broad-leaved planted forests, evergreen coniferous planted forests, and deciduous needleleaf planted forests. To enhance the accuracy of the proposed method and to achieve better robustness, research regions that included most of the forest types were selected from around China. The tree species that were planted in these forest farms were also the most commonly planted species in China. The Wangyedian Forest Farm, Gaofeng Forest Farm, and Mengjiagang Forest Farm were chosen for this study because they represent the typical plantation types in the north and south regions of China ( Figure 1).
The Mengjiagang Forest Farm (MJG) is located in Heilongjiang Province, which is located in the northeast part of China. It is located at 130 • 32 -130 • 52 E, 46 • 20 -46 • 30 N. The relative height of the mountain on which it is located is 168 to 575 m, and the slope is 10-20 • . The total area of the forest farm is 1.7 × 10 4 km 2 , and the plantation accounts for 76.7% of the total area and mainly includes larch (Larix gmelinii (Rupr.) Kuzen.), mongolica (Pinus sylvestris var. mongolica Litv.), and Korean pine (Pinus koraiensis Sieb. et Zucc.), with these tree species accounting for about 80% of the plantation population.
The Gaofeng Forest Farm (GF) is located in Guangxi Province, which is located in the south of China. It is located at 108 • 08 -108 •  area is about 2.5 × 10 4 km 2 , and the forest coverage rate is over 80%. The elevation is between 800 m to 1890 m. The relative height of the mountain is around 200 to 400 m and has a slope of 15-35°. The area of the forest farm is composed of 47% planted forest and 53% natural forest, and the main tree species include oil pine (Pinus tabuliformis Carrière) and larch (Larix gmelinii (Rupr.) Kuzen.). Oil pine is present in a proportion of 51%, and the proportion of larch is 47%. The natural forests in this area are mainly composed of white birch (Betula platyphylla), black birch (Betula dahurica), mountain apricot (Armeniaca sibirica), aspen (Populus davidiana), elm (Ulmus pumila), hazelnut (Corylus heterophylla), and other tree species.

Field Data
A total of 102 plots were set up in the WYD, GF, and MJG. When setting up the plots, it was necessary to ensure that the plots were fully representative of the stand, that they were not scattered across different forest types, and that they were evenly distributed at different levels of the slope and aspect in each forest farm. A total of 30 plots were set up in the WYD with plot areas of 25 × 25 m, of which 18 plots were oil pine forests and 12 were larch forests. Additionally, there were 43 plots in the GF with plot areas of 20 × 20 m, with eucalyptus forests comprising most of the plots. There were 29 plots in the MJG

Field Data
A total of 102 plots were set up in the WYD, GF, and MJG. When setting up the plots, it was necessary to ensure that the plots were fully representative of the stand, that they were not scattered across different forest types, and that they were evenly distributed at different levels of the slope and aspect in each forest farm. A total of 30 plots were set up in the WYD with plot areas of 25 × 25 m, of which 18 plots were oil pine forests and 12 were larch forests. Additionally, there were 43 plots in the GF with plot areas of 20 × 20 m, with eucalyptus forests comprising most of the plots. There were 29 plots in the MJG with plot Remote Sens. 2022, 14, 1983 5 of 22 areas of 30 × 20 m. All of the plots were larch forests. The data were normalized to a minimum size of 20 × 20 m to ensure the consistency of the data.
The measured forest parameters included diameter at breast height, tree height, stem height, crown width, and tree species. All of the trees with a diameter at breast height above 5 cm in the plot were measured and recorded in one plot. Then, the GPS coordinates of each plot were recorded. The statistical information of the three research areas can be found in Table 1. The canopy closure, leaf area index (LAI), and clumping index of each plot were measured at the same time. The canopy closure was measured using transects and fisheye camera photos. Two transect lines were laid along the diagonals of the plot, and the vertical projection lengths of the crowns along each transect were recorded. The average ratio of the total projection lengths of the crowns along the two diagonals to the diagonal length was the canopy closure.
The LAI was measured using fisheye camera photos. All of the photos were taken with a Nikon Coolpix 4500. The photos were taken at the four corners and in the center of each plot. The locations in which the images were collected can be found in Figure 2a. The images were analyzed using digital hemispherical photography (DHP). The edges of the photos with large amounts of distortion were removed and processed via binarization. Then, the canopy pixel value was recorded as 1, the gap was recorded as 0, and the ratio of the number of pixels with a value of 1 to the total number of pixels was the canopy closure. Finally, the average canopy closure of the five images was the canopy closure of the plot obtained using the fisheye camera method. At the same time, the effective LAI was also able to be automatically calculated using DHP. Details of this method can be found in [62]. with plot areas of 30 × 20 m. All of the plots were larch forests. The data were normalized to a minimum size of 20 × 20 m to ensure the consistency of the data. The measured forest parameters included diameter at breast height, tree height, stem height, crown width, and tree species. All of the trees with a diameter at breast height above 5 cm in the plot were measured and recorded in one plot. Then, the GPS coordinates of each plot were recorded. The statistical information of the three research areas can be found in Table 1. The canopy closure, leaf area index (LAI), and clumping index of each plot were measured at the same time. The canopy closure was measured using transects and fisheye camera photos. Two transect lines were laid along the diagonals of the plot, and the vertical projection lengths of the crowns along each transect were recorded. The average ratio of the total projection lengths of the crowns along the two diagonals to the diagonal length was the canopy closure.
The LAI was measured using fisheye camera photos. All of the photos were taken with a Nikon Coolpix 4500. The photos were taken at the four corners and in the center of each plot. The locations in which the images were collected can be found in Figure 2a. The images were analyzed using digital hemispherical photography (DHP). The edges of the photos with large amounts of distortion were removed and processed via binarization. Then, the canopy pixel value was recorded as 1, the gap was recorded as 0, and the ratio of the number of pixels with a value of 1 to the total number of pixels was the canopy closure. Finally, the average canopy closure of the five images was the canopy closure of the plot obtained using the fisheye camera method. At the same time, the effective LAI was also able to be automatically calculated using DHP. Details of this method can be found in [62]. The clumping index of each plot was measured by using a TRAC instrument. Two line transects that were 20 m in length were measured. Then, the TRAC-based PPFD gradient values along the transects perpendicular to the incident directions of the solar beams The clumping index of each plot was measured by using a TRAC instrument. Two line transects that were 20 m in length were measured. Then, the TRAC-based PPFD gradient values along the transects perpendicular to the incident directions of the solar beams were collected. In the end, the clumping index of each plot was calculated using TRACWin software. Details can be found in [63].

The GOST Model and Canopy Gap Fraction Simulation
The GOST model is a geometric-optical model for sloping terrains developed based on the four-scale model [64]. The four-scale model is one of the most popular geometric-optical models and can be used to simulate the bidirectional reflectance distribution function of forest canopies on flat surfaces [65]. By considering the structure of the canopy at four scales, viz., tree groups, tree crowns, branches, and shoots, the bidirectional reflectance characteristics of the forest canopies can be simulated. The four-scale model defines how the canopy reflectance in one pixel is a linear combination of the signals from four components: the sunlit and shaded foliage and the sunlit and shaded backgrounds. The total canopy reflectance is as follows: where R T is the reflectivity of the sunlit foliage; R G is the reflectivity of the sunlit background; R ZT is the reflectivity of the shaded foliage, R ZG is the reflectivity of the shaded background; P T , P G , Z T , and Z G are the sensor-viewing probabilities of the four scene components, respectively ( Figure 3).
were collected. In the end, the clumping index of each plot was calculated TRACWin software. Details can be found in [63].

The GOST Model and Canopy Gap Fraction Simulation
The GOST model is a geometric-optical model for sloping terrains develope on the four-scale model [64]. The four-scale model is one of the most popular geo optical models and can be used to simulate the bidirectional reflectance distributio tion of forest canopies on flat surfaces [65]. By considering the structure of the ca four scales, viz., tree groups, tree crowns, branches, and shoots, the bidirectiona tance characteristics of the forest canopies can be simulated. The four-scale model how the canopy reflectance in one pixel is a linear combination of the signals fro components: the sunlit and shaded foliage and the sunlit and shaded backgroun total canopy reflectance is as follows: where is the reflectivity of the sunlit foliage; is the reflectivity of the sunl ground; is the reflectivity of the shaded foliage, is the reflectivity of the background; , , , and are the sensor-viewing probabilities of the fou components, respectively ( Figure 3). In the four-scale model, the probability of seeing the ground represents the that can be seen between tree crowns. It is a function of the projected tree crown a the spatial distribution characteristics of trees in a quadrat. It can be calculated as based on the method in [65]: where _ is the probability of seeing the ground, i is the number of trees in a with the area A, and is the ground surface not seen by the viewer because of o ( ) is the probability of having i trees in A, and it is determined by the Neyma bution. In this equation, the tree crowns are assumed to be opaque, and the gaps the tree crowns are not considered. In the four-scale model, the probability of seeing the ground represents the ground that can be seen between tree crowns. It is a function of the projected tree crown area and the spatial distribution characteristics of trees in a quadrat. It can be calculated as follows, based on the method in [65]: where P vg-c is the probability of seeing the ground, i is the number of trees in a region with the area A, and V g is the ground surface not seen by the viewer because of one tree. P N (i) is the probability of having i trees in A, and it is determined by the Neyman distribution. In this equation, the tree crowns are assumed to be opaque, and the gaps within the tree crowns are not considered.
If the gaps in the crowns and the overlaps in the crowns are considered, then Equation (2) can be written as follows: This is the gap probability inside the trees.
The gap probability P gap (θ v ). is calculated based on the method developed by Li and Strahler [53], but the foliage clumping effect is considered in the four-scale model. G(θ v ). is the function of the foliage angle distribution G(θ v ) = 0.5 [65] in this research. L 0 is the LAI. Ω E . is the clumping index of the shoots within the tree crowns. γ E is the ratio of the needle to shoot area.
In Equation (2), P tj V g is the probability of having j trees intercepting the view line. Additionally, it can be calculated using a negative binomial function.
In the case of j = 0, P t0 is equal to Equation (2). The canopy gap fraction determines the contribution of the under surface to the reflectance measured above the canopy. Additionally, the canopy gap fraction can be calculated using the tree crown projection made on the ground surface. It is based on the method described above for calculating the shadow area. After calculating the shadow area projected by a single tree onto the ground, the sunlit crown proportion seen by the viewer can be computed using the total surface area of the tree visible to the viewer projected to a plane perpendicular to the view line. It is not difficult to determine that the canopy gap fraction is highly related to canopy closure according to the method in the four-scale model.
In addition, topography is an important factor that has a serious effect on the bidirectional reflectance distribution function of forest canopies. Therefore, in order to make the four-scale model more suitable to simulate the bidirectional reflectance distribution function of the forest canopies on slopes, the GOST model, which is a geometric-optical model that considers the oblique topographical factors, was used in this study [64].
The GOST model was used to simulate the canopy reflectance for two different gap fractions. In one case, the gap fraction (P vg-c ) was calculated by Formula (2) during the assumption of the opaque canopy bodies. In the other case, the gap fraction (P vg ) was calculated using Formula (3), and gaps and overlaps were observed in the crowns. The output of the GOST model was the canopy reflectance and probability of seeing the four scene components under different view angles (Equation (1)). The four components were the sunlit foliage, the sunlit ground, the shaded foliage, and the shaded ground, and the mixed pixel decomposition problem was able to be solved as well. A database of the reflectance under various gap fractions was established, allowing the relationships between the canopy gap fraction and the stand canopy closure to be discussed. When the canopy was an opaque body, the estimated canopy closure was the true stand canopy closure, and the canopy closure measured by the traditional transects was used to verify the estimated results. When the gaps in the crowns were considered, the estimated canopy closure was compared to those measured results obtained via fisheye camera photos. This study discusses the feasibility of estimating the canopy closure using the GOST model and the spatial resolution effects of the remote sensing images.

Canopy Closure Estimation Based on the GOST Model
The GOST model was selected and used to simulate the canopy reflectance for a complex canopy structure. When the canopy reflectance was able to be estimated accurately, it was found to be closely related to the gap fraction function. This was because the canopy gap fraction was used to describe the light passing through the canopy. As such, the canopy gap fraction can describe canopy closure. A higher canopy gap fraction means less canopy closure, and vice versa. Additionally, it is also the reason why the canopy gap fraction can be used to estimate the canopy closure.
Two gap fractions were included in the GOST model. In one case, the canopy was assumed to be an opaque canopy body. Additionally, the gaps within the tree crown were not considered. Since canopy closure does not consider the gaps in the canopy, if the viewing direction was vertical, then the gap fraction of the viewing canopy gap fraction is the percentage of light passing through the canopy and projected to the ground surface. Additionally, the stand canopy closure should be 1 − P vg-c .
where CC 1 is the canopy closure of the opaque canopy bodies and P vg-c is the gap fraction of the opaque canopy bodies.
In the other case, the gap fraction was assumed to be a tree crown with gaps. Generally, it was assumed that there were gaps in the canopy for the canopy closure values that were obtained by the fisheye camera images or estimated from the remote sensing data; therefore, the estimated canopy closure should be 1 − P vg . As such, the canopy closure could be calculated as follows: where CC 2 and P vg were the canopy closure and gap fraction when there were tree crown gaps, respectively. The input parameters of the GOST model were determined to be {x 1 , x 2 , · · · x n }, and then the P vg-c and P vg could be simulated under various input parameters using the GOST model, and the database for the P vg-c and P vg and the input parameters was established. The relationships among the P vg-c , P vg , and the inputs were able to be established using a statistical method. P vg-c = F(x 1 , x 2 , · · · x n ) (9) where F(x 1 , x 2 , · · · x n ) and G(x 1 , x 2 , · · · x n ) were the functions of the P vg-c and P vg related inputs, respectively. Once the relationship between the gap fraction and the inputs was established, then the gap fraction could be estimated, and the canopy closure was able to be inverted based on Equations (7) and (8).

Sensitivity Analysis of the GOST Model Parameters
Geometric-optical models have a large number of input variables, but we only need to know the key input variables. It is crucial to carry out a sensitivity analysis on the parameters for simulation in a complex model, especially for a model with multiple input parameters. The inputs of the GOST model included plot parameters (size of plot, number of trees in plot, LAI, slope, aspect, solar zenith angle, solar azimuth angle, view zenith angle, view azimuth angle), tree structural parameters (radius of the crowns, stem height, crown height, half apex angle, clumping index for shoots in the crown), and spectral parameters (leaf reflectivity, leaf transmittance, and ground reflectance). Therefore, the sensitivity analysis was an important process in this study. The function (Equation (11)) of the sensitivity analysis was used to determine the sensitivity of the inputs affecting the gap Remote Sens. 2022, 14,1983 9 of 22 fraction of the canopy. It was very helpful to improve the accuracy and robustness of the model for the canopy closure estimation [66,67].
where P j 0 is the P vg-c , and P vg is calculated based on the measured parameters (reference). P j pert is the P vg-c and P vg when the parameter was disturbed, and F is the sensitivity of the input parameters [68].
The sensitivity of the tree structure parameters, leaf reflectivity, leaf transmittance, and ground reflectivity was analyzed in this study. Additionally, the topography effects were corrected in the GOST model, so there was no need for a slope sensitivity analysis.
According to the results of the sensitivity analysis, the input parameter dataset, which was sensitive to canopy closure, was determined as {x 1 , x 2 , · · · x n }. The insensitive inputs were set as the average values of the measurements in the plot. Then, the P vg-c and P vg were able to be simulated under various input parameters using the GOST model, and the database between the P vg-c and P vg and the input parameters was established. The relationship between the P vg-c , P vg , and the inputs was able to be established, and the canopy closure could then be inverted based on Equations (7) and (8).

Validation
In this study, the R 2 and RMSE were selected to evaluate the precision of the canopy closure estimation model. These were calculated using Equations (12) and (13).
where y i was the measured value of the sample i,ŷ i was the predicted value of the sample i, y was the average of all the samples, and n was the number of samples.

Results of Sensitivity Analysis of the Parameters in the GOST Model
Since the canopy closure was the ratio of the vertical projection area of the crowns to the area of forest land and was not affected by the slope, the incident angle of the view direction was equal to 0 degrees. Without considering the influence of the slope, the sensitivity of the slope was 0. The sensitivity of the radius of the crown, stem height, crown height, half apex angle, clumping index for shoots in the crown, shape of crown, LAI, tree number, leaf reflectivity, leaf transmittance, ground reflectivity, solar azimuth angle, and view azimuth angle were calculated, and the results are shown in Table 2.
The radius of the crown and the number of plants had a greater influence on the P vg-c (the gap fraction between the crowns), and the other parameters had no influence. This was because the canopy was considered to be an opaque canopy. Additionally, the canopy closure was only affected by the radius of the crown and individual trees, and it was not related to the leaf area index and clumping index. Therefore, when the radius of the crown and the number of trees were determined, the vertical projection area of the canopy on the forest could be estimated. However, when the gaps in the crown were considered, the LAI had the greatest effect on the P vg (total gap fraction of the canopy). The clumping index had the next largest effect, followed by the radius of the crown and the number of trees.

Estimation of the Canopy Gap Fraction Based on the GOST Model
The results of the parameter sensitivity analysis show that the number of trees (n), radius of the crown (d), and LAI (l) had a greater impact on P vg-c and P vg simulation. While n, d, and l were set as the measured values, other insensitive parameters were set as the average value of all of the sample plots. Then, the P vg-c and P vg of all of the plots were simulated by the GOST model. The max of the simulated P vg-c was 0.6671, the min value of the P vg-c was 0.2051, and the mean value was 0.4589. Additionally, the max, min, and mean values of the P vg were 0.7439, 0.3011, and 0.5641, respectively. Additionally, the value of P vg-c was obviously smaller than the simulated P vg value.
When the P vg-c and P vg of all of the plots were simulated, the relationship between the P vg-c and P vg and the field-measured data was analyzed. In our statistical analysis, we found that the linear relationships among the P vg-c and P vg and the single parameters of the number of trees and the radius of the crown were quite insignificant, respectively. In contrast, this combination of n × d 2 variables had a better linear relationship with the canopy gap fraction. Additionally, the n × d 2 variable was a variable because it represented the projected area of the forest canopy and was physically significant for P vg-c and P vg estimation compared to the single variables d or n.
As such, the relationships among the P vg-c , P vg , and n × d 2 were established ( Table 3). The established models passed the significance test. Table 3. Statistical models of P vg and P vg-c estimation using n × d 2 .

Dependent Variable
The results show that the linear relationship between n × d 2 and the P vg-c was significant, with the coefficient of determination being 0.5619 ( Figure 4). Additionally, the reason for this could be that the P vg-c only considered the gaps between the canopies and ignored the gaps in the crown. The variable n × d 2 represented the projected area of all of the canopies in the plot, taking the canopy as an opaque body without considering the gaps in the crown. Therefore, n × d 2 , which represents the projected area of the forest canopy, can describe the canopy closure well. Therefore, the linear relationship between n × d 2 and the P vg-c was more significant, and the coefficient of determination was larger than the one in Table 3.  The relationship between n × d 2 and the Pvg can be seen in Figure 5. The coef determination and the RMSE of the canopy closure estimation model were 0. 0.0813, respectively. Compared to the model established with the Pvg-c, these resu seem to be as good as the results above. This is because the gap between the can represented by the Pvg in the GOST model. Additionally, n × d 2 represented the pr of the canopy without considering the gaps in the canopy. Therefore, the linear ship between these two variables was weakened. At the same time, the results t estimated using this model underestimate the canopy due to the gaps in the can resented by the Pvg.  The relationship between n × d 2 and the P vg can be seen in Figure 5. The coefficient of determination and the RMSE of the canopy closure estimation model were 0.3138 and 0.0813, respectively. Compared to the model established with the P vg-c , these results do not seem to be as good as the results above. This is because the gap between the canopy was represented by the P vg in the GOST model. Additionally, n × d 2 represented the projections of the canopy without considering the gaps in the canopy. Therefore, the linear relationship between these two variables was weakened. At the same time, the results that were estimated using this model underestimate the canopy due to the gaps in the canopy represented by the P vg . The relationship between n × d 2 and the Pvg can be seen in Figure 5. The c determination and the RMSE of the canopy closure estimation model were 0.0813, respectively. Compared to the model established with the Pvg-c, these re seem to be as good as the results above. This is because the gap between the represented by the Pvg in the GOST model. Additionally, n × d 2 represented the of the canopy without considering the gaps in the canopy. Therefore, the line ship between these two variables was weakened. At the same time, the resul estimated using this model underestimate the canopy due to the gaps in the c resented by the Pvg.

Estimation of the Canopy Gap Fraction Based on the LAI
The canopy gap fraction can be estimated based on the number of trees and the radius of the crown. However, it was difficult to estimate the number of trees and the radius of the crown from most of the middle-resolution remote sensing data, and it turned out to be impossible. To solve this shortage, the LAI was used to solve this problem. Based on the equation to calculate the canopy gap fraction (Equation (5)), the LAI had a significant relationship with the canopy gap fraction. Additionally, the LAI was a key parameter for solving this problem. The LAI is a commonly used forest structural parameter, and there are many LAI estimation methods that are highly accurate [69,70]. This means that the LAI can be obtained from remote sensing images easily and that it can be used to estimate the canopy closure for middle-resolution remote sensing images. The relationships among the P vg-c , P vg , and LAI (l) were established and are shown in Table 4. These models passed the significance test. The results show that the linear relationships among the P vg-c , P vg , and LAI were significant, with the coefficients of determination being 0.2597 and 0.5467, respectively (Table 4). An exponential function between the P vg-c and LAI was observed, with the coefficient of determination being 0.2597, shown in Figure 6. Although this model passed the significance test, the accuracy of the simulation was not as good as expected. Additionally, these fitting results are also not as good as the results that were obtained using the variable n × d 2 .
Remote Sens. 2022, 14, x FOR PEER REVIEW solving this problem. The LAI is a commonly used forest structural paramet are many LAI estimation methods that are highly accurate [69,70]. This m LAI can be obtained from remote sensing images easily and that it can be use the canopy closure for middle-resolution remote sensing images. The relation the Pvg-c, Pvg, and LAI (l) were established and are shown in Table 4. These m the significance test. Pvg-c l y = 0.5762e −0.071x 0.2597 The results show that the linear relationships among the Pvg-c, Pvg, and L nificant, with the coefficients of determination being 0.2597 and 0.5467, respec 4). An exponential function between the Pvg-c and LAI was observed, with th of determination being 0.2597, shown in Figure 6. Although this model pass icance test, the accuracy of the simulation was not as good as expected. Addit fitting results are also not as good as the results that were obtained using the d 2 . The linear relationship between the Pvg and LAI was significant, with th of determination being 0.5467 and RMSE of 0.0654, respectively (Table 4). An function between the Pvg and LAI is also visible in Figure 7. These results a with Equation (5), but the relationship between the gap fraction and the LAI by other parameters of Equation (5). The reason for this is that the real expo tion was not as significant as Equation (5) The linear relationship between the P vg and LAI was significant, with the coefficients of determination being 0.5467 and RMSE of 0.0654, respectively (Table 4). An exponential function between the P vg and LAI is also visible in Figure 7. These results are consistent with Equation (5), but the relationship between the gap fraction and the LAI was affected Remote Sens. 2022, 14,1983 13 of 22 by other parameters of Equation (5). The reason for this is that the real exponential function was not as significant as Equation (5). In addition, compared to the results of the P vg-c , the P vg had a better linear relationship with the LAI. This could be related to the LAI, which was defined as half of the total leaf area of all of the canopies in the plot per unit of ground area [71]. As such, the gaps in and between the canopies were considered. Therefore, the linear relationship between the P vg and the LAI was better than that between the P vg-c and the LAI. At the same time, compared to the variable n × d 2 , the LAI was more suitable for P vg estimation, achieving R 2 and RMSE values of 0.5467 and 0.0654, respectively.

Verification of Estimation of Canopy Closure Based on Pvg-c
Based on the definition of canopy closure, only the gaps between considered, so the value of 1 − Pvg-c should be equal to the canopy closu sion model was established between the canopy closure and 1 − Pvg-c. values were 0.5216 and 0.0832, respectively (Figure 8), and the preci closure estimation was 86.69%. At the same time, the established mode icance test. It can be seen that the canopy closure measured by the l significantly good relationship with 1 − Pvg-c. Among the traditional m ing canopy closure, the line transect method was more accurate. The effective and feasible to estimate canopy closure using the GOST mod

Verification of Estimation of Canopy Closure Based on P vg-c
Based on the definition of canopy closure, only the gaps between the canopies were considered, so the value of 1 − P vg-c should be equal to the canopy closure. A linear regression model was established between the canopy closure and 1 − P vg-c . The R 2 and RMSE values were 0.5216 and 0.0832, respectively (Figure 8), and the precision of the canopy closure estimation was 86.69%. At the same time, the established model passed the significance test. It can be seen that the canopy closure measured by the line transects had a significantly good relationship with 1 − P vg-c . Among the traditional methods for measuring canopy closure, the line transect method was more accurate. Therefore, it was very effective and feasible to estimate canopy closure using the GOST model.

Verification of the Estimation of the Canopy Closure Based on P vg
The canopy closure calculated based on the fisheye camera photos included the gaps in the canopy. Therefore, calculating the canopy closure using the fisheye camera photos can verify the canopy closure estimated based on the P vg . The R 2 and RMSE values of the linear regression model were 0.1418 and 0.3295, respectively (Figure 9). The average accuracy of the canopy closure estimation was 73.2%. The model also passed the significance test. The results show that the relationship does not appear to be good. This could be because the fisheye camera images were taken at a 180 degree angle, resulting in the deformation increasing as the view angle became larger, while the GOST model calculated the P vg when the view angle was 90 degrees [72]. The reason for this was that the projected area of the crown between these two styles was different, so the linear relationship was not as significant.
sion model was established between the canopy closure and 1 − Pvg-c. Th values were 0.5216 and 0.0832, respectively (Figure 8), and the precisio closure estimation was 86.69%. At the same time, the established model p icance test. It can be seen that the canopy closure measured by the line significantly good relationship with 1 − Pvg-c. Among the traditional meth ing canopy closure, the line transect method was more accurate. Theref effective and feasible to estimate canopy closure using the GOST model.

Verification of the Estimation of the Canopy Closure Based on Pvg
The canopy closure calculated based on the fisheye camera photos in in the canopy. Therefore, calculating the canopy closure using the fisheye can verify the canopy closure estimated based on the Pvg. The Sens. 2022, 14, x FOR PEER REVIEW linear regression model were 0.1418 and 0.3295, respectively ( Figure 9). curacy of the canopy closure estimation was 73.2%. The model also passed test. The results show that the relationship does not appear to be good. T cause the fisheye camera images were taken at a 180 degree angle, result mation increasing as the view angle became larger, while the GOST mod Pvg when the view angle was 90 degrees [72]. The reason for this was th area of the crown between these two styles was different, so the linear not as significant.

Discussion
Canopy closure is a multipurpose ecological indicator that is used to assess the light conditions and forest floor microclimate as well as to distinguish types of plant and animal habitats [13]. Canopy closure also determines light interception and is crucial for understanding forest carbon fixation and responses and feedback to climate change [73]. Unfortunately, the in situ measurement of canopy closure is a time-consuming and laborious process, whilst advances in non-destructive, indirect techniques have been made [15]. Indirect optical methods derive canopy closure by measuring the canopy gap fraction or via transmittance [74]. These in situ measurements are able to capture the canopy at the plot (and site) scale, but these methods cannot adequately characterize spatiotemporal dynamics [32].
More and more remote sensing applications involve the estimation of canopy closure. The least squares method is one of the most commonly used methods for canopy closure estimation. However, the linear relationship is data dependent. This means that accuracy will be affected by the data quality and that applications will be limited by the area of the research region [75]. Other studies have used the random forest and Cubist models to estimate the canopy closure. However, the precision of these models is not accurate enough because a single variable is used for regression analysis [10]. In contrast, non-parametric models have always had better estimation results, and the accuracy of the modeling and estimation would be limited by the number of the samples, as small samples weaken the predictive ability of the model [76,77]. To make up for these shortages, we established a canopy closure estimation method that uses the gap fraction based on the geometric-optical GOST model, and the results were discussed.
A canopy closure estimation method that uses the gap fraction was represented in this study. The canopy closure P vg-c was highly correlated with the average projected area of the canopy. The P vg-c was estimated based on the number of plants in the plot and the average radius of the crowns, and the accuracy of the estimated canopy closure was 86.69%, while the accuracy of the canopy closure estimated by the P vg was 7%. In the residual analysis of the canopy closure model, the P vg-c , P vg , and n × d 2 (Figures 10 and 11) show that the residuals of the two models were distributed within a reasonable range and relatively evenly. The results show that the accuracy of the P vg-c estimation of the model was higher than that of the residual plots for P vg-c and n × d 2 . The residuals were distributed between ±0.2. However, the P vg was not estimated as well as the P vg-c was ( Figure 11). The residuals of the P vg were distributed between −0.2 and 0.3. It can be seen that estimating the canopy closure P vg-c between the canopies using the number of plants and the average radius of the crowns produced better results.
For most of the remote sensing images with a high spatial resolution, the number of plants and the radius of the crown were able to be extracted easily by means of image segmentation technology, and the parameters of these two variables can be considered an approach to canopy closure estimation using remote sensing images with a high spatial resolution. The relationship between the P vg-c and n × d 2 was the best optimal method, obtaining high-precision canopy closure estimation results. The line transect method was more accurate compared to other methods. As such, the scattering plot for the estimated canopy closure determined and measured by the gap fraction can be found in Figure 12. The canopy closure that was estimated by the P vg-c was consistent with the measured canopy closure. However, the canopy closure that was estimated by the P vg deviated from the line y = x. Additionally, the canopy closure was underestimated. between ±0.2. However, the Pvg was not estimated as well as the Pvg-c was (F residuals of the Pvg were distributed between −0.2 and 0.3. It can be seen th the canopy closure Pvg-c between the canopies using the number of plants an radius of the crowns produced better results.  For most of the remote sensing images with a high spatial resolution, t plants and the radius of the crown were able to be extracted easily by me segmentation technology, and the parameters of these two variables can be approach to canopy closure estimation using remote sensing images with resolution. The relationship between the Pvg-c and n × d 2 was the best optima taining high-precision canopy closure estimation results. The line transec between ±0.2. However, the Pvg was not estimated as well as the Pvg-c wa residuals of the Pvg were distributed between −0.2 and 0.3. It can be see the canopy closure Pvg-c between the canopies using the number of plant radius of the crowns produced better results.  For most of the remote sensing images with a high spatial resolutio plants and the radius of the crown were able to be extracted easily by segmentation technology, and the parameters of these two variables can approach to canopy closure estimation using remote sensing images w resolution. The relationship between the Pvg-c and n × d 2 was the best opt taining high-precision canopy closure estimation results. The line tran more accurate compared to other methods. As such, the scattering plot Compared to traditional statistical models, the model in the present study was not affected by the study area or the dataset, and there was no need to establish models in different study areas. The canopy closure estimation method simulated the probability of being able to see the ground by considering the mechanism of the light transmission process in the canopy and the gaps or opaque canopy in the viewing direction, which can also be considered individually using the geometric-optical model (the GOST model), and the estimated results also have high accuracy. In other words, this means that this estimated method was more robust compared to statistical methods. If remote sensing data with medium or low spatial resolution are used, then the parameters used for canopy estimation cannot be inverted accurately because of the mixed pixel effect [44]. Additionally, most of the parameters, such as the number of the trees and the radius of the crown could be obtained, so this method was difficult to implement. Instead, some medium-level variables, such as the LAI, should be used for canopy closure estimation even though the accuracy varies from study to study [78][79][80][81][82]. In this study, we evaluated the canopy closure estimation efficiency when using the LAI. The values of the residuals were more discrete from the residual P vg-c and LAI plots ( Figure 13). Additionally, the figure shows that the accuracy of the P vg-c estimation was better when the canopy closure was greater than 0.4 and when the value of residuals was lower than 0.1. When the canopy closure was greater than 0.4, the accuracy of the P vg-c prediction was obviously lower.
te Sens. 2022, 14, x FOR PEER REVIEW 16 of canopy closure determined and measured by the gap fraction can be found in Figure 1 The canopy closure that was estimated by the Pvg-c was consistent with the measured ca opy closure. However, the canopy closure that was estimated by the Pvg deviated from t line y = x. Additionally, the canopy closure was underestimated. Compared to traditional statistical models, the model in the present study was n affected by the study area or the dataset, and there was no need to establish models different study areas. The canopy closure estimation method simulated the probability being able to see the ground by considering the mechanism of the light transmission pr cess in the canopy and the gaps or opaque canopy in the viewing direction, which can al be considered individually using the geometric-optical model (the GOST model), and t estimated results also have high accuracy. In other words, this means that this estimat method was more robust compared to statistical methods. If remote sensing data wi medium or low spatial resolution are used, then the parameters used for canopy estim tion cannot be inverted accurately because of the mixed pixel effect [44]. Additional most of the parameters, such as the number of the trees and the radius of the crown cou be obtained, so this method was difficult to implement. Instead, some medium-level va iables, such as the LAI, should be used for canopy closure estimation even though t accuracy varies from study to study [78][79][80][81][82]. In this study, we evaluated the canopy cl sure estimation efficiency when using the LAI. The values of the residuals were more d crete from the residual Pvg-c and LAI plots ( Figure 13). Additionally, the figure shows th the accuracy of the Pvg-c estimation was better when the canopy closure was greater th 0.4 and when the value of residuals was lower than 0.1. When the canopy closure w greater than 0.4, the accuracy of the Pvg-c prediction was obviously lower.
The situation with the Pvg was similar to that of the Pvg-c. The residuals of the Pvg b came significantly discrete from the residual Pvg and LAI plots when the canopy closu was larger than or less than 0.45 ( Figure 14). The most possible reason for this is that mo gaps in the crown could be detected when using remote sensors with a low or mediu Sens. 2022, 14, x FOR PEER REVIEW the pixel binary model, and the robustness and accuracy of the model w In addition, the results also indicate that the LAI can be used to estimate t accuracy. Once high-quality LAI data were obtained, the canopy closu estimated, especially for remote sensing images with a medium or low s This method provided a way to estimate the canopy closure using the LA this relationship was derived from geometric-optical models and was not which is the case using statistical theory. Additionally, this method can the shortcomings of canopy estimation models that are based on statistic The situation with the P vg was similar to that of the P vg-c . The residuals of the P vg became significantly discrete from the residual P vg and LAI plots when the canopy closure was larger than or less than 0.45 ( Figure 14). The most possible reason for this is that more gaps in the crown could be detected when using remote sensors with a low or medium spatial resolution, and the efficiency of the mixed pixels enhanced the signals of the forest floor or those of other features. Inversely, the larger amount of canopy closure meant smaller gaps in the crown, and the detected probability of the signal outside of the canopy was lower compared to larger gaps, weakening the mixed pixel effect. As such, it was necessary to consider the efficiency of the mixed pixels in the images when the LAI was used to estimate the canopy closure. In comparison, a geometric-optical model was used to decompose the mixed pixels into four components, allowing the forest components in the pixels to be more accurately distinguished. That was the advantage of the method using a physical model. At the same time, the physical canopy closure estimation procedure using a geometric-optical model was clearer than the traditional statistical model or the pixel binary model, and the robustness and accuracy of the model were self-evident. In addition, the results also indicate that the LAI can be used to estimate the P vg with good accuracy. Once high-quality LAI data were obtained, the canopy closure could also be estimated, especially for remote sensing images with a medium or low spatial resolution. This method provided a way to estimate the canopy closure using the LAI. Additionally, this relationship was derived from geometric-optical models and was not data dependent, which is the case using statistical theory. Additionally, this method can effectively avoid the shortcomings of canopy estimation models that are based on statistical methods.

Conclusions
Canopy closure is an important forest inventory parameter, an role in forestry production and management and forest health ev the GOST model was used to simulate the characteristics of fores by considering an opaque canopy (Pvg-c) and gaps in the tree cro exponential models for estimating the canopy closure of a plant

Conclusions
Canopy closure is an important forest inventory parameter, and it plays an important role in forestry production and management and forest health evaluation. In this study, the GOST model was used to simulate the characteristics of forest canopy gap fractions by considering an opaque canopy (P vg-c ) and gaps in the tree crown (P vg ). Additionally, exponential models for estimating the canopy closure of a plantation based on the gap fraction were established. The results show the following: (1) It was feasible to estimate canopy closure based on the GOST model, and the feasible method was proved with sample data measured from three different regions in China. (2) Compared to the LAI, n × d 2 had a better relationship with the gap fractions simulated using the GOST model. Therefore, when remote sensing images or LiDAR data of the study area with high spatial resolution were available, the crown recognition method could be used to obtain the number of plants and the average radius of the crowns in the plot, so the gap fraction P vg-c and the forest canopy closure could be accurately estimated and predicted in the research area. (3) When the number of plants and the average radii of the crowns in the plot could not be extracted using remote sensing images, especially when only medium-or low-spatial resolution remote sensing data were available, the LAI, a medium parameter, could be used to estimate the canopy closure with an acceptable level of accuracy. This also provided a new a canopy closure estimation approach using medium-or low-spatial resolution remote sensing data. This study can provide a reference for canopy closure estimation using geometric-optical models.

Conflicts of Interest:
The authors declare no conflict of interest.