# ABSTRACT:

Keywords:
RTK; SBAS; accuracy; differential correction

# Introduction

The use of Global Navigation Satellite Systems (GNSS) in agriculture has evolved, enabling a revolution in georeferenced data collection, which is performed faster, more accurately, and less costly. Precise positioning is demanded in several agricultural applications (Kabir et al., 2016Kabir, M.S.N.; Song, M.Z.; Sung, N.S.; Chung, S.O.; Kim, Y.J.; Noguchi, N.; Hong, S.J. 2016. Performance comparison of single and multi-GNSS receivers under agricultural fields in Korea. Engineering in Agriculture, Environment and Food 9: 27-35.). Accuracy of commercial GNSS receivers is one of the requirements from farmers and depends on operations to be performed in the field. Thus, the accuracy of its signals and its behavior in agricultural settings should be further investigated. Some challenges in agricultural operations involve the selection of an adequate source of differential correction signals, depending on the accuracy required, and the choice for the right option, regarding availability, practical aspects, and costs.

Some agricultural operations, such as auto-guidance technology for precision inter-row cultivation, require high accuracy, normally below 0.10 m, and it is only possible with the use of differential correction signals (Machado and Molin, 2011Machado, T.M.; Molin, J.P. 2011. Static and cinematic tests of GPS receivers. Revista Brasileira de Engenharia Agrícola e Ambiental 15: 981-988 (in Portuguese, with abstract in English).). Real-Time Kinematic (RTK) and Satellite-based Augmentation Systems (SBAS) with signal distributed by geostationary satellites stand out as the main options for users of real-time differential correction. The use of Network RTK can be an option to share costs without degrading positional accuracy (Bae and Kim, 2018Bae, T.S.; Kim, M. 2018. Performance analysis of Network-RTK techniques for drone navigation considering ionospheric conditions. Journal of Sensors 1: 1-8.). However, for remote agricultural areas in Brazil, there is little or no access to Internet networks, making Network RTK unavailable or costly.

Machado et al. (2010)Machado, T.M.; Molin, J.P.; Povh, F.P.; Salvi, J.V. 2010. Methodology for performance evaluation of GPS receiver for agricultural use in kinematic condition. Engenharia Agrícola 30: 121-129 (in Portuguese, with abstract in English). and Machado and Molin (2011)Machado, T.M.; Molin, J.P. 2011. Static and cinematic tests of GPS receivers. Revista Brasileira de Engenharia Agrícola e Ambiental 15: 981-988 (in Portuguese, with abstract in English). evaluated the performance of GNSS receivers installed on an agricultural vehicle, with circular paths in a constant trajectory, using a GNSS receiver with RTK signal as a reference for accuracy and precision calculations. Carballido et al. (2014)Carballido, J.; Perez-Ruiz, M.; Emmi, L.; Aguera, J. 2014. Comparison of positional accuracy between RTK and RTX GNSS based on the autonomous agricultural vehicles under field conditions. Applied Engineering in Agriculture 30: 361-366. and Kabir et al. (2016)Kabir, M.S.N.; Song, M.Z.; Sung, N.S.; Chung, S.O.; Kim, Y.J.; Noguchi, N.; Hong, S.J. 2016. Performance comparison of single and multi-GNSS receivers under agricultural fields in Korea. Engineering in Agriculture, Environment and Food 9: 27-35. evaluated the performance of multi-GNSS receivers in static and dynamic conditions in different agricultural sites.

Due to limited information on the performance of GNSS receivers and signals for static and dynamic tests, this study aims to advance knowledge on the methodological approach for evaluating the static and dynamic performance of GNSS receivers commonly used in agricultural operations, including statistical issues related to the number of replications.

# Materials and Methods

A supporting frame was allocated at the geographic coordinates 22°42’47.8” S, 47°37’44.9” W, altitude of 540 m. Ten different GNSS receivers, 0.50 m apart, were placed on top of the support frame with the geographic coordinates properly established (Figure 1).

