Construction and Accuracy Analysis of a BDS/GPS-Integrated Positioning Algorithm for Forests

The objective of this study was to construct a BeiDou navigation satellite system (BDS)/global positioning system (GPS)-integrated positioning algorithm that meets the accuracy requirement of forest surveys and to analyze its accuracy to provide theoretical and technical support for accurate positioning and navigation in forests. The Quercus variabilis broad-leaved forest in Jiufeng National Forest Park and the Sabina Coniferous forest in Dongsheng Bajia forest farm were selected as the study area. A Sanding T-23 multi-frequency three-constellation receiver and a u-blox NEO-M8T multi-constellation receiving module were used for continuous observation under the forest canopy. Compared with T-23, the u-blox NEO-M8T is much lighter and more flexible in the forest. The BDS/GPS-integrated positioning algorithm for forests was constructed by temporally and spatially unifying the satellite systems and using a reasonable observed value weighting method. Additionally, the algorithm is also written into the RTKLIB software to calculate the three-dimensional (3D) coordinates of the forest observation point in the World Geodetic System 1984 (WGS-84) coordinate system. Finally, the results were compared with the positioning results obtained using GPS alone. The experimental results indicated that, compared with GPS positioning, there were 13–27 visible satellites available for the BDS/GPS-integrated positioning algorithm for forests, far more than the satellites available for the GPS positioning algorithm alone. The Position Dilution of Precision (PDOP) values for the BDS/GPS-integrated positioning ranged from 0.5 to 1.9, lower than those for GPS positioning. The signal noise ratio (SNR) of the BDS/GPS-integrated satellite signals and GPS satellite signals were both in the range of 10–50 dB-Hz. However, because there were more visible satellites for the BDS/GPS-integrated positioning, the signals from the BDS/GPS-integrated satellites were stronger and had a more stable SNR than those from the GPS satellites alone. The results obtained using the BDS/GPS-integrated positioning algorithm for forests had significantly higher theoretical and actual accuracies in the X, Y and Z directions than those obtained using the GPS positioning algorithm. This suggests that the BDS/GPS-integrated positioning algorithm can obtain more accurate positioning results for complex forest environments.

Construction and Accuracy Analysis of a BDS/GPS-Integrated Positioning Algorithm for Forests

Fei Yan, Xueqi Hu, Lidan Xu, Yongrui Wu

Abstract

The objective of this study was to construct a BeiDou navigation satellite system (BDS)/global positioning system (GPS)-integrated positioning algorithm that meets the accuracy requirement of forest surveys and to analyze its accuracy to provide theoretical and technical support for accurate positioning and navigation in forests. The Quercus variabilis broad-leaved forest in Jiufeng National Forest Park and the Sabina Coniferous forest in Dongsheng Bajia forest farm were selected as the study area. A Sanding T-23 multi-frequency three-constellation receiver and a u-blox NEO-M8T multi-constellation receiving module were used for continuous observation under the forest canopy. Compared with T-23, the u-blox NEO-M8T is much lighter and more flexible in the forest. The BDS/GPS-integrated positioning algorithm for forests was constructed by temporally and spatially unifying the satellite systems and using a reasonable observed value weighting method. Additionally, the algorithm is also written into the RTKLIB software to calculate the three-dimensional (3D) coordinates of the forest observation point in the World Geodetic System 1984 (WGS-84) coordinate system. Finally, the results were compared with the positioning results obtained using GPS alone. The experimental results indicated that, compared with GPS positioning, there were 13–27 visible satellites available for the BDS/GPS-integrated positioning algorithm for forests, far more than the satellites available for the GPS positioning algorithm alone. The Position Dilution of Precision (PDOP) values for the BDS/GPS-integrated positioning ranged from 0.5 to 1.9, lower than those for GPS positioning. The signal noise ratio (SNR) of the BDS/GPS-integrated satellite signals and GPS satellite signals were both in the range of 10–50 dB-Hz. However, because there were more visible satellites for the BDS/GPS-integrated positioning, the signals from the BDS/GPS-integrated satellites were stronger and had a more stable SNR than those from the GPS satellites alone. The results obtained using the BDS/GPS-integrated positioning algorithm for forests had significantly higher theoretical and actual accuracies in the X, Y and Z directions than those obtained using the GPS positioning algorithm. This suggests that the BDS/GPS-integrated positioning algorithm can obtain more accurate positioning results for complex forest environments.

Keywords: GNSS positioning, forest, accuracy assessment

1. Introduction

Forestry, as a principal part of ecological civilization construction, plays important ecological, economic, social and cultural roles. To better realize the functions of forestry, maintain healthy development of forests and informatize forest management, precision is of particular importance. Global navigation satellite systems (GNSS) play a tremendous role in fine forest monitoring and management. GNSS can provide real-time, accurate positioning information, occurrence area information and navigation services for core forestry operations, such as forest resource surveys, forest fire prevention and emergency planning, and plant disease and insect pest monitoring, prevention and control. As a result, GNSS are a central topic in global forest survey research (Bettinger et al. 2019, Becker et al. 2017, Brach et al. 2014, Liu et al. 2017).

