Effects of Canopy and Multi-Epoch Observations on Single-Point Positioning Errors of a GNSS in Coniferous and Broadleaved Forests

: Global navigation satellite systems (GNSS) can quickly, efﬁciently, and accurately provide precise coordinates of points, lines, and surface elements, plus complete surveys and determine various boundary lines in forest investigations and management. The system has become a powerful tool for dynamic forest resource investigations and monitoring. GNSS technology plays a unique and important role in estimating timber volume, calculating timber cutting area, and determining the location of virgin forest roads and individual trees in forests. In this study, we quantitatively analyzed the inﬂuence of crown size and observation time on the single-point positioning accuracy of GNSS receivers for different forest types. The GNSS located single points for different forest types and crown sizes, enabling the collection of data. The locating time for each tree was more than 10 min. Statistical methods were used to analyze the positioning accuracy of multi-epoch data, and a model was developed to estimate the maximum positioning errors under different forest conditions in a certain positioning time. The results showed that for a continuous positioning time of approximately 10 min, the maximum positioning accuracies in coniferous and broadleaf forests were obtained, which were 12.13 and 15.11 m, respectively. The size of a single canopy had no obvious inﬂuence on the single-point positioning error of the GNSS, and canopy density was proven to be closely related to the positioning accuracy of a GNSS. The determination coefﬁcients (R 2 ) in the regression analysis of the general model, coniferous forest model, and broadleaved forest model that were developed in this study were 0.579, 0.701, and 0.544, respectively. These results indicated that the model could effectively predict the maximum positioning error in a certain period of time under different forest types and crown conditions at middle altitudes, which has important guiding signiﬁcance for forest resource inventories and precise forest management. the satellite distribution in the sky at the time of positioning and the occlusion of the tree crown or other obstacles. The residuals of the carrier phase and pseudo-range did not decrease signiﬁcantly with increasing positioning time, which indicates that the satellite signals received by the receiver in the positioning process did not increase with increasing positioning time. This information is consistent with the description in the results. The positioning accuracy of the GNSS for the coniferous forest was higher than that for the broadleaved forest, which is consistent with the results of Deckert and Bolstad [4], but also differed from those of Deckert and Bolstad [4], i.e., GNSS positioning error was not affected by crown type. In this study, the minimum positioning error, 12.13 m in the coniferous forest, is signiﬁcantly higher than that in a coniferous forest determined by Næsset [1] of 1.98 m. Many factors contributed to this large difference; among them, geographical location may be one of the most important.