Figure 1
Supporting frame to allocate GNSS receivers for the static test.

Field test procedures are used for evaluating satellite-based auto-guidance systems in agriculture (Sama and Stombaugh, 2014Sama, M.P.; Stombaugh, T.S. 2014. Performance evaluation of a tracking total station as a position reference for dynamic GNSS accuracy testing. Applied Engineering in Agriculture 30: 557-563.; Kim et al., 2016Kim, J.Y.; Kim, H.J.; Shim, S.B.; Park, S.H.; Kim, J.H.; Kim, Y.J. 2016. Development of automated guidance tracking sensor system based on laser distance sensors. Journal of Biosystems Engineering 41: 319-327.; Santos et al., 2017Santos, A.F.; Silva, R.P.; Tavares, T.O.; Ormond, A.T.S.; Rosalen, D.L.; Assis, L.C. 2017. Parallelism error in peanut sowing operation with auto-steer guidance. Revista Brasileira de Engenharia Agrícola e Ambiental 21: 731-736.). The American Society of Agricultural and Biological Engineers (ASABE) suggests ASABE X605 (2008American Society of Agricultural and Biological Engineers [ASABE]. 2008. X605 Satellite-Based Auto-Guidance Systems Testing During Straight and Level Travel. ASABE, St. Joseph, MI, USA.) for automated satellite guidance system testing during direct and level travel. ISO 12188-2 (2012)International Organization for Standardization [ISO]. 2012. ISO12188-2: Tractors and Machinery for Agriculture and Forestry – Test Procedures for Positioning and Guidance Systems in Agriculture. Part 2. Testing of Satellite-Based Auto-Guidance Systems During Straight and Level Travel. ISO, Geneve, Switzerland. is a standard for testing automatic satellite guidance systems based on ASABE X605, which defines a tracking sensor as an instrument or instrumented system designed to perform repeated horizontal distance measurements required for cross-track error calculations.

To perform the dynamic tests, a circular metal rail with 9.55 m radius, centered at the same point, similar to the static test, was placed at level ground, with a platform driven by an electric motor to perform the circular path (Figure 2A) of the GNSS receivers at a constant speed (1.0 m s−1). The platform stayed 2.0 m behind the electric motor (Figure 2B) to avoid possible electromagnetic interference to the antennas of GNSS receivers (Figure 2C).

Figure 2
Rail for the radial displacement of receivers (A); a platform with an electric motor (B) and platform for the allocation of GNSS receivers (C).

Figure 3
Calculation of GNSS receiver signal path error (A); location of receivers all aligned in the same distance from the reference coordinate (B).

The site where the tests were performed was not completely free of obstacles; however, it is considered adequate to conduct the tests, since it replicates real field conditions, especially the boundary of agricultural fields, with the presence of trees. The shading profile is represented by showing that obstacles did not exceed 30.0° (Figure 4).

Figure 4
Static shadowing profile with a starting point in the true North.

We used receivers with different specifications and accuracy levels, commonly applied in agriculture (Table 1). Group A comprises receivers without differential correction and Group B has receivers with SBAS or RTK correction. SBAS is a geostationary satellite system that improves the accuracy, integrity, and availability of GNSS signals. It does not require local infrastructure and offers competitive accuracy as a private service, commonly used for agricultural purposes in Brazil (Rovira-Más et al., 2015Rovira-Más, F.; Chatterjee, I.; Sáiz-Rubio, V. 2015. The role of GNSS in the navigation strategies of cost-effective agricultural robots. Computers and Electronics in Agriculture 112: 172-183.; Banville et al., 2018Banville, S.; Collins, P.; Lahaye, F. 2018. Model comparison for GLONASS RTK with low-cost receivers. GPS Solutions 22: 1-12.; Li et al., 2018Li, M.; Yuan, Y.B.; Wang, N.B.; Liu, T.; Chen, Y.C. 2018. Estimation and analysis of the short-term variations of multi-GNSS receiver differential code biases using global ionosphere maps. Journal of Geodesy 92: 889-903.).