However, to date, most research on the application of satellite navigation in forestry has relied on the U.S. Navigation System with Timing and Ranging/Global Positioning System (GPS). GPS technology is relatively mature with an extensive user base and has established an industrial leader role globally (Yang et al. 2016). Wing et al. (2006) investigated the accuracy and reliability of the GPS in a variety of forest conditions in the forest region of Oregon, a western state of the U.S. Sawaguchi et al. (2005) studied the relations between cedar trees, needle pines and broad-leaved forests and the signal-to-noise ratio (SNR) of GPS-receiving satellites. Lai of Sichuan Agricultural University (2005) studied the factors affecting the accuracy of the GPS, concluding that the accuracy of the GPS is relatively closely related to the forest stand canopy density, the dominant forest stand tree species, the observation period and the date. Feng et al. (2004) and Tan et al. (2008), of Beijing Forestry University, extensively investigated positioning results obtained using differential GPS technology in various forest canopy and mountainous terrain conditions, summarizing the effects of land type and terrain on the accuracy of the GPS. Zhang et al. (2014) studied the effects of site condition and positioning mode on the accuracy of the GPS. In 2015, Zhang et al. (2015) conducted an experiment on the use of a satellite-ground-combined positioning approach for high-canopy-density forest stands. Kaartinen et al. (2015) tested the accuracy of various instruments utilizing global satellite navigation systems (GNSS) in motion under forest canopies of varying densities. Overall, because of the unique satellite constellation network of the single GPS and the complexity of forestlands, research on the application of satellite positioning in forest regions has met with problems, such as poor satellite signal, loss of lock on the signal and strong multipath effects in forest environments, resulting in relatively low positioning accuracy or an inability to conduct satellite positioning in forest regions (Feng et al. 2017, Zhang et al. 2015, Li et al. 2014, Xu et al. 2013). Consequently, the positioning accuracy requirement for some forest operations cannot be met.

In addition to GPS, there are Galileo, GLONASS, and Beidou navigation systems, QZSS and IRNSS satellite signal enhancement systems. Improving positioning accuracy through compatibility and interoperability between systems is a current research hot spot. However, more than two system choices could greatly increase the complexity of data solution, i.e. the system options are very important (Zrinjski et al. 2019, Montenbruck et al. 2017, Solarić et al. 2018). In the Asia-Pacific region, with the completion of Beidou, the combination of BDS and GPS is a suitable choice (Zhu et al. 2018). The BeiDou navigation satellite system (BDS) was independently developed by China. It is capable of providing global, continuous, all-weather, high-accuracy, autonomous navigation and positioning services. The BDS consist of 46 satellites operating in three types of satellite orbits, differing in altitude, i.e., medium Earth orbits, geostationary orbits (GEO) and inclined geosynchronous orbits (IGSO) (Liu 2011, 2013, Li et al. 2005). The BDS network was completed in June 23, 2020. Because of the unique high-orbit (GEO and IGSO) satellites of the BDS and the resultant local enhancement for the Asia-Pacific region, there has been a significant increase in the number of satellites and in the quality of satellite signals that can be received within the entire Asia-Pacific region. Even in inclement environments, where tall buildings, trees and canyons cause obstructions, more high-elevation-angle visible satellites in the BDS can be received compared to those received by the GPS (Yan et al. 2019, Liu 2003, 2015). The gaps left by a relatively small number of satellites have been eliminated and multipath effects have been weakened. As a result, the accuracy and precision of positioning under the forest canopy have improved. The BDS provides an important technical means for in-depth research on satellite positioning and navigation for forestry.

In view of the above, a BDS/GPS-integrated positioning algorithm that meets the accuracy requirement of forest resource surveys is constructed using the BDS/GPS-integrated single-point positioning algorithm for forests as the main technical method, the Jiufeng National Forest Park and Dongsheng Bajia forest farm in the Haidian District of Beijing, China as the study area, and a T-23 multi-frequency three-constellation receiver and a u-blox NEO-M8T multi-constellation receiving module as data acquisition equipment. Additionally, the integrated positioning algorithm is also written into the RTKLIB software to calculate the three-dimensional (3D) coordinates of the forest observation point in the World Geodetic System 1984 (WGS-84) coordinate system. Moreover, the results are compared with the positioning results obtained using GPS alone to investigate the effects of high-canopy density forest environments on the number of visible satellites (NVS), the position dilution of precision (PDOP) value and the SNR. Furthermore, the accuracy of the results is also evaluated to verify the advantages of multi-system positioning algorithms in complex forest environment conditions. The results obtained from this study provide a theoretical basis and technical support for research on the application of the BDS in forestry and contribute to the industrialization of the BDS in the forestry field.