Introduction
Many forest resource surveys require highly accurate field measurements of the spatial locations of trees. According to the different requirements of practical applications, the accuracy of locating a single tree or sample site varies from several centimeters to several meters [1]. With increasing attention focused on the spatial locations of trees and sample plots in forestry surveys, exploring which factors are the dominant factors affecting spatial positioning and determining how to quantitatively express the impact of those factors on positioning accuracy are becoming necessary [2,3].
The development of satellite positioning is rapidly providing new options for forest resource managers to obtain timely and highly accurate location data [4]. Global navigation satellite systems (GNSSs) mainly include GPS in the United States, Russia's GLONASS global satellite navigation system, the European Union's Galileo system, and China's Beidou satellite navigation system, BDS. Among the four navigation and positioning systems, GPS and GLONASS are available globally, and the other two systems remain under construction. At the end of June 2020, China successfully launched the 55th navigation satellite of the Beidou system and the last global networking satellite of Beidou-3. So far, the constellation deployment of the Beidou global navigation satellite system has been completed. Recently, airborne laser scanning to create forest resource inventories has become a practical method for mapping large areas of forest resources in detail [5,6]. The practical application of airborne laser scanning relies on the widespread use of millimeter-scale receivers and the differential positioning of field sample maps. The acquisition cost of these receivers is high, and managing differential positioning data is time-consuming in forest resource inventories, especially in remote areas of reference stations that require time to collect the data and are located away from the forest [7]. In many countries, continuously operating reference stations are used for differential positioning. However, in remote areas, the distance from a reference station reduces the accuracy of this positioning. If a precise position can be obtained by a technique that does not rely on differential correction, the data processing of the receiver can be simplified, and the cost can be reduced. GNSSs have increasingly been playing roles in positioning forest fires [8], mapping forest resources [9], inventorying forest resources [10], positioning fixed sample sites [11], mapping roads [12], and measuring other natural resources [13].
Positioning with a GNSS in a forest can be blocked by trunks and branches, or the satellite signal can be reflected to produce a multipath effect; therefore, the receiver cannot be located in the forest, or the positioning accuracy is low [14,15]. The canopy is an important variable affecting the statistical characteristics of GNSSs, such as path loss, attenuation, and polarization [16]. The canopy density is a key factor affecting GNSS positioning accuracy, and many studies have been conducted on the influence of density indexes on positioning accuracy [11,[17][18][19]. In general, the smaller the canopy density, the wider the sky, the higher the positioning accuracy and vice versa. The forest type is also an important factor affecting the positioning of GNSSs. Although many previous experimental studies were conducted, these studies were not able to effectively quantify the extent to which forest types affect the positioning of GNSSs [20]. The positioning time of GNSSs under the tree crowns is also an important factor affecting the positioning accuracy of GNSSs. Reasonable positioning time and appropriate positioning accuracy are important considerations for forest resource inventories and work efficiency [21]. The season is another important factor. It is generally thought that the presence of foliage has a strong impact on GNSS signals, reducing the accuracy of positioning, whereas positioning accuracy is higher during seasons without leaves [22][23][24].
To explore the quantitative effects of these factors on the positioning errors of GNSSs, a GNSS equipment system was used to conduct experiments in different forest types and to quantitatively analyze and express the results. The main objectives of this study were: (1) to explore the impact of coniferous and broadleaved tree species on GNSS positioning errors; (2) to quantitatively analyze the potential relationship between positioning time and positioning accuracy, and to determine the optimal positioning time for forest inventories;

Study Area
For the study area, we chose the Changpoling National Forest Park in Baiyun District, Guiyang city, Guizhou province, China. The selected study area is a typical artificial forest in the southwestern mountainous area of China, with a total area of 1075 hectares; this area is an important part of the green belt around the city of Guiyang (Figure 1). The altitude of the area ranges from 2000 to 2500 m; the long-term annual minimum temperature and maximum daily air temperature averages are 12 and 19 • C, respectively; and the forest coverage rate is as high as 82.96%. The forest types include coniferous forest, broadleaved forest, and mixed coniferous and broadleaved forest. The main species are Populus L., Pinus, Pinus armandii Franch, Pinus massoniana Lamb, Ginkgo biloba L., Liriodendron chinensis, and Eucommia ulmoides. Additional shrubs and herbs are found in the forest. In this study, because the selected sampling data area was small, the topographical terrain, slope direction, elevation, latitude, and longitude did not vary considerably. Therefore, the impact of these factors on the single-point positioning accuracy of the Galaxy 6 GNSS among different trees was not considered in this study. However, due to the relatively large forest density, the multipath effect will have a certain impact on the positioning accuracy. In order to reduce the influence of the multipath effect on the positioning accuracy, we extended the observation time appropriately.

Study Area
For the study area, we chose the Changpoling National Forest Park in Baiyun Distr Guiyang city, Guizhou province, China. The selected study area is a typical artificial for in the southwestern mountainous area of China, with a total area of 1075 hectares; t area is an important part of the green belt around the city of Guiyang (Figure 1). The a tude of the area ranges from 2000 to 2500 m; the long-term annual minimum temperat and maximum daily air temperature averages are 12 and 19 °C, respectively; and the for coverage rate is as high as 82.96%. The forest types include coniferous forest, broadleav forest, and mixed coniferous and broadleaved forest. The main species are Populus L., nus, Pinus armandii Franch, Pinus massoniana Lamb, Ginkgo biloba L., Liriodendron chinen and Eucommia ulmoides. Additional shrubs and herbs are found in the forest. In this stu because the selected sampling data area was small, the topographical terrain, slope dir tion, elevation, latitude, and longitude did not vary considerably. Therefore, the imp of these factors on the single-point positioning accuracy of the Galaxy 6 GNSS amo different trees was not considered in this study. However, due to the relatively large for density, the multipath effect will have a certain impact on the positioning accuracy. order to reduce the influence of the multipath effect on the positioning accuracy, we tended the observation time appropriately.

Data Acquisition
We used a Galaxy 6 RTK ( Figure 2a) (South Surveying & Mapping Technology C Ltd.; http://www.southinstrument.com/ (15/3/2021), Guangzhou city, China) to acqu the location data of the coniferous and broadleaf forests in the field. The positioning o put frequency of the Galaxy 6 RTK GNSS is 1-50 Hz, and the initial reliability is grea than 99.99%, which greatly improves the accuracy of carrier positioning by relying highly reliable carrier tracking technology. The lithium battery capacity of the devic 6800 mAh, and the continuous working time is greater than 30 h in the static stand mode. The static measuring accuracy of this device is theoretically 2.5 mm. With ever km increase in measuring distance between two points, the positioning accuracy creases by 1 mm. The horizontal positioning accuracy of the single-point code differen GNSS is 0.25 m + 1 ppm RMS, which can provide relatively more stable data by rely on mature CORS technology. The technical specifications of Galaxy 6 RTK instrument shown in Table 1. In this experiment, we set the output frequency of the Galaxy 6 R GNSS to 1 Hz; that is, the sampling interval was 1 s. We conducted static single-po positioning for each tree in the experiment for at least 10 min to obtain tree location inf mation.

Data Acquisition
We used a Galaxy 6 RTK ( Figure 2a) (South Surveying & Mapping Technology Co., Ltd.; http://www.southinstrument.com/ (15 March 2021), Guangzhou city, China) to acquire the location data of the coniferous and broadleaf forests in the field. The positioning output frequency of the Galaxy 6 RTK GNSS is 1-50 Hz, and the initial reliability is greater than 99.99%, which greatly improves the accuracy of carrier positioning by relying on highly reliable carrier tracking technology. The lithium battery capacity of the device is 6800 mAh, and the continuous working time is greater than 30 h in the static standard mode. The static measuring accuracy of this device is theoretically 2.5 mm. With every 1 km increase in measuring distance between two points, the positioning accuracy decreases by 1 mm. The horizontal positioning accuracy of the single-point code differential GNSS is 0.25 m + 1 ppm RMS, which can provide relatively more stable data by relying on mature CORS technology. The technical specifications of Galaxy 6 RTK instrument are shown in Table 1. In this experiment, we set the output frequency of the Galaxy 6 RTK GNSS to 1 Hz; that is, the sampling interval was 1 s. We conducted static single-point positioning for each tree in the experiment for at least 10 min to obtain tree location information.     in Figure 2b) because we could not directly locate a tree's core coordinates when locating a tree's position in practice; therefore, the coordinate values were measured a certain distance from the tree's center (30 cm to the south). In addition, we used the total station to measure the coordinates of five open area points (no trees or other obstacles nearby) based on the coordinates of the known points, A and B, and labeled them C ( The five points without shelter (C, D, E, F, and G) were positioned using Galaxy 6 RTK GNSS equipment to obtain the GNSS single point positioning error, σ 0 , without a canopy or other shelter. The formulas are defined as follows: where ∆x i and ∆y i are the differences between the measured values obtained by the Galaxy 6 RTK GNSS in the x and y directions and the reference values obtained by the total station, respectively; x 0 and y 0 are the positioning system errors of the Galaxy 6 RTK system in the x and y directions, respectively, and k is the number of the 5 points in the open area. The coordinate measurements of 29 coniferous and broadleaved trees in the plot were obtained using the Galaxy 6 RTK system 30 cm in the southward direction of each tree ( Figure 2b). To obtain accurate tree position coordinate information and account for the influence of positioning time on work efficiency, the single-point positioning time of each tree position was at least 10 min or longer to ensure that sufficient epoch data were obtained.
Field data were collected in April 2019 ( Figure 3). The tree structure parameters obtained in the field included diameter at breast height (DBH), tree height (h), canopy (D), and tree species (coniferous or broadleaved). The DBH of each tree was measured using a measuring tape with millimeter accuracy. The heights of the trees were measured using a mobile RGBD-SLAM tree measurement system developed by the Beijing Key Laboratory of Precision Forestry [6]; the tree height accuracy measured by this system was between 10 and 50 cm and was obtained using two optical cameras and a field camera with infrared depth. A tape measure was used to measure the size of the canopy, and the species were identified by a forestry professional. The structural characteristics of coniferous and broadleaved forests based on field forest inventory data are shown in Table 2.

Data Processing
In this study, we only considered the positioning accuracy of the RTK GNSS plane coordinates and did not analyze the elevation positioning accuracy. The data obtained by the Galaxy 6 RTK GNSS for single-point positioning were in .sth format. Information re-

Data Processing
In this study, we only considered the positioning accuracy of the RTK GNSS plane coordinates and did not analyze the elevation positioning accuracy. The data obtained by the This software, which has relatively high accuracy, is the most commonly used tool for calculating epoch data. Therefore, we input the converted data into the RTKPOST application program, which is a part of RTKLIB. Because the Galaxy 6 RTK GNSS only receives carrier information from the GPS and GLONASS satellite systems when acquiring single-point positions, the singlepoint positioning information of trees in the sample plot only included positioning data from the GPS and GLONASS satellite systems. By analyzing the configuration file of the remote sensing satellite data set, the spatial continuity test parameters of the corresponding remote sensing satellite data set were obtained; Based on the remote sensing satellite data set and the spatial continuity test parameters, we observed that the remote sensing satellite data set did not have serious missing data, which indicated that the satellite signal had good continuity [27]. The coordinates of the epoch data calculated from singlepoint positioning data are the geodesic latitude and longitude coordinates in the WGS-84 coordinate system, i.e., two sets of data with completely different coordinate systems, the WGS-84 and CGCS2000 coordinate systems, which were obtained by the total station. COORD (an open-source program package for transforming coordinates, version 4.2) is a commonly used coordinate conversion software program that converts coordinate information among the Beijing-54 coordinate system, the Xi'an-80 coordinate system, the WGS-84 coordinate system, and the WGS-72 coordinate system using seven or four parameters. We used this software to convert the geodetic coordinates and plane coordinate points in the WGS-84 coordinate system. Before the coordinate transformation, we first determined the central meridian of the study area. A 3-degree zoning projection ensured the necessary accuracy for a large-scale topographic map. The central meridian of the area was determined with a 3-degree banding tape number and the longitude information of the study area.
Although the basic definition of the WGS-84 coordinate system is consistent with that of the CGCS2000 coordinate system, their reference ellipsoids are similar, and only the oblateness of the ellipsoid constant is slightly different. However, to ensure strict consistency with the coordinates measured by the total station, the plane coordinates of WGS-84 were converted to the plane coordinates of CGCS2000. The calculated data of each epoch were compared with the corresponding data measured by the total station. The positioning error between the measured data and the reference value is defined as follows: where ∆x i and ∆y i are the differences between the measured values obtained by the Galaxy 6 RTK GNSS in the x and y directions and the reference values obtained by the total station, respectively. The epoch data of the single-point positioning of each tree were sorted according to the positioning error, ∆p i , from large to small. Then, in accordance with 1, 2, 5, 10...400, 500, and all epochs, the epoch data were grouped to obtain the average positioning error and the positioning errors in the x and y directions were analyzed separately. By analyzing and comparing the mean values of different multi-epoch data, the quantity of multi-epoch data with minimum positioning error was determined.

Statistical Analysis
To determine the maximum positioning error of different amounts of multi-epoch data at a single positioning time point under the tree canopy using the Galaxy 6 RTK GNSS system, we calculated bias and determined the root mean square error (RMSE) of the coordinates to compare the positioning accuracy of the RTK GNSS system with the reference values measured by the total station. BIAS and RMSE are defined as follows: where ∆x i and ∆y i are the differences between the measured values obtained by the Galaxy 6 RTK GNSS in the x and y directions and the reference values obtained by the total station, respectively; and m is the total number of trees measured. Root mean square coordinate error (RMSE xy ) is one of the most commonly used indicators for evaluating measurement accuracy in geodetic surveying, and this measurement is also a characteristic of point set accuracy [24]. RMSE xy is calculated as follows:

Construction of the Positioning Error Model
To further explore the relationships between the RTK GNSS positioning accuracy and tree crown size, tree species, and multi-epoch observation data, we built a general model of the RTK GNSS positioning error under different tree crowns: where σ i is the error of the ith tree positioning, n is the number of epochs, i indicates the ith tree, a 2 is a coefficient related to the number of epochs, b 1 and b 2 are coefficients associated with the crown diameters of coniferous and broadleaved forests, respectively; and D 1i and D 2i are the crown diameter values of coniferous and broadleaved forest species, respectively. In our constructed positioning error model, when the ith tree was a coniferous tree, the coefficient b 2 of D 2i was 0, and then the form of the error model was as follows: Similarly, when the ith tree was a broadleaved tree, the form of the error model was as follows: σ i = σ 2 0 + a 2 n + b 2 2 * D 2 2i . We analyzed the correlations between independent variables n and D and the dependent variable σ using IBM SPSS Statistics 23 software (International Business Machines Corporation, https://www.ibm. com/cn-zh, 12 June 2021), which offers advanced statistical analysis, a vast library of machine-learning algorithms, text analysis, open-source extensibility, integration with big data, and seamless deployment in applications. We also used the software to conduct linear regression analysis on our established general model to verify the reliability and accuracy of the model. In addition, coniferous and broadleaved forest species were modeled and analyzed separately to determine the impact of different forest types on positioning accuracy.

Results
For a single navigation and positioning system, a GNSS needs to track at least four satellites for three-dimensional positioning; for a dual system, the system needs to track at least five satellites. Therefore, there is a case when positioning epoch data are missing in actual GNSS positioning data processing due to satellite signal loss. The satellite distribution at a certain moment during the acquisition of single-point positioning data in the field is shown in Figure 4. We found that of 29 trees, only 1 had fewer than 500 epoch data points (489) when calculating the multi-epoch data, which did not affect our calculation and analysis of the epoch data.
error model was as follows: = + + * . We analyzed the correlations between independent variables n and D and the dependent variable σ using IBM SPSS Statistics 23 software (International Business Machines Corporation, https://www.ibm.com/cn-zh,2019), which offers advanced statistical analysis, a vast library of machine-learning algorithms, text analysis, open-source extensibility, integration with big data, and seamless deployment in applications. We also used the software to conduct linear regression analysis on our established general model to verify the reliability and accuracy of the model. In addition, coniferous and broadleaved forest species were modeled and analyzed separately to determine the impact of different forest types on positioning accuracy.

Results
For a single navigation and positioning system, a GNSS needs to track at least four satellites for three-dimensional positioning; for a dual system, the system needs to track at least five satellites. Therefore, there is a case when positioning epoch data are missing in actual GNSS positioning data processing due to satellite signal loss. The satellite distribution at a certain moment during the acquisition of single-point positioning data in the field is shown in Figure 4. We found that of 29 trees, only 1 had fewer than 500 epoch data points (489) when calculating the multi-epoch data, which did not affect our calculation and analysis of the epoch data

Positioning Accuracy of Different Multi-Epoch Data
The epoch data of the single-point positioning of 29 trees in the sample were grouped by different quantities, and the mean positioning error of each group of epoch data was calculated. The positioning error results of different groups of multi-epoch data are shown in Table 3. From the table, with increases in the amount of multi-epoch data, the positioning errors of bias and RMSE on the x and y axes generally showed a decreasing trend, and

Positioning Accuracy of Different Multi-Epoch Data
The epoch data of the single-point positioning of 29 trees in the sample were grouped by different quantities, and the mean positioning error of each group of epoch data was calculated. The positioning error results of different groups of multi-epoch data are shown in Table 3. From the table, with increases in the amount of multi-epoch data, the positioning errors of bias and RMSE on the x and y axes generally showed a decreasing trend, and the total positioning error RMSE xy also decreased gradually. The GNSS single-point positioning error was greatly reduced from the 1 epoch group to the 90 epoch groups.  All trees in the sample plot were also precisely located using the CSRS-PPP service provided by the Canadian Geodetic Survey of Natural Resources Canada, which has high computational accuracy. The error plane results of the precise single-point calculation of positioning data in the RTK GNSS are shown in Figure 5. As shown in the figure, the positions of all trees in the y-axis direction (i.e., longitude direction) were overestimated, and the positions in the x-axis direction (i.e., latitude direction) were both overestimated and underestimated, which were evenly distributed on both sides of the x = 0 line. the accuracy of single-point positioning improved, but the rate of improvement was lower than that of the previous multi-epoch data group. The positioning time was approximately 10 min; the longer the time allowed for determining a single-point location, the more accurate the location of the tree. All trees in the sample plot were also precisely located using the CSRS-PPP service provided by the Canadian Geodetic Survey of Natural Resources Canada, which has high computational accuracy. The error plane results of the precise single-point calculation of positioning data in the RTK GNSS are shown in Figure 5. As shown in the figure, the positions of all trees in the y-axis direction (i.e., longitude direction) were overestimated, and the positions in the x-axis direction (i.e., latitude direction) were both overestimated and underestimated, which were evenly distributed on both sides of the x = 0 line.  Table 4. The maximum bias was 31.41 m, the minimum bias was 10.61 m, and the median of the positioning error was 19.81 m. Positioning errors less than 15 m in the x-axis direction and y-axis direction accounted for 44.83% and 68.97% of the total, respectively. From the data in Table 4, with increasing tree crown cover, the error of single-point positioning did not increase or decrease significantly for either the broadleaved or coniferous forest species.

Comparison of Single-Point Positioning Errors between Coniferous and Broadleaved Forest Species
To quantitatively analyze the impact of forest species (broadleaved or coniferous) on the positioning of the GNSS, the positioning accuracy of different multi-epoch data was evaluated and analyzed in more detail ( Table 5). The RMSE xy of the broadleaved forest species ranged from 15.85 to 24.58 m, and that of the coniferous forest species ranged from 14.02 to 22.40 m. The position error range of the coniferous forest species was slightly narrower than that of the broadleaved forest species. For the same amount of multi-epoch data, the RMSE xy of the coniferous forest species was lower than that of the broadleaved forest species, and the positioning accuracy was higher. The positioning errors of the GNSS in the x-and y-axis directions under the broadleaved forest were greater than those under the coniferous forest, which were −9.26 and 11.94, and −4.95 and 11.07 m, respectively. The highest positioning accuracy was obtained in the total positioning time, with values of 15.11 and 12.13 m, respectively. Table 5. Accuracy evaluation of single-point positioning error in the coniferous forest and broadleaf forest.

Verification and Analysis of the Positioning Error Model
The construction and accuracy of the maximum positioning error model were evaluated for the different forest types. Table 6 provides the model regression analysis results of the general model, the coniferous forest model, and the broadleaved forest model. From the analysis results of the three models, the regression analysis coefficient (R 2 adj = 0.701) of the coniferous forest model was found to be significantly higher than those of the general model (R 2 adj = 0.579) and the broadleaved forest model (R 2 adj = 0.544). The independent variable factors of the three models were multi-epoch data and tree canopy. The model analysis results showed that the positioning error with the GNSS was very significant. The parameters of the different models are also shown in Table 6.

Discussion
In this study, we quantitatively analyzed and determined the influences of positioning time, forest type, and tree canopy on positioning errors using single-point positioning data in a GNSS under different forest types, and further built the maximum error model of different forest types within a certain positioning time. The results showed that with an increase in positioning time, the positioning accuracy of a single point under the forest improved considerably and gradually stabilized. The positioning accuracy under the tree canopy increased obviously as the amount of multi-epoch data increased, but the overall impact was an increase in positioning time. Single-point positioning in the different forest types achieved the highest positioning accuracy when the positioning time was approximately 10 min. Figure 6 shows the changes in GNSS positioning accuracy in the longitude and latitude directions with positioning time for the different forest types (coniferous and broadleaved forests). As shown in the figure, the positioning error ranges of the GNSS in the longitude and latitude directions under the coniferous forest were narrow and tended to more quickly achieve stable states. Compared with the results for the coniferous forest, the GNSS had larger positioning error deviations in the longitude and latitude directions under the broadleaved forest, and single-point positioning tended to achieve stability more slowly. However, the positioning errors in both the coniferous and broadleaved forests decreased as the positioning time increased, which also indicates that increases in positioning time help improve the accuracy of single-point positioning. Figure 7 depicts the distributions of the residuals of the carrier phase and pseudo-range received by different satellites for the coniferous and broadleaved forests with positioning time. As shown in the figure, the carrier phase and pseudo-range information received by the receiver from different satellites differed noticeably, which was closely related to the satellite distribution in the sky at the time of positioning and the occlusion of the tree crown or other obstacles. The residuals of the carrier phase and pseudo-range did not decrease significantly with increasing positioning time, which indicates that the satellite signals received by the receiver in the positioning process did not increase with increasing positioning time. This information is consistent with the description in the results. The positioning accuracy of the GNSS for the coniferous forest was higher than that for the broadleaved forest, which is consistent with the results of Deckert and Bolstad [4], but also differed from those of Deckert and Bolstad [4], i.e., GNSS positioning error was not affected by crown type. In this study, the minimum positioning error, 12.13 m in the coniferous forest, is significantly higher than that in a coniferous forest determined by Naesset [1] of 1.98 m. Many factors contributed to this large difference; among them, geographical location may be one of the most important.
positioning time help improve the accuracy of single-point positioning. Figure 7 depicts the distributions of the residuals of the carrier phase and pseudo-range received by different satellites for the coniferous and broadleaved forests with positioning time. As shown in the figure, the carrier phase and pseudo-range information received by the receiver from different satellites differed noticeably, which was closely related to the satellite distribution in the sky at the time of positioning and the occlusion of the tree crown or other obstacles. The residuals of the carrier phase and pseudo-range did not decrease significantly with increasing positioning time, which indicates that the satellite signals received by the receiver in the positioning process did not increase with increasing positioning time. This information is consistent with the description in the results. The positioning accuracy of the GNSS for the coniferous forest was higher than that for the broadleaved forest, which is consistent with the results of Deckert and Bolstad [4], but also differed from those of Deckert and Bolstad [4], i.e., GNSS positioning error was not affected by crown type. In this study, the minimum positioning error, 12.13 m in the coniferous forest, is significantly higher than that in a coniferous forest determined by Naesset [1] of 1.98 m. Many factors contributed to this large difference; among them, geographical location may be one of the most important.  The results of this study showed that the positioning error of a single point was relatively less affected by the size of a single crown. Previous studies have found that the positioning accuracy of a GNSS is strongly affected by stand density and canopy size [18,[28][29][30][31]. The GNSS performed single-point positioning under a coniferous forest and a broadleaved forest with different canopy sizes, but the positioning accuracy did not improve as the tree canopy decreased, and there was no strong correlation between the factors. However, the regression analysis of the constructed model showed that measurement time, tree canopy size, and positioning error were extremely significant, which agreed with our hypothesis. The determination coefficients of the regression analysis for the general model, coniferous model and broadleaf model were 0.579, 0.701, and 0.544, respectively, indicating that the coniferous forest model had the most accurate model prediction effect. One possible reason for this finding is that conifers have more space between their crowns, which helps pick up stable satellite signals to improve positioning accuracy. The model's prediction accuracy was slightly lower than those of different models constructed by Holden et al. [31] using differential GNSS data (R 2 = 0.81, R 2 = 0.74, and R 2 = 0.78). However, previous studies mainly focused on the impact of forest canopy on GNSS positioning and did not establish an effective GNSS positioning error prediction model [29,32,33]. The maximum error model established in this study can effectively predict the maximum positioning error of a particular positioning time to determine whether the positioning error meets the positioning accuracy requirements in advance.
In this study, a quantitative analysis of positioning errors of a GNSS system in coniferous and broadleaved forests was conducted. The relationship between positioning errors and positioning time length was described in detail. In future studies, we will select different geographic locations, terrain conditions, satellite systems, and types of receivers to test and verify the impact of these factors on positioning errors and consider how to quantify these effects. The results of this study showed that the positioning error of a single point was relatively less affected by the size of a single crown. Previous studies have found that the positioning accuracy of a GNSS is strongly affected by stand density and canopy size [18,[28][29][30][31]. The GNSS performed single-point positioning under a coniferous forest and a broadleaved forest with different canopy sizes, but the positioning accuracy did not improve as the tree canopy decreased, and there was no strong correlation between the factors. However, the regression analysis of the constructed model showed that measurement time, tree canopy size, and positioning error were extremely significant, which agreed with our hypothesis. The determination coefficients of the regression analysis for the general model, coniferous model and broadleaf model were 0.579, 0.701, and 0.544, respectively, indicating that the coniferous forest model had the most accurate model prediction effect. One possible reason for this finding is that conifers have more space between their crowns, which helps pick up stable satellite signals to improve positioning accuracy. The model's prediction accuracy was slightly lower than those of different models constructed by Holden et al. [31] using differential GNSS data (R 2 = 0.81, R 2 = 0.74, and R 2 = 0.78). However, previous studies mainly focused on the impact of forest canopy on GNSS positioning and did not establish an effective GNSS positioning error prediction model [29,32,33]. The maximum error model established in this study can effectively predict the maximum positioning error of a particular positioning time to determine whether the positioning error meets the positioning accuracy requirements in advance.
In this study, a quantitative analysis of positioning errors of a GNSS system in conif-

Conclusions
In this study, single-point positioning coordinates were obtained with the Galaxy 6 GNSS under different conditions. The results were compared with reference coordinates measured with a total station, and the influences of different forest types, positioning time, and crown size on the single-point positioning accuracy were determined. Finally, the maximum positioning error model of single-point positioning with a GNSS for a certain period of time was developed. The results showed that the different forest types had certain impacts on the positioning accuracy of the GNSS. The highest positioning accuracies of the GNSS under the coniferous and broadleaved forests were 15.11 and 12.13 m, respectively. This result showed that the positioning error of the GNSS was less affected by the coniferous forest, and the single-point positioning accuracy obtained was higher. For a certain positioning time, the positioning accuracy increased as the positioning time increased; the more the epoch data, the higher the positioning accuracy, which is important for obtaining the locations of trees in forests. Our study also showed that single crown size had little effect on the GNSS positioning error. In the future, we will obtain more experimental data to verify the results of this experiment and reduce unintentional errors. The maximum positioning error model developed in this study can effectively predict the positioning error for a forest in a forest resource inventory and determine in advance whether the positioning accuracy meets the accuracy requirements of the forest resource survey. In the future, we will conduct in-depth experiments under different terrain conditions and in different geographical locations to verify the impact of these factors on the positioning accuracy of the GNSS.