Table 1
Equipment used in static and dynamic tests.

Software C7 (CR Campeiro) was installed on the cell phone to capture and save data. For data collection of L1 and L1/L2 receivers, portable computers were used with SST Field Rover II 7.13 software. For C/A code receivers, data was stored by the internal memory of the equipment.

All GNSS receivers were set to collect data at 1.0 s interval (1.0 Hz), logging for 4,200 s with five replications in the static test and seven replications for the dynamic test, keeping 3,600 s intervals between replications to allow complete reconfiguration of satellite constellations. GNSS involves signals from satellites subject to interference and provide positioning results with a range of common errors over time (Souza and Machado, 2016Souza, S.V.R.; Machado, A.M.L. 2016. Systematic Errors Calibration for MEMS Accelerometers. Boletim de Ciências Geodésicas 22: 835-850 (in Portuguese, with abstract in English).; Silva and Marques, 2016Silva, E.G.S.; Marques, H.A. 2016. PPP with GPS/GLONASS data Integration: Fundamentals Involved and Accuracy Assessment. Boletim de Ciências Geodésicas 22: 217-231 (in Portuguese, with abstract in English).; Ye et al., 2018Ye, S.; Zhao, L.; Song, J.; Chen, D.; Jiang, W. 2018. Analysis of estimated satellite clock biases and their effects on precise point positioning. GPS Solutions 22: 1-14.).

Before analyzing the data, the initial and the final 300 s were eliminated, using 3,600 s (3,600 points collected). Data were generated in decimal geographical coordinates in Datum WGS 84 and converted to the metric UTM system using software QGIS 2.18 (OSGEO, version 2.18.26).

For the static test, North-South (ENS) (Equation 1) and East-West (EEW) (Equation 2) errors were calculated, resulting in the radial error (E) (Equation 3), generated by the receiver in relation to the reference point, as described by Machado et al. (2010)Machado, T.M.; Molin, J.P.; Povh, F.P.; Salvi, J.V. 2010. Methodology for performance evaluation of GPS receiver for agricultural use in kinematic condition. Engenharia Agrícola 30: 121-129 (in Portuguese, with abstract in English)., using MS Excel.

(1) $E N S = | Y r e f − Y r e a l |$
(2) $E E W = | X r e f − X r e a l |$
(3) $E = E n s 2 + E e w 2$

where: (Xref, Yref) is the reference point coordinate (m); (Xreal, Yreal) is the coordinate of the GNSS receiver (m).

An algorithm based on Java language was developed for analyzing data obtained by the dynamic test, using software NetBeans IDE (version 8.1) to create the reference line and calculate the path error. Both were expressed as an error with 95 % probability (E95) and for each receiver, a root mean squared (RMS) error was calculated (Equation 4) (Pérez-Ruiz and Upadhyaya, 2012Pérez-Ruiz, M.; Upadhyaya, K.S. 2012. GNSS in Precision Agricultural Operations. IntechOpen, Rijeka, Croatia.).

(4) $R M S t = 1 N t ∑ i = 1 N t E 2$

where: RMSt is the RMS error for the tth receiver; Nt is the total number of measurement points for the tth receiver; E is the error of the point ith to the tth receiver.

For the statistical analysis, errors were calculated for each point and the SAS general linear models procedure was used to test for significant differences between receivers using the ANOVA. The Duncan test was used to compare the performance of the different receivers for their E95 mean accuracy, at 5 % error probability level. The probability transformation suggested by Box and Cox (1964)Box, G.E.P.; Cox, D.R. 1964. An analysis of transformations. Journal of the Royal Society 26: 211-252. was applied in all distributions, transforming E values to stabilize or reduce variability between replications. To compare the performance of different receivers with SBAS and RTK signal correction, we did not compare replications with a loss of signal correction. The number of replications may increase accuracy (Zhang et al., 2018Zhang, L.; Kong, Y.; Guo, Y.; Yan, J.; Wang, Z. 2018. Survey on network flow watermarking: model, interferences, applications, technologies, and security. IET Communications 12: 1639-1648.) and improve the capacity of a statistical test to detect smaller differences between estimates of the means of treatments (Danilogorskaya et al., 2017Danilogorskaya, E.A.; Zernov, N.N.; Gherm, V.E.; Strangeways, H.J. 2017. On the determination of the effect of horizontal ionospheric gradients on ranging errors in GNSS positioning. Journal of Geodesy 91: 503-517.).