2. Research Materials

2.1 General Information on the Study Area

The selected broad-leaved forest experimental area (Table 1) is located in the Jiufeng National Forest Park in the Haidian District of Beijing, China. There are rolling mountain ranges and a complex terrain in the park, where the highest peak has an elevation of 1153 m and the main peak is 465 m. Broad-leaved deciduous forests are the primary vegetation type. Below the low-mountain level (800 m), there are large areas of man-made forests (Platycladus orientalis, Chinese pine, Quercus variabilis, Robinia pseudoacacia, etc.) dotted with Prunus armeniaca and Vitex negundo shrubs. In the mid-mountain range (800–1100 m), oak-pine forests are the native vegetation and Chinese pine and larch are the primary tree species; the forest coverage is up to 96.4%. Situated on the east side of the park (40°03'49"N, 116°05'51"E), the selected study area has an elevation of approximately 130 m, a slope gradient of approximately 22°, and a southeastern slope aspect, as shown in Fig. 1 (left). Data were collected in mid- and late June 2018. Before the experiment began, it had rained multiple times in the forest region. As a result, the soil was soft, which helped to set up the equipment. Quercus variabilis, mostly over 30 years old, is the primary tree species in the study area. Other stand parameters are detailed in Table 1.

The Coniferous forest experimental area (Table 1) is located in the Dongsheng Bajia forest farm, also in the Haidian District of Beijing, China, covering an area of 101.45 ha. There are abundant vegetation species in the experimental area. The coniferous forest is mainly composed of Sabina, Pinus tabulaeformis and Chinese fir, and the broad-leaved forest is dominated by poplar and willow. Situated on the north side of the park (40°01′19″N，116°20′30″E), there are no high-rise buildings around the selected study area, and satellite signals can be effectively prevented from being blocked, as shown in Fig. 1 (right). Data were collected in mid- and late April 2019. Sabina chinensis is the primary tree species in this study area and the stand parameters are also detailed in Table 1.

Fig. 1 The study area of broad-leaved forest and Coniferous forest

Table 1 Mean stand characteristics of the study area

2.2 Experimental Equipment and Calculation Software

The equipment used in the study included a Sanding T-23 multi-frequency three-constellation receiver (Fig. 2a), a u-blox NEO-M8T multi-constellation receiving module (Fig. 2b), a tripod, a DBH girthing tape, and a wooden post. The software used in the study included GAMIT (Version 10.61) and RTKLIB (Version 2.4.3). GAMIT was used in the high-accuracy calculation of the 3D coordinates of the forest observation point in the WGS-84 coordinate system to obtain their true values with positioning accuracy in the centimeter range. RTKLIB was used in the programming of the BDS/GPS-integrated algorithm, the real-time single-point positioning and the accuracy analysis. Both software packages support multiple GNSS standards and construction of precise positioning algorithms, can achieve real-time GNSS positioning and post-processing calculations, and are commonly used to study GNSS-integrated algorithms (Paziewski et al. 2017).

Fig. 2 Experimental instrument and solving software

3. Methods

3.1 Observation Point Positioning Data Collection

An observation point at the center of the forest stand was set up to sufficiently represent its high-canopy density. A wooden post was driven into the forest ground to a depth of 10–15 cm. An observation point was set up and marked at the top of the wooden post, which was ensured to be firm. The antennas of a u-blox satellite receiver were placed very close to the observation point. The phase center of the antenna and the observation point were ensured to be on the same plumb line. In addition, a Sanding T-23 receiver was set up at the observation point and started to receive the signal after being strictly centered and leveled (Fig. 3). No less than 200 observation epochs were measured, and the epoch interval was 15 seconds. To obtain the high-accuracy calculation of the 3D coordinates of the forest observation point，the observation lasted for at least 4 hours.

Fig. 3 Observation point setting

3.2 Study and Construction of a BDS/GPS-Integrated Algorithm for Forests

The established BDS/GPS-integrated positioning and observation model for forests in this study is as follows:

(1)

Where:

ri observed pseudo range between satellite and receiver

B, G BDS and GPS, respectively

i satellite number

si geometric distance between satellite and receiver phase center

c speed of light

receiver clock error, consisting of the clock error of the receiver of each BDS and GPS

satellite clock error

ionospheric delay error

tropospheric delay error

and the satellite ephemeris error, multipath error and measurement noise are omitted in the equation (Li et al. 2016, Gao et al. 2014). The geometric distance between the satellite and the receiver is as follows:

(2)

Where:

coordinates of satellite i

coordinates of observation point

A linearized observation equation is obtained by Taylor series expansion at the approximate coordinates of the observation point. By reorganizing this equation, an indirect adjustment error equation is obtained as follows:

(3)

(4)

(5)

Where:

observed residual

observation point coordinate correction values

constant term

distance between the approximate location of the forest observation point and satellite i (Li et al. 2016, Qu et al. 2015)

The BDS/GPS-integrated forest positioning error equation contains 5 unknown parameters. According to the method of least squares, it is necessary to observe no fewer than 5 satellites simultaneously to calculate the location. Assuming that the satellite receiver receives signals from m number of GPS satellites and n number of BDS satellites at a certain observation time t, according to equation (3), the integrated positioning error equation can be expressed as follows:

(6)

Where:

V observed residual

B coefficient matrix

x corrected coordinates and receiver clock error

L constant matrix

The corrected 3D coordinates of the observation point can be calculated using the classic method of least squares as follows:

(7)

Where:

P weight of BDS and GPS satellite-observed values

Finally, the coordinates of the forest observation point in the WGS-84 coordinate system are obtained as follows:

(8)

3.2.1 Temporal and Spatial Unification

The time-space references used in satellite navigation consist of a time reference and a coordinate reference. The GPS uses the GPS time (GPST) and the WGS-84 coordinate system, whereas the BDS uses the BeiDou time (BDT) and the China Geodetic Coordinate System 2000 (CGCS-2000) coordinate system. If these two systems are integrated for positioning, their time-space references must be unified.

With respect to the spatial coordinate reference, the WGS-84 and the CGCS-2000 coordinate systems are the same in the definitions of the origin of the coordinates, the scale, the coordinate axes and changes in the directions of the coordinate axes with time. They only differ slightly in oblateness, whose effect is on the sub-millimeter level. Therefore, the difference in oblateness between the two systems is negligible. According to the accuracy requirement of forest surveys, the spatial coordinate references for the two systems can be viewed as equivalent (Paziewski et al. 2017).

With respect to the time reference, the BDT and GPST are both atomic time and count in seconds and numbers of weeks. However, there is a difference in the start of the time counting between the BDT and GSPT. There is a starting time difference of 1356 weeks between the BDT and GPST. The conversion equation for the BDT and GPST is thus as follows (Yang et al. 2016):

(9)

Where:

tGPS GPS time

tBDS BeiDou time

3.2.2 Observation Weight Determination

When positioning in a forest, the elevation angle between the observation point and each satellite has a relatively significant impact on the positioning accuracy of the observed value. Elevation angle weight determination refers to the determination of the weight of the observed value based on the elevation angle of each satellite used in observation. Different elevation angles lead to different ranging code transmission paths, which in turn lead to different errors. When positioning in a forest, if the elevation angle of a satellite is too small, its signals can be very easily blocked by the tree canopy. Correspondingly, the number of signal transmission paths in the forest increases and the ionospheric delay, tropospheric delay and multipath effects on the ranging codes transmitted by the satellite may increase, resulting in an increase in the measurement error in the distance between the station and the satellite. Therefore, the use of a reasonable elevation angle weight can help to obtain more accurate positioning results (Yang et al. 2016). The following elevation angle weight determination model is used in this study:

(10)

Where:

E elevation angle corresponding to the observed value

Observation weight determination was the key of this research and was well integrated and written in RTKLIB software.

3.3 Accuracy Evaluation

3.3.1 Analysis of the Characteristic Factors for Positioning Accuracy

The number of visible satellites (NVS), the PDOP value and SNR are important metrics that characterize satellite positioning accuracy. Under normal circumstances, a larger number of visible satellites forms a better geometric pattern formed by the satellites and the observation station, resulting in a lower PDOP value and a higher SNR, i.e., more accurate final positioning results. Therefore, when comparing BDS/GPS-integrated and GPS positioning in forests, it is necessary to first compare the aforementioned three factors. In this study, the RTKLIB software was used to extract the NVS, PDOP value, SNR data and patterns in the data.

3.3.2 Accuracy Analysis of Positioning Results

The experimental results are comparatively analyzed from the perspectives of theoretical and actual accuracies.

3.3.2.1 Theoretical Accuracy

Based on equations (6)–(8), the error expression and the equation for calculating accuracy are obtained as follows:

(11)

Where:

Q inverse matrix of the unknown coefficient matrix of the normal equation, which is also referred to as the cofactor matrix or the accuracy matrix

D medium error of the corresponding unknown parameters (X, Y, Z and receiver clock error), i.e., theoretical accuracy

unit weight MSE

n number of observed values

t = 5

V residual error of observed values

3.3.2.2 Actual Accuracy

The actual conditions involved in positioning in the field are complex; therefore, positioning results are often affected by a combination of accidental and residual system errors. This means that there may be a certain difference between the actual and theoretical accuracies. Therefore, it is necessary to estimate the actual accuracy of the positioning results obtained using the BDS/GPS-integrated algorithm. By comparing the 3D coordinates of the observation point obtained using the BDS/GPS-integrated algorithm and the GPS algorithm, with the true coordinates of the observation point obtained from a high-accuracy calculation of the carrier phase using the GAMIT software (observed value – true value), changes in error for each epoch can be obtained. In addition, the actual positioning accuracy is evaluated using equation (12) as follows (Balenović et al. 2018, Xu et al. 2020):

(12)

Where:

, , actual positioning accuracies in X, Y and Z directions, respectively

n number of observed values

, , true coordinates of the receiver in X, Y and Z directions, respectively

, , coordinates calculated using BDS/GPS-integrated positioning algorithm or GPS positioning algorithm

4 Results and Discussion

4.1 NVS

The number of visible satellites (NVS) that could be observed during the whole period of positioning under the forest canopy and the corresponding signal quality were recorded using the RTKLIB software. In the broad-leaved forest area, for the BDS/GPS-integrated positioning algorithm, there was a total of 15–23 satellites available for positioning and observation during the observation period, including 11 GPS satellites and 12 BDS satellites (Figs. 4 a and b). In the Coniferous forest, there was a total of 13–27 satellites available for positioning and observation including 11 GPS satellites and 16 BDS satellites (Figs. 4 c and d).

Figs. 5 a–d show the continuity of the received signals, a–b was the signals received situation in broad-leaved forest, while c–d was the signals received situation in Coniferous forest. A continuous reception of signals (e.g. those from the G04, G14, C09 in Fig. 5 a–b and G26, C06, C09 satellites, as shown in Fig. 5 c–d indicates relatively high signal quality. The reception of signals from the G23, C03 satellites in a–b and the G23, C01, C21 satellites shown in c–d were intermittent, suggesting that losses of the satellite lock might have occurred during the observation period. In addition, there were more than four visible satellites for the BDS/GPS-integrated positioning and the GPS positioning, meeting the requirement of the necessary number of observed values for the calculation of coordinates. However, there were more visible satellites for the BDS/GPS-integrated positioning than for the GPS positioning, suggesting that the integrated positioning algorithm can have an excess number of observed values for the calculation when used for positioning in forests and, consequently, produce more accurate positioning results.

Fig. 4 Satellite visible number comparison

Fig. 5 Signal continuity comparison

4.2 PDOP Value

The PDOP value refers to the 3D PDOP factor, which is an important metric that reflects the strength of the geometric pattern formed by the observation point and the satellites at the observation time. It is also used to evaluate the positioning capacity of a satellite navigation system (Cai et al. 2011, Li et al. 2014). As demonstrated in Figs. 6 a–d, the PDOP values for BDS/GPS-integrated positioning and observation ranged from 0.7 to 1.8 in broad-leaved forest as shown in Fig. 6 a–b, while the PDOP value changed steadily ranging from 1.1 to 1.9 in coniferous forest, whereas the PDOP values for GPS positioning ranged from 1.0 to 2.1 in broad-leaved forest and from 1.6 to 2.5 in coniferous forest. The PDOP values for BDS/GPS-integrated positioning were lower than those for GPS positioning, suggesting that the use of BDS/GPS-integrated positioning in forests of different vegetation types significantly increased the strength of the geometric pattern formed by the observation station and the satellites and thereby improved the accuracy of the positioning results.

Fig. 6 PDOP comparison

4.3 SNR

SNR refers to the ratio of the strength of a signal received by a receiver to the strength of the noise. Generally speaking, a higher SNR implies less noise in the signal, a higher satellite signal quality, and more accurate positioning results (Khodabandeh et al. 2016, Gao et al. 2015). As demonstrated in Figs. 7 a–d, in the two forest experimental areas of two different types, the SNRs of the satellite signals received during the observation period were in the range of 10–50 dB and mostly in the range of 30–45 dB, suggesting relatively strong satellite signals. In addition, as the BDS/GPS-integrated positioning algorithm significantly increased NVS for positioning under the forest canopy, the SNRs of the signals from the BDS/GPS satellites were more stable than those from the GPS satellites alone and the accuracy of the BDS/GPS-integrated positioning algorithm was higher than that of the GPS positioning algorithm.

Fig. 7 SNR comparison

4.4 Theoretical Accuracy

Two hundred continuous epochs with stable signals were selected for the accuracy calculation. The theoretical accuracies of the BDS/GPS-integrated positioning algorithm for forests and the GPS positioning algorithm were obtained. Table 2 summarizes the theoretical accuracies of the results obtained using the two positioning algorithms, based on equation (11).

The theoretical accuracy analysis shows that, no matter in which vegetation types, there was no significant difference between the results obtained using the BDS/GPS-integrated positioning and GPS positioning algorithms in the X direction, but the theoretical accuracy of the results obtained using the BDS/GPS-integrated positioning algorithm were smaller than those of the results obtained using the GPS positioning algorithm in the Y and the Z directions. These results are in agreement with the network characteristics of the BDS.

Table 2 Theoretical accuracy of positioning results

4.5 Actual Accuracy

The absolute error of the results for each epoch in the X, Y and Z directions were determined by comparing the calculated data used in Section 4.4 with the true coordinates of the observation point obtained by Sanding T-23 receiver and GAMIT software, as shown in Figs. 8 to 11.

Fig. 8 Absolute error of BDS/GPS-integrated positioning in X, Y and Z directions in broad-leaved forest

Fig. 9 Absolute error of GPS positioning in X, Y and Z directions in broad-leaved forest

Fig. 10 Absolute error of BDS/GPS-integrated positioning in X, Y and Z directions in coniferous forest

Fig. 11 Absolute error of GPS positioning in X, Y and Z directions in coniferous forest

The data show that, in the broad-leaved forest, the error of the results obtained using the BDS/GPS-integrated positioning algorithm for forests fluctuated around 0 m and its extreme value did not exceed 8 m. Its absolute value in the Z direction exceeded 20 m. In the coniferous forest, the error of the results obtained using the BDS/GPS-integrated positioning algorithm also fluctuated around 0 m and its extreme value in X direction did not exceed 4 m, in Y and Z directions, the extreme value was lower than 10 m. It showed that, in forests of different vegetation types, the error of the results obtained using the GPS positioning algorithm alone was significantly greater than that of the results obtained using the BDS/GPS-integrated positioning algorithm.

The actual accuracies of the two algorithms for positioning under the forest canopy were analyzed using equation (12). Table 3 summarizes the actual accuracies in the X, Y and Z directions.

Table 3 Actual accuracy of positioning results

A comparison of the actual accuracies of the four sets of data shows that the error in positioning in the high-canopy density forest stand using the GPS algorithm alone was relatively large; the largest error in the Z direction was 7.261 m. The comparison showed that the results obtained using the BDS/GPS-integrated positioning algorithm were closer to the true values. The components of the error of the results obtained using the BDS/GPS-integrated positioning algorithm in the X directions were relatively close to one another, whereas the component of the error in the Y and Z direction was relatively large. The components of the error of the results obtained using the BDS/GPS-integrated positioning algorithm in the three directions were significantly smaller than those obtained using the GPS positioning algorithm.

5. Conclusion

Using the forestlands of typical vegetation types in Jiufeng National Forest Park and Dongsheng Bajia forest farm in the Haidian District of Beijing, China as the study area, a BDS/GPS-integrated positioning algorithm was developed by creating interoperability between the BDS and the GPS and setting the system weights. The algorithm was also written into the RTKLIB software to obtain positioning results by direct calculation, while receiving signals from the two satellite systems. The results obtained from this study show the following. Under the obstructing conditions in the high-canopy dense forest, the NVS of the BDS/GPS-integrated system was significantly higher than that of the GPS alone. The PDOP values for integrated positioning were better than those for GPS positioning. The SNRs of the signals from the BDS/GPS-integrated positioning satellites were higher and more stable than those from the GPS positioning satellites. The comparison of the theoretical and actual accuracies of the results of positioning under the forest canopy shows that the accuracy of the BDS/GPS-integrated single-point positioning algorithm was notably superior to that of the GPS positioning algorithm. Additionally, it can be inferred that since the BDS/GPS-integrated positioning algorithm can significantly increase NVS and the high-orbit (GEO and IGSO) satellites of the BDS can reduce the multipath error, the geometric patterns of the satellites is strengthened when compared to the GPS. Therefore, the accuracy advantage of the BDS/GPS-integrated positioning algorithm should become more pronounced in regions with higher forest canopy densities or more complex terrains.

This study investigated two high-canopy density forest regions in Beijing. Further experimental verification is required for other complex environments, such as forestlands and canyons with various vegetation types, terrain and topographical conditions. As the BDS is continuously being improved, it can be expected that multi-system-integrated positioning approaches could provide more stable and accurate positioning, time and navigation services for monitoring and management of forest resources in forest environments with higher canopy densities and more complex terrains.

Acknowledgements

Financial support for this study was supported by National Natural Science Foundation of China (Project No. 31800468). We are grateful to the undergraduate students and staff of the Laboratory of Forest Management and »3S« technology and Beijing Forestry University.

6. References

Balenović, I., Gašparović, M., Milas, A.S., Berta, A., Seletković, A., 2018: Accuracy assessment of digital terrain models of lowland pedunculate oak forests derived from airborne laser scanning and photogrammetry. Croatian Journal of Forest Engineering 39(1): 117–128.

Becker, R.M., Keefe, R.F., Anderson, N.M., 2017: Use of real-time GNSS-RF data to characterize the swing movements of forestry equipment. Forests 8(2): 44. https://doi.org/10.3390/f8020044

Bettinger, P., Merry, K., Bayat, M., Tomaštík, J., 2019: GNSS use in forestry – A multi-national survey from Iran, Slovakia and southern USA. Computers and Electronics in Agriculture 158: 369–383. https://doi.org/10.1016/j.compag.2019.02.015

Brach, M., Zasada, M., 2014: The effect of mounting height on GNSS receiver positioning accuracy in forest conditions. Croatian Journal of Forest Engineering 35(2): 245–253.

Chen, G.Y., 2004: Application of global positioning in forestry. Journal of Beijing Agricultural College 19(3): 46–48.

Cai, C.S., Dai, W.J., Kuang, C.L., 2011: Calculation and analysis of PDOP for combined GPS/GLONASS systems. Bulletin of Surveying and Mapping (11): 5–7.

Feng, Z.K., Yan, F., Mohammad, R.U., Dang, Y., 2017: Developing a volume model using south NTS-372R total station without tree felling in a Populus canadensis moench plantation in Beijing, China. Croatian Journal of Forest Engineering 38(1): 127–150.

Gao, W., Gao, C.F., Pan, S., Wang, D., Wang, S., 2015: Single-epoch positioning method in network RTK with BDS triple-frequency. Acta Geodaetica et Cartographica Sinica 44(6): 641–648. https://doi.org/10.11947/j.AGCS.2015.20140308

Isao, S., Yuko, S., Shiro, T., 2005: A study of the effects of stems and canopies on the signal to noise ratio of GPS signals. Journal of Forest Research 10(5): 395–401. https://doi.org/10.1007/s10310-005-0158-7

Jia, H.X., Xing, H.Y., Jia, H.S., 2009: Application of GPS technology in forestry. Inner Mongolia Forestry Investigation and Design 32(3): 78–80．

Jiang, Z.X., Zhang, H., Guo, Z.H., 2010: The method for calculation of internal accord accuracy in airborne gravity survey. Geophysical & Geochemical Exploration 34(5): 672–676.

Kaartinen, H., Hyyppä, J., Vastaranta, M., Kukko, A., Jaakkola, A., Yu, X., Pyörälä, J., Liang, X., Liu, J., Wang, Y., Kaijaluoto, R., Melkas, T., Holopainen, M., Hyyppä, H., 2015: Accuracy of kinematic positioning using global satellite navigation systems under forest canopies. Forests 6(9): 3218–3236. https://doi.org/10.3390/f6093218

Khodabandeh, A., Teunissen, P., 2016: PPP-RTK and inter-system biases: the ISB look-up table as a means to support multi- system PPP-RTK. Journal of Geodesy 90(9): 837–851. https://doi.org/10.1007/s00190-016-0914-9

Lai, J.M., 2005: Application research of GPS technology in positioning and area measurement in forest resources survey in low mountain areas of western Sichuan. Chengdu, Sichuan Agricultural University.

Li, B., Li, Z.H., Liu, W.K., 2013: Simulation analysis of combined system of COMPASS/ GPS / GLONASS in china region. Journal of Geodesy and Geodynamics 33(6): 94–97.

Li, P., Zhang, X., 2014: Modeling and performance analysis of GPS/GLONASS/BDS precise point positioning. China Satellite Navigation Conference (CSNC) 2014 Proceedings, Volume III. Springer Berlin Heidelberg, 251–263.

Li, Y.N., Liu, L.M., Cui, L. Huang, L., Zhao, Z., Huang, X., 2014: Effects of canopy openness on positioning availability and initialization time of GNSS RTK in forests. Scientia Silvae Sinicae 50(2): 78–84. https://doi.org/10.11707/j.1001-7488.2014021

Li, Z.H., Huang, J.S., 2016: GPS Measurement and Data Processing. Wuhan, Wuhan University Press.

Liu, J., Hyyppä, J., Yu, X., Jaakkola, A., Kukko, A., Kaartinen, H., Zhu, L., Liang, X., Wang, Y., Hyyppä, H., 2017: A novel GNSS technique for predicting boreal forest attributes at low cost. IEEE Transactions on Geoscience and Remote Sensing 55(9): 4855–4867. https://doi.org/10.1109/TGRS.2017.2650944

Liu, J.Y., 2003: Principle and method of GPS satellite navigation and positioning. Beijing, Science Press, 399–401.

Liu, J.Y., 2011: Some suggestion for enhancing compass navigation performance. Digital Communication World 2: 31–35.

Liu, J.Y., 2013: Status and development of the Beidou navigation satellite system. Journal of Telemetry, Tracking and Command 34(3): 1–8.

Liu, J.Y., 2015: Development thinking on GNSS signal receivers-transmitting/receiving issue (14) of GNSS navigation signals. Digital Communication World 10: 1–9.

Michael, G.W., Aarone, E., 2006: Standard and real-time accuracy and reliability of a mapping-grade GPS in a coniferous western Oregon forest. Western Journal of Applied Forestry 21(4): 222–227. https://doi.org/10.1093/wjaf/21.4.222

Montenbruck, O., Steigenberger, P., Prange, L., Deng, Z., Zhao, Q., Perosanz, F., Romero, I., Noll, C., Stürze, A., Weber, G., Schmid, R., MacLeod, K., Schaer, S., 2017: The multi-GNSS experiment (MGEX) of the international GNSS service (IGS) – achievements, prospects and challenges. Advances in space research 59(7): 1671–1697. https://doi.org/10.1016/j.asr.2017.01.011

Yuzao, N., Zhongke, F., Xiangrui, D., Xiaokun, W., 2004: Effects of crown cover and hill country on RTD GPS positioning precision and corresponding reducing countermeasures. Journal of Beijing Forestry University 26(2): 1–5.

Paziewski, J., Sieradzki, R., 2017: Integrated GPS+BDS instantaneous medium baseline RTK positioning: signal analysis, methodology and performance assessment. Advances in Space Research 60(12): 2561–2573. https://doi.org/10.1016/j.asr.2017.04.016

Qu, L., Zhao, Q., Guo, J., Wang, G., Guo, X., Zhang, Q., Jiang, K., Luo, L., 2015: BDS/GNSS real-time kinematic precise point positioning with un-differenced ambiguity resolution. China Satellite Navigation Academic Annual Meeting, 13–29.

Solarić, N., Solarić, M., Bogdanovski, Z., Dimeski, S., 2018: Three days before the earthquake in skopje there was a compression of the earth's crust. Geodetski list 72(1): 15–35.

Tan, W., Wang, K.L., Luo, X., Wang, Z.-J., 2008: Positioning precision with handset GPS receiver in different stands. Journal of Beijing Forestry University 30(3): 163–167.

Xu, W., Su, Z., Feng, Z., Xu, H., Jiao, Y., Yan, F., 2013: Comparison of conventional measurement and LiDAR-based measurement for crown structures. Computers and Electronics in Agriculture 98: 242–251. https://doi.org/10.1016/j.compag.2013.08.015

Xu, W., Huang, S., Wu, C., Xiong, Y., Wang, L., Lu, N., Kou, W., 2020: The pruning phenological phase-based method for extracting tea plantations by field hyperspectral data and Landsat time series imagery. Geocarto International: 1–21. https://doi.org/10.1080/10106049.2020.1801859

Yan, F., Wang, C.B., Wu, Y.R., Lin., L., Feng, Z., 2019: Algorithm Implementation and Precision Analysis of Forest BDS/GPS Combined Positioning. Transactions of the Chinese Society for Agricultural Machinery 50(4): 221–227. https://doi.org/10.6041/j.issn.1000-1298.2019.04.025

Yang, Y.X., Lu, M.Q., Han, C.H., 2016: Some notes on interoperability of GNSS. Acta Geodaetica et Cartographica Sinica 45(3): 253–259. https://doi.org/10.11947/j.AGCS.2016.20150653

Zhang, H.C., Zhou, H.P., Zheng, J.Q., 2014: Effects of the site condition and positioning mode on GPS performance for precision forestry. Journal of Nanjing Forestry University (Natural Sciences Edition) 38(6): 72–76.

Zhang, J.Q., Pan, L., Wang, S.G., 2012: Photogrammetry. Wuhan, Wuhan University Press, 90–93.

Zhang, X.Q., Zhang, X., Yu, X.W., Deng, G., Ma, Y., 2015: Study on satellite-ground combining positioning in high canopy stands. Journal of Nanjing Forestry University (Natural Sciences Edition) 39(3): 25–28.

Zhao, Y.D., Tu, J.Y., 2018: Intelligent surveying and mapping system of forest inspection based on Beidou Satellite. Transaction of the Chinese Society for Agricultural Machinery 49(7): 177–185. https://doi.org/10.6041/j.issn.1000-1298.2018.07.022

Zhu, N., Marais, J., Bétaille, D., Berbineau, M., 2018: GNSS position integrity in urban environments: A review of literature. IEEE Transactions on Intelligent Transportation Systems 19(9): 2762–2778. https://doi.org/10.1109/tits.2017.2766768

Zrinjski, M., Barković, Đ., Matika, K., 2019: Development and modernization of GNSS. Geodetski list 73(1): 45–65.

© 2020 by the authors. Submitted for possible open access publication under the
terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Authors’ addresses:

Fei Yan *

e-mail: yf-perfect@163.com

Xueqi Hu

e-mail: 969616662@qq.com

Lidan Xu

e-mail: xlidan@bjfu.edu.cn

Yongrui Wu

e-mail: 1427663547@qq.com

Precision Forestry Key Laboratory of Beijing

College of Forestry

Beijing Forestry University

100083 Beijing

P.R. CHINA

* Corresponding author

Received: July 21, 2020

Accepted: October 10, 2020

Original scientific paper