Sensitivity of a statistical test is largely influenced by sample size. The test power was calculated to verify if a supposed difference between accuracies of receivers was genuine or subject to sample error for the static and dynamic tests. The test power corresponds to 1-β, the probability of rejecting the null hypothesis when it is false, and indicates the correct decision probability based on the alternative hypothesis. It is usually interpreted as the chance to detect a real difference between means or proportions. The simulations were made for various combinations of replication numbers with the nominal level of significance α equals to 5 %, admitting a difference between accuracies of receivers equals to the standard error of the mean.

# Results and Discussion

In the static test, the error dispersion for C/A receivers was relatively larger than for other receivers (Figure 5). The L1/L2 receivers without differential correction (ID 6) resulted in a low error dispersion. However, there was a shift of the data towards the south of the reference, generating greater ENS that decreased the accuracy of the receiver.

Figure 5
Dispersion of North-South (ENS) and East-West (EEW) errors in the static test. Cellular with GPS and GLONASS (ID 1); GNSS with signal C/A with only GPS (ID 2) and with GPS+GLONASS (ID 3, ID 4); GNSS L1 with GPS+GLONASS (ID 5); GNSS L1/L2 without correction (ID 6), with geo-satellite RTXSA connection (ID 7), geo-satellite AORW connection (ID 8 and ID 9) and Real-Time Kinematic correction (ID 10).

Receiver ID 5 also showed greater dispersion of errors in the 2nd and 3rd quadrants. Receivers with SBAS differential correction have similar dispersion of ENS and EWE. A large number of observations in the 2nd quadrant of receiver ID 9 may be attributed to problems with the SBAS signal correction, observed during the tests. Receiver ID 10, with RTK differential correction, displayed less dispersion of errors when compared to receivers with SBAS.

The L1/L2 receivers without differential correction showed low variation amplitude of E during data collection, within its replications (Figure 6). However, there was great variation in the temporal stability of the error. In the 3rd replication of receiver ID 7, with SBAS differential correction, there was a visible increase of E. In contrast, in the 4th replication, E started with a high value and decreased over time. Lim et al. (2019)Lim, C.; Yoon, H.; Cho, A.; Yoo, C.S.; Park, B. 2019. Dynamic performance evaluation of various GNSS receivers and positioning modes with only one flight test. Electronics 8: 1518. described it as a factor associated with the need of providing the same signal environment to all receivers.

Figure 6
Variation of the parallelism error (E) as a function of the time of data collection for all GNSS receivers with L1/L2 signal in the static test. GNSS with geo-satellite RTXSA connection (ID 7), geo-satellite AORW connection (ID 8 and ID 9), and Real-Time Kinematic correction (ID 10).

Figure 7 shows the dispersion of errors during the dynamic tests. There was a small error dispersion for the receiver with RTK differential correction in relation to the other receivers. In general, the dynamic test conditions resulted in greater dispersion of errors when compared to the static test. The L1/L2 receiver without differential correction (ID 6) showed high error variation during the time of data collection within replications, which did not occur during the static test, with error values up to 1.45 m.

Figure 7
Variation of the parallelism error (E) as a function of the time of data collection for all GNSS with L1/L2 signal in the dynamic test. GNSS with geo-satellite RTXSA connection (ID 7), geo-satellite AORW connection (ID 8 and ID 9), and Real-Time Kinematic correction (ID 10).

In addition, a variation of the position errors of receivers with SBAS leads to a greater differential correction in the dynamic tests. Similar to the static tests, there was a lack of correction signal for receivers with SBAS signal correction, which can be visualized in replication number 4 of receiver ID 8 and replication number 1 of receiver ID 9 when there is more variation of positioning errors. Although these two replications presented a large amplitude of errors, the other replications presented variation below 0.250 m. However, when compared to the static test, there is greater error variation in the dynamic tests.

Table 2
Error-values calculated for static and dynamic tests.

According to Table 2, the coefficient of variation (CV) between the replication for Group A did not demonstrate any standard. Some CV values were higher (> 50 %) in the static test in relation to the dynamic test (19 %). For Group B, except for receiver ID 7, there was an increase in error variation between replications in the dynamic test. It is considered that the lower the CV estimates, the greater the precision of the experiment and vice versa, and the greater the experimental precision.

It was expected an L1/L2 receiver without differential correction to be more accurate than C/A code and L1 band receivers. This is due to better reception capacity of the signal, as it can receive weak signals from satellites without self-interference, and it has a better structure for positioning data recovery. Considering that ionospheric delay is a relevant error source in multi-GNSS positioning, Su et al. (2019)Su, K.; Jin, S.; Hoque, M.M. 2019. Evaluation of ionospheric delay effects on multi-GNSS positioning performance. Remote Sensing 11: 171. and Wang et al. (2019)Wang, A.; Chen, J.; Zhang, Y.; Meng, L.; Wang, J. 2019. Performance of selected ionospheric models in multi-global navigation satellite system single-frequency positioning over China. Remote Sensing 11: 2070. verified that the combination of multi-GNSS could effectively improve the positioning performance.

Receivers of Group A, except receiver ID 6, had greater accuracy in the dynamic test. It was expected greater error variation in the dynamic test in relation to the static test, due to the continuous changes in the configuration of the satellite constellation and response time of the receiver, causing greater accuracy variation of GNSS positioning (ASABE, 2008American Society of Agricultural and Biological Engineers [ASABE]. 2008. X605 Satellite-Based Auto-Guidance Systems Testing During Straight and Level Travel. ASABE, St. Joseph, MI, USA.). When performing the static and dynamic tests in different moments, these factors also change, implying greater dataset variability.

For GNSS L1/L2 with differential correction, there was greater variability of E95 between replications in the dynamic test in relation to the static test (Table 2). In a static test, GNSS receivers with differential correction had an E95 lower than 0.100 m and E95 of receiver ID 7 was significantly different from GNSS receiver ID 10, with RTK differential correction. E95 of all receivers with correction used in the static test was closer to the mean errors, indicating low variation in positioning.

GNSS receivers with SBAS resulted in E95 greater than receiver ID 10, with RTK differential correction in the dynamic test and SBAS receivers had higher E95 in the dynamic test when compared to the static test. Despite greater variation in the positioning of GNSS receivers in the dynamic test, E95 of receiver ID 10 was higher in the dynamic test.

Comparing L1/L2 GNSS receivers, E95 of the receiver without differential correction is greater than the other receivers with differential correction for both tests (static and dynamic). This means that the use of an L1/L2 GNSS receiver without differential correction could generate an E95 of 0.843 m higher than a receiver with SBAS signal correction in a static condition, and 0.976 m higher than a receiver with SBAS differential correction in the dynamic condition. Therefore, a lack of differential correction signals could provide a considerable increase in GNSS signal error during agricultural operations.

The effect of the replication number indicated that data from receivers with differential correction had a power of 81 % in the static and of 44 % in the dynamic test (Figure 8), which means that the tests were performed with a small number of replications. The low value of the test power for the dynamic test in Group B could be explained by a greater variety of error values in the dynamic than in the static test. Considering that higher CV implies less experimental accuracy and, consequently, more replications are required to represent GNSS positioning error.

Figure 8
Power of the test versus the number of replications on static and dynamic tests for receivers from Group A and Group B.

To test receivers with a differential correction to have at least a power of the 80 % (β = 20 %), we carried out seven replications for the static test and 13 replications for the dynamic condition. For receivers with no differential correction (Group A), the test power was 46 % for the static test and 96 % for the dynamic test conditions, showing that nine replications are necessary for the static and five replications for the dynamic test conditions to achieve a minimum test power of 80 %. The results demonstrate that the number of replications suggested in the power test to compare different GNSS with differential correction complies with the number of replications indicated by ION STD 101 (1997)Institute of Navigation [ION]. 1997. STD 101: Recommended Test Procedures for GPS Receivers, Revision C, January 27, 1997. Institute of Navigation, Alexandria, VA, USA..

# Conclusion

This study developed a methodology that allowed the behavior analysis of GNSS receivers of different accuracy levels under static and dynamic conditions. The test under dynamic conditions measuring error perpendicular to the path simulates agricultural demands, for example, the use of GNSS for auto-guidance on agricultural operations. The signal of receivers presented greater variability of positioning errors in the dynamic than in the static condition, consequently, we observed the need to increase the number of replications to improve the power of the statistical test. A minimum of five replications is necessary for the dynamic tests to C/A receivers and 13 replications for L1/L2 receivers with differential correction signals. The static tests require a minimum of nine replications for C/A receivers and five replications for L1/L2 receivers with differential correction signals.

# Acknowledgments

To the Brazilian National Council for Scientific and Technological Development (CNPq) for the Ph.D. scholarship granted to the first author. To the Coordination for the Improvement of Higher Level Personnel (CAPES) - Finance Code 001, for the Ph.D. scholarship granted to the second author.

# References

• American Society of Agricultural and Biological Engineers [ASABE]. 2008. X605 Satellite-Based Auto-Guidance Systems Testing During Straight and Level Travel. ASABE, St. Joseph, MI, USA.
• Bae, T.S.; Kim, M. 2018. Performance analysis of Network-RTK techniques for drone navigation considering ionospheric conditions. Journal of Sensors 1: 1-8.
• Banville, S.; Collins, P.; Lahaye, F. 2018. Model comparison for GLONASS RTK with low-cost receivers. GPS Solutions 22: 1-12.
• Bažec, M.; Dimc, F.; Pavlovčič-Prešeren, P. 2020. Evaluating the vulnerability of several geodetic GNSS receivers under chirp signal L1/E1 jamming. Sensors 20: 814.
• Box, G.E.P.; Cox, D.R. 1964. An analysis of transformations. Journal of the Royal Society 26: 211-252.
• Carballido, J.; Perez-Ruiz, M.; Emmi, L.; Aguera, J. 2014. Comparison of positional accuracy between RTK and RTX GNSS based on the autonomous agricultural vehicles under field conditions. Applied Engineering in Agriculture 30: 361-366.
• Danilogorskaya, E.A.; Zernov, N.N.; Gherm, V.E.; Strangeways, H.J. 2017. On the determination of the effect of horizontal ionospheric gradients on ranging errors in GNSS positioning. Journal of Geodesy 91: 503-517.
• Institute of Navigation [ION]. 1997. STD 101: Recommended Test Procedures for GPS Receivers, Revision C, January 27, 1997. Institute of Navigation, Alexandria, VA, USA.
• International Organization for Standardization [ISO]. 2012. ISO12188-2: Tractors and Machinery for Agriculture and Forestry – Test Procedures for Positioning and Guidance Systems in Agriculture. Part 2. Testing of Satellite-Based Auto-Guidance Systems During Straight and Level Travel. ISO, Geneve, Switzerland.
• Kabir, M.S.N.; Song, M.Z.; Sung, N.S.; Chung, S.O.; Kim, Y.J.; Noguchi, N.; Hong, S.J. 2016. Performance comparison of single and multi-GNSS receivers under agricultural fields in Korea. Engineering in Agriculture, Environment and Food 9: 27-35.
• Kim, J.Y.; Kim, H.J.; Shim, S.B.; Park, S.H.; Kim, J.H.; Kim, Y.J. 2016. Development of automated guidance tracking sensor system based on laser distance sensors. Journal of Biosystems Engineering 41: 319-327.
• Li, M.; Yuan, Y.B.; Wang, N.B.; Liu, T.; Chen, Y.C. 2018. Estimation and analysis of the short-term variations of multi-GNSS receiver differential code biases using global ionosphere maps. Journal of Geodesy 92: 889-903.
• Lim, C.; Yoon, H.; Cho, A.; Yoo, C.S.; Park, B. 2019. Dynamic performance evaluation of various GNSS receivers and positioning modes with only one flight test. Electronics 8: 1518.
• Machado, T.M.; Molin, J.P.; Povh, F.P.; Salvi, J.V. 2010. Methodology for performance evaluation of GPS receiver for agricultural use in kinematic condition. Engenharia Agrícola 30: 121-129 (in Portuguese, with abstract in English).
• Machado, T.M.; Molin, J.P. 2011. Static and cinematic tests of GPS receivers. Revista Brasileira de Engenharia Agrícola e Ambiental 15: 981-988 (in Portuguese, with abstract in English).
• Pérez-Ruiz, M.; Upadhyaya, K.S. 2012. GNSS in Precision Agricultural Operations. IntechOpen, Rijeka, Croatia.
• Rovira-Más, F.; Chatterjee, I.; Sáiz-Rubio, V. 2015. The role of GNSS in the navigation strategies of cost-effective agricultural robots. Computers and Electronics in Agriculture 112: 172-183.
• Sama, M.P.; Stombaugh, T.S. 2014. Performance evaluation of a tracking total station as a position reference for dynamic GNSS accuracy testing. Applied Engineering in Agriculture 30: 557-563.
• Santos, A.F.; Silva, R.P.; Tavares, T.O.; Ormond, A.T.S.; Rosalen, D.L.; Assis, L.C. 2017. Parallelism error in peanut sowing operation with auto-steer guidance. Revista Brasileira de Engenharia Agrícola e Ambiental 21: 731-736.
• Silva, E.G.S.; Marques, H.A. 2016. PPP with GPS/GLONASS data Integration: Fundamentals Involved and Accuracy Assessment. Boletim de Ciências Geodésicas 22: 217-231 (in Portuguese, with abstract in English).
• Souza, S.V.R.; Machado, A.M.L. 2016. Systematic Errors Calibration for MEMS Accelerometers. Boletim de Ciências Geodésicas 22: 835-850 (in Portuguese, with abstract in English).
• Su, K.; Jin, S.; Hoque, M.M. 2019. Evaluation of ionospheric delay effects on multi-GNSS positioning performance. Remote Sensing 11: 171.
• Wang, A.; Chen, J.; Zhang, Y.; Meng, L.; Wang, J. 2019. Performance of selected ionospheric models in multi-global navigation satellite system single-frequency positioning over China. Remote Sensing 11: 2070.
• Ye, S.; Zhao, L.; Song, J.; Chen, D.; Jiang, W. 2018. Analysis of estimated satellite clock biases and their effects on precise point positioning. GPS Solutions 22: 1-14.
• Zhang, L.; Kong, Y.; Guo, Y.; Yan, J.; Wang, Z. 2018. Survey on network flow watermarking: model, interferences, applications, technologies, and security. IET Communications 12: 1639-1648.

# Publication Dates

• Publication in this collection
24 Aug 2020
• Date of issue
2021

# History

17 Sept 2019
• Accepted
12 Apr 2020
Escola Superior de Agricultura "Luiz de Queiroz" USP/ESALQ - Scientia Agricola, Av. Pádua Dias, 11, 13418-900 Piracicaba SP Brazil, Phone: +55 19 3429-4401 / 3429-4486 - Piracicaba - SP - Brazil
E-mail: scientia@usp